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ABSTRACT 

A procedure for the numerical solution of the complete, iso- 
thermal elastohydrodynamic lubrication problem for point contacts 
is given. This procedure calls for the simultaneous solution of che 
elasticity and Reynolds equations. By using this theory the influ- 
ence of the ellipticity parameter and the dimensionless speed, load, 
and material parameters on the minimum and central film thicknesses 
was investigated. Thirty-four different cases were used in obtaining 
fully flooded minimum- and central-film-thickness formulas. Lubri- 
cant starvation was also studied. From the results it was possible 
to express the minimum film thickness for a starved condition in 
terms of the minimum film thickness for a fully flooded condition, 
the speed parameter, and the inlet distance. Fifteen additional 
cases plus three fully flooded cases were used in obtaining this 
formula. Contour plots of pressure and film thickness in and around 
the contact have been presented for both fully flooded and starved 


lubrication conditions. 
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CHAPTER I 


INTRODUCTION 

1-1 Statement of Problem 

In many contacts between machine elements, forces are trans- 
mitted through thin, but continuous, fluid films. One of the basic 
problems is to accurately describe the fluid film thickness between 
these machine elements. The provision of an adequate fluid film 
thickness will reduce wear and increase fatigue life and therefore 
avoid early damage of the machine elements. These fluid films as re- 
lated to hydrodynamic lubrication in journal and thrust bearings have 
been well understood for some time, and experimental work confirms 
the theory. In the early 1900's it was recognized by Martin (1916) 
that many loaded contacts of low geometrical conformity, commonly re- 
ferred to as nonconforming contacts, such as gears and rolling- 
contact bearings, behaved as though they were hydrodynamically lu- 
bricated. As opposed to the journal and thrust bearing counterpart, 
the original hydrodynamic lubrication theory of gears and rolling- 
contact bearings differed substantially from experimental findings. 
Only in recent years has the consideration of elastic deformation of 
contacts been coupled to hydrodynamics to yield a closer agreement 
of theory with experiments. 

Elastohydrodynamic lubrication (EHL) then deals with the lubri- 
cation of elastic contacts. The analysis requires the simultaneous 
solution of the elasticity and Reynolds equations. The EHL theory 
differs from conventional hydrodynamic theory in the following way: 

(1) In defining the film thickness in EHL theory, elastic de- 


formation of the contact is considered. 

(2) The viscosity is no longer independent of pressure, as is 
assumed in conventional hydrodynamic theory. 

(3) Hydrodynamic lubrication is characterized by surfaces that 
are conforming, but elastohydrodynamic lubrication is usually char- 
acterized by nonconforming surfaces. 

Because of this last point the load in hydrodynamic lubrication is 

usually carried over a relatively large area. Also, typical maximum 

pressures for elastohydrodynamic lubrication are of the order of 

1. 4xl0 9 N/m 2 (200 000 lb/in. 2 ), and the usual hydrodynamic pressures 
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generated in journal bearings are of the order of only 7x10 N/m 

2 

(1000 lb/in. ). 

1. 2 Historical Developments 

When two solids are in contact under zero load condition, one 

of two types of contact is experienced: 

(1) Point contact, in which two solids touch at a single point, 

as in ball bearings 

(2) Line contact, in which two solids touch along a straight or 

curved line, as in roller bearings 
After a load is applied, the point expands to an ellipse and the 
line to a rectangle. Although we are concerned with loaded contacts, 
it is convenient to distinguish between these situations by referring 
to them as being either point or line contact. Lubricant removal o 
the side of the contact is ignored in the line-contact problem and 
the problem becomes two dimensional. This is a considerable simpli- 
fication from the point-contact, three-dimensional problem. 

1.2.1 Line Contact 

One of the earliest solutions of the lubrication of a line con- 
tact was presented in 1916 by Martin. By assuming rigid solids and 
an incompressible constant-viscosity lubricant he was able to deter- 
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mine a lubricating film thickness. Figure 1.1(a) shows typical pres 
sure and film thickness curves from Martin's solution as well as his 
formula for minimum film thickness. Martin’s formula greatly under- 
estimates the film thickness; however, it was a useful beginning to 
the theoretical study of elastohydrodynamic lubrication of line 
contacts. 

Some 30 years transpired before any significant accomplishments 
were made in solving the EHL line-contact problem. Grubin (1949) 
obtained the first satisfactory solution to this problem by taking 
account of elastic distortion and viscosity-piessure effects. In 
Grubin s analysis it was assumed that the shape of the elastically 
deformed solids in a highly loaded lubricated contact is the same as 
the shape produced in a dry Hertzian contact. A pressure distribu- 
tion in a Hertzian contact is shown in figure 1.1(b). This stipula- 
tion (assuming the shape of the elastically deformed solids is the 
same as is produced in a Hertzian contact) facilitated the solution 
of the Reynolds equation in the inlet region of the contact and en- 
abled the separation of the solids in the central region of the con- 
tact to be determined with commendable accuracy. Grubin’s approach 
produced an excellent account of the physical mechanism of the lubri- 
cation process in highly loaded EHL line contacts, and it marked a 
very important development in the history of elastohydrodynamic 
lubrication. 

Dowson and Higginson (1961) produced an empirical formula for 
isothermal EHL line contacts. This formula shows the effect of 
speed, load, and material properties on minimum film thickness and 
is based on their earlier theoretical solutions (1959). Figure 
1.1(c) gives the pressure and film thickness for an elastohydrody- 
namic lubricated line contact. Also shown in figure 1.1(c) is the 
KHL line-contact minimum film thickness obtained from the results of 
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Dowson and Higginson (1961). In the Dowson-Higginson theory (1959) 
for EHL line contacts, a new approach of introducing a solution of 
the inverse hydrodynamic lubrication problem was presented. Nor- 
mally a solution oi the Reynolds equation calls for the determina- 
tion of a pressure distribution corresponding to a given film thick- 
ness. In the inverse problem the film shape responsible for the 
generation of a given pressure distribution is determined. In the 
procedure adopted by Dowson and Higginson (1959) the computed film 
shape was compared with the shape of the elastically deformed solids, 
and the pressure curve was then modified to improve the agreement 
between the uwo shapes. By using this approach they were able to 
obtain satisfactory solutions of the elastic and hydrodynamic equa- 
tions after a small number of numerical iterations. 

The Dowson-Higginson minimum-film-thickness formula for EHL 
line contacts agrees quite well with experimental observations. In 
particular, Sibley and Orcutt's (1961) X-ray method, Christensen's 
(1964) slide displacement method, and Dyson, Naylor, and Wilson's 
(1966) capacity method all give good agreement with the Dowson- 
Higginson formula. 

1.2.2 Point Contact 

Most of the work on elastohydrodynamic lubrication has dealt 
with line contacts. Furthermore, the majority of the work done on 
the EHL point-contact problem has been experimental. Two good ex- 
amples of this experimental work are Cameron and Gohar's (1966) ob- 
servation of film thickness between a steel sphere and a glass plate 
using interference rings and Archard and Kirk's (1961) observation 
of film thickness between two crossed cylinders with the same diam- 
eter makin ; .in angle of 90° in a simulated point contact. 


The first step toward a theoretical solution of the EHL point- 
contact problem was presented by Archard and Cowking (1966). They 
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adopted an approach similar to that used by Crubin (1949) for line- 
contact conditions. The Hertzian contact zone was assumed to form 
a parallel film region, and the generation of high pressure in the 
approaches to the Hertzian zone was considered. 

Cheng (1970) also used a Grubin type of approach in determin- 
ing a film thickness formula for the EHL of point contacts. He 
evaluated the deformation by using the Hertz equation and then ap- 
plying the Reynolds equation to this geometry. 

Recently an interesting numerical solution of the EHL point- 
contact problem for a sphere near a plane was put forth by Ranger, 
et al. (1975). This solution is presented in dimensional terms, which 
thereby limits its general usage. A puzzling feature of the Ranger, 
et al. work is the fact that his resulting equation for the minimum 
film thickness has a positive load exponent, which contradicts exper- 
iments (e.g. .Cameron and Gohar r s(1966) and Archard and Kirk (1961)). 

1. 3 Approach to the Problem 

In the literature the EHL line-contact problem is completely 
solved, with theory and experiment agreeing well with each other. 

Tor the EHL point-contact problem quite a bit of experimental work 
has been done, but the theoretical solution to the complete isother- 
mal e las to hydrodynamic lubrication problem for point contacts has 
not emerged. The reason for this is the extreme difficulty of the 
numerical coupling of the elasticity and Reynolds equations and the 
Vc.st increase in computing involved in a transfer from line- to 
point-contact conditions. The work presented in this dissertation 
is an attempt to solve this problem. A brief description of the ap- 
proach to the problem follows. 

The radii of curvature of the contacting solids ^re used to de- 
fine the ratio of the semimajor and semiminor axes of the contact 
ellipse. From this ratio and the normal applied load, the semimajor 
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and semiminor axes of the contact ellipse are determined. In the 
elasticity analysis the computing zone is divided into equal rec- 
tangular areas, and a uniform pressure is applied over each area. 
Elasticity studies are performed to determine how finely the semi- 
major and semiminor axes need to be divided to achieve a given ac- 
curacy and how far from the center of the contact deformation be- 
comes insignificant compared with the separation of the solids. 
These equations were investigated for light and heavy loads and 
for geometries ranging from a ball on a plate to a line-contact 
configuration. The answers to these questions determined what the 
computing zone should be in and around the contact ellipse. 

In the numerical analysis of the Reynolds equation a "phi" 

(<£) substitution is used to aid the relaxation process, where <j> 
is equal to the pressure times the film thickness to the 3/2 power. 
The pressure-viscosity analysis of Roelands (1966) is used. The 
numerical coupling of the elasticity and Reynolds equations results 
in a converged solution for the pressure profile. This pressure 
profile Is then integrated over the computing zone to give the value 
of the corresponding normal applied load. This load is then com- 
pared with the input load, and corrections are made to the film 
thickness until these two loads are in agreement. 

The most important practical aspect of the EHL point-contact 
theory is the determination of the minimum film thickness within the 
contact. By using the variables resulting from the theory, the di- 
mensionless minimum film thickness can be written as a function of 
the dimensionless load, speed, material, and ellipticity parameters. 

The influence of these parameters on minimum film thickness was 
investigated. The ellipticity parameter was varied from 1 (a ball 
on a plate) to 8 (a cont iguration approaching a line contact). The 
dimensionless speed parameter was varied over a range of nearly two 


orders of magnitude. The dimensionless load parameter was varied 
over a range of one order of magnitude. Solid materials of bronze, 
steel, and silicon nitride and lubricants of paraffinic and naph- 
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thenic mineral oils were varied in obtaining the exponent in the 
dimensionless material parameter. Thirty-four different cases were 
used in obtaining the fully flooded minimum-f ilm-thickness formula. 

In addition to the minimum-! ilm-thickness formula a central-film- 
thickness formula was developed. Contour plots are also shown that 
indicate in detail the pressure spike and the two side lobes in which 
the minimum film thickness occurs. These theoretical solutions for 
film thickness have all the essential features of the previously re- 
ported experimental observations based upon optical interferometry. 

The effect of lubricant starvation was also investigated. This 
study of lubricant starvation was achieved simply by reducing the 
inlet distance, which is the distance from the center of the contact 
to the edge of the computing area. A fully flooded condition exists 
when the dimensionless inlet distance ceases to influence in any 
significant way the minimum film thickness. Starting from a fully 
flooded condition and decreasing the inlet distance, the value at 
which the minimum film thickness first starts to change is called 
the fully flooded - starved boundary. Simple expressions for the 
fully flooded - starved boundary were obtained as a function of the 
fully flooded central or minimum film thickness. Simple expressions 
defining the central and minimum film thicknesses for a lubricant 
starvation condition were also obtained. Fifteen different cases 
from those presented in the fully flooded results were used in ob- 
taining these formulas. Furthermore, the effects of lubricant 
starvation are clearly shown in contour plots of press. e and film 
thickness . 
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CHAPTER 2 


GEOMETRY OF CONTACTING SOLIDS 
2.1 Curvature Sum and Difference 

Two solids having different radii of curvature in a pair of 
principal planes (x and y) passing through the contact between ^he 
solids make contact at a single point under the condition of no ap 
piied load. Such a condition is called "point contact" and is 
shown in figure 2.1, where the radii of curvature are ienoted by 
r's. In the analysis that follows it was assumed that for convex 
surfaces as shown in figure 2.1 the curvature is positive but that 
for concave surfaces the curvature is negative. 

The curvature sum and difference are defined as 



where 


1 

R 

x 




1 

R 

y 




(2.3) 


(2.4) 


2-2 Geometric Separation o f Ellipsoidal Solids 

Figure 2.2 shows how the geometric separation between two ellip- 
soidal solids can be made equivalent to that between a single ellip- 
soidal solid near a plane. The geometric requirement is that for any 
values of x and y in figure 2.2(a) the geometric separation must 


be equivalent to the separation at the same 


x ami y shown in f ig- 


ure 2.2(b). From figure 2.2(a) the mathematical expression for the 
separation of the two ellipsoidal solids can be written as 
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s = s. + s_ + S A + s„ 

Ax Bx Ay By 


(2.5) 


From figure 2.2(b) the mathematical expression for the separation of 
a single ellipsoidal solid near a plane can be written as 

S = S x + S y (2.6) 

Therefore, for the two expressions to be equivalent, the following 
must be true: 


s = S. + S D 
x Ax Bx 


s = S A + S_ 
y Ay By 


From figure 2.2(a-l) the following can be written: 


r 2 = x 2 + (r - S ) 2 
Ax Ax Ax 


(2.7) 

( 2 . 8 ) 

(2.9) 


or 


x = S. (2r. - S. ) 

Ax Ax Ax 


( 2 . 10 ) 


But for the problem being considered, 2r >> S ; so that equation 

AX AX 

(2.10) can be rewritten as 


C '"b 11 

Ax ^ 2r 


( 2 . 11 ) 


Ax 


This is the well-known parabolic approximation to the circular sec- 
tion of the solid. Similarly, by making use of figure 2.2 the fol- 
lowing can be written: 


% 


Bx 2r 


S, 


Bx 

2 


'Ay ^ 2r 


Ay 
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s. 


'By * 2r 


By 


where 2r„ >> S_ 

Bx Bx 


where 2r A >> S k 
Ay Ay 


where 2r_ >> 

By By 


( 2 . 12 ) 


(2.13) 


(2.14) 
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Substituting equations (2.11) 
gives 


where 2R >> s 
x x 


where 2R >> S 

y y 

to (2.16) into equations 


(2.15) 


(2.16) 
(2.7) and (2.8) 



But equations (2.17) and (2.18) are exactly equations (2.3) and (2.4), 
respectively. Therefore the equivalency shown in figure 2.2 is sat- 
isfied. Henceforth the geometry of an ellipsoidal solid near a plane 
as shown in figure 2.2(b) will be used. 

2* 3 Ellipticitv Parameter 

When the ellipsoidal solid just touches the plane shown in fig- 
ure 2.2(b), contact is made at a single point when no load is ap- 
plied. When a normal load is applied to the ellipsoidal solid, the 
point expands to an ellipse with a as the semimajor axis and b 
as the semiminor axis. It is assumed that the plane remains rigid 
while the equivalent elastic ellipsoid is pressed against it. The 
normal applied load lies along the axis that passes through the cen- 
ter of the solid and through the point of contact and is perpendicu- 
lar to the plane shown in figure 2.2(b). For the special case where 

R x - R y> the resu lting contact is a circle rather than an ellipse. 

The ellipticity parameter (k) is defined as 


w - £ 

k b (2.19) 

where 

a semimajor axis of contact ellipse 

b semiminor axis of contact ellipse 
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Harris (1966) has shown that the ellipticity parameter can be written 
to relate the curvature difference and tne elliptic integrals of the 
first and second kind as 


where 






J(k) = m. ± r j 

Ju; \ <m - r) 


/ 

/ 


1 - /l “j singer 


- 1/2 


dep 


!-(■• 


1/2 


1 \ . 2 

-z- sin (p 

k / 


d<p 


( 2 . 20 ) 


( 2 . 21 ) 


( 2 . 22 ) 


A one-point iteration method that has been used successfully in the 
past by Hamrock and Anderson (1973) was used, where 

k„ +1 - J(k n > (2.23) 

The iteration process is continued until k , , differs from k 

n+1 n 

by less than 1x10 Note that the ellipticity parameter is a func- 
tion of the radii of curvature of the contacting solids only 

k = f (r ,r_ ,r. ,r_ ) 

Ax’ Bx Ay By 

When the ellipticity parameter (k) , the normal applied load 
(F) , Poisson’s ratio (v) , and modulus of elasticity (E) of the con- 
tacting solids are known, the semimajor axis of the contact ellipse 
can be written as 


a = 


/ 2 \l/3 

(6k F#R\ 

\ ” e ' ) 


(2.24) 


where 


F' = 


i 2 , 2 

1 ~ V A 1 ~ V fl 

e a e b 


(2.25) 


By making use of equation (2.19) the semiminor axis of the contact 
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ellipse can be written as 

b = f ( 2 . 26 ) 

k 

Therefore, from the geometry of the contacting solids the ellipticity 
parameter, as well as the semimajor and seraiminor axes, has been de- 
fined. These parameters form the foundation of the EHL point-contact 
analysis. 
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CHAPTER 3 

NUMERICAL EVALUATION OF THE ELASTIC DEFORMATION OF 
SOLIDS SUBJECTED TO A HERTZIAN CONTACT STRESS 

Elastohydrodynamic lubrication is defined as the study of 
situations in which elastic deformation of the surrounding solids 
plays a significant role in the hydrodynamic lubrication process. 

This chapter is not concerned with the hydrodynamic lubrication 
process, but only with deformation due to the pressure of one 
elastic solid upon another. 

Dowson (1965) distinguishes between two modes of deformation 
that may exist in machine elements. In one mode, the contact 
geometry may be affected by overall distortion of the elastic ma- 
chine element resulting from applied loads, as shown in figure 
3.1(a). In the other, the normal stress distribution in the vi- 
cinity of the contact zone may produce local elastic deformations 
that are significant when compared with the lubricant film thickness, 
as shown in figure 3.1(b). This is the mode of deformation with 
which this present investigation is concerned. The important dis- 
tinction is that the first form of deformation is relatively insen- 
sitive to the distribution and magnitude of the stresses in the con- 
tact zone, whereas the second mode of deformation is intimately 
linked to the local stress conditions. 

The correct evaluation of elastic deformation on the surface of 
a solid depends upon an adequate representation of the applied nor- 
mal pressures. The simplest procedure is to divide the actual pres- 
sure distribution into rectangular blocks of uniform pressure and to 

; , i r)[ i if ;i i i i s •. n lit e, 

UKlulNAi. i Mu j i;> 1\X)R 
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permit each rectangle to be of such small dimensions that adequate 
predictions of elastic displacements ensue. More complex represen- 
tations of two-dimensional pressure distributions within each rec- 
tangle would generally permit larger rectangles to be used, but the 
additional complexity of the expressions and added computation time 
make it desirable to exploit the simpler representation to the full- 
est. This chapter is devoted to a study of the adequacy of repre- 
senting the applied normal pressures by rectangular blocks of uni- 
form pressure. 

The deformation analysis is developed in general form since it 
will be used in later chapters to calculate elastic deformation in 
the elastohydrodynamic lubrication of point contacts. To evaluate 
the influence of mesh or clock size upon accuracy, the numerical in- 
vestigation of this chapter con ; ders only Hertzian contact stress 
distributions. 

The deformation analysis it. elf assumed that the contact zone 
can be divided into rectangular areas and that the pressure is uni- 
form within each rectangular area. Once the elastic deformation 
had been formulated, investigations were performed to answer the 
following questions: 

(1) How fine do the semimajor and semiminor axes need to be 
divided to achieve a given accuracy in deformation prediction? 

(2) How far from the center of the contact does deformation be- 
come insignificant compared with the separation of solids? 

These questions were investigated for both light and heavy applied 
loads and for both equal spheres in contact and a contact that is 
common to the outer race of a ball bearing. 

3. 1 Elastic Deformation Analysis 

In chapter 2 the general geometry of two ellipsoidal solids in 
elastic contact was described. In the subsequent analysis it will 
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be convenient to consider the deformation of an equivalent elastic 
half-space subjected to a Hertzian pressure distribution over the el- 
lipse of semimajor and semiminor axes, a and b, as previously de- 
fined. The resulting elastic deformation can be considered to be 
equivalent to the total deformation of two elastic ellipsoids having 
elastic constants E^,v^ and E g> v g> respectively, if the half- 
space is allocated the equivalent elastic parameter (E') defined by 
equation (2.25). 

Once the semimajor and semiminor axes of the contact ellipse 
have been defined, the elastic deformation that occurs inside and 
outside the contact zone can be evaluated. Figure 3.2 shows a rec- 
tangular area of uniform pressure with the coordinate system to be 
used. From Timoshenko and Goodier (1951) the elastic deformation at 
a point (X, Y) of a semi-infinite solid subjected to a pressure (p) 
at the point (X^,Y^) can be written as 


dw = 


2p dX dY 
ttE'7 


The elastic deformation at a point (X,Y) due to the uniform pressure 
over the rectangular area 2a x 2b is thus 

>b 


w = 


2P 



dX 2 dY 1 


/- a 


/(Y - Y : ) 2 + (X - X x ) 2 


where 



Integrating the preceding equation gives 


w^ere 



7T 


(3.1) 
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D = (X + b) In 


(Y -f a) + 
(Y - a) + 


V (Y + a) 2 + (X + b) 2 

V(Y - I ) 2 + (X + b) 2 



(3.2) 


As a check on the validity of equation (3.1) the following two 
cases were evaluated: 

Case 1: For b = a and X = Y = 0, equation (3.1) reduces to 

W = ~ Pa ln(l + -/2) ( 3 . 3 ) 

Equation (3.3) represents the elastic deformation at the center of a 
square of uniform pressure. This equation is in agreement with that 
shown by Timoshenko and Goodier (1951). 

Case 2: For b = a and X = Y = a, equation (3.1) reduces to 

w = ® Pa ln(l + -/2) (3.4) 

Equation (3.4) represents the elastic deformation at the corner of a 
square of uniform pressure. This equation is also in agreement with 
that of Timoshenko and Goodier (1951). From equations (3.3) and (3.4) 
we find the corner deformation to be one-half the deformation at the 
center of a square block of pressure. 

Now the elastic deformation (w) in equation (3.1) represents 
the elastic deformation at a point (X,Y) due to a rectangular area 
2a x 2b of uniform pressure (p) . tf the contact ellipse is divided 
into a number of equal rectangular areas, the total deformation at a 
point (X,Y) due to the contributions of the various rectangular areas 
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of uniform pressure in the contact ellipse can be evaluated numer- 
ically. Figure 3.3 shows how the area inside and outside the contact 
ellipse may be divided into a number of equal rectangular areas. 

For purposes of illustration the contact was divided into a grid of 
6 x 6 rectangular areas. The effects of the fineness of this grid 
are discussed in section 3.4 of this chapter entitled Conditions In - 
v estigated . Figure 3.3 can be used to write the total elastic de- 
formation, caused by the rectangular areas of uniform pressure 

within the contact ellipse, at any point inside or outside the con- 
tact ellipse as 


w. n 
k, £ 



6 


E 

i=l , 2 , . . . 


P. .D 
i » J m,s 


where 


(3.5) 


hi = |lc - i | + l (3.6) 

8 = U - j I + 1 (3.7) 

Note that D ltl would be D in equation (3.2) evaluated at 
while D 2 ^ would be evaluated at X=2b\Y=4a. 

Equation (3.8) points out more explicitly the meaning of equa- 
tion (3.5). The elastic deformation at the center of the rectangu- 
lar area w^ (shown in fig. 3.3) caused by the pressure of the 
various rectangular areas in the contact ellipse can be written as 


“9.5 ’ V {P 1,1 D 9.5 + P 2,1 D 8,5 + • • • + 

+ P 1,2 D 9,4 + P 2,2 D 8,« + ’ ' ' + P 6,2 D 4,4 


+ P 1,6 D 9,2 + P 2.6°8,2 + 


+ P 6,6 D 4,2' 


(3.8) 


3 - 2 Hertzian Pressure DistrihnHnn 

In this chapter it is assumed that within the contact ellipse 
the pressure is described by the theory of Hertz (1882). Hertz 
(1882), using a purely elastomechanical process, was able to describe 
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the pressure distribution developed during contact of an ellipsoidal 
and a plane solid (as those described in fig. 2.2(b)) that are un- 
lubricated. Hertz found this pressure distribution to be semiellip- 
soidal, and for any arbitrary point within the contact ellipse it 
can be described by the following equation: 

v 2 


P(X,Y) = 


— 3F_ / / Y - a ' 

2irE'ab V 1 ' [ ~T~ y 


b 
b 


(3.9) 


Note that the coordinate system common to figure 3.3 is used in equa- 
tion (3.9). Inside the contact area the pressure was assumed to be 
described by equation (3.9); outside the contact area it was assumed 
to be zero. Therefore, for example, from figure 3.3, would be 

equivalent to the dimensionless pressure (P) from equation (3.9) 
evaluated at X = 5b and Y = 7a. 

3 ■ 3 Film Thi ckness 

By making use of equations (2.5) and (2.9) developed in chap- 
the exact geometric separation between the ellipsoidal solid 
and the plane shown in figure 2.2(b) for the coordinate system 
developed in figure 3.3 can be written as 

IZ 2 9 

(3.10) 


s ■ r * ■ 


(X - b) 2 + R - ^R 2 _ ( Y - a) 2 


y * y 

As discovered in chapter 2, if 2 R x » S x and 2R y » Sy , an ap- 
proximate expression for the separation between the solid and plane 
shown in figure 2.2(b) can be written as 


IX -..bl 2 (Y - a) 2 

2R + 2R 

X y 


(3.11) 


The degree to tfhlch the approximate equation (3.11) represents the 

exact separation of the ellipsoidal solid and plane Is determined by 
the following ratio: 


R i - 


( 3 . 12 ) 


The film thickness in an elastohydrodynamic lubricated point 
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contact can be written as 


where 
h 


o 

S(X,Y) 


h(X,Y) = h Q + S (X,Y) + w(X,Y) 


(3.13) 


constant 

the approximate geometric separation of the ellipsoidal 
solid and plane 
w(X,Y) elastic deformation 

The significance of the elastic deformation relative to the geometric 
separation of the ellipsoidal solid and plane can be expressed as 

R 2 = s (3.14) 


3.4 Conditions Investigated 

Figure 3.3 shows that we need to be concerned with the following 
questions : 

(1) How fine must the divisions of a and b be? In this 
chapter we assume that the number of divisions of a and b will 
be the same. Therefore, we can define the number of divisions as 


d = 


(3.15) 


2a 2b 

In this thesis we let d equ..l 3, 4, and 5. 

(2) How far from the semimajor and semiminor axes does 
(eq. (3.14)) become insignificant? In this study was evaluated 

at distances from the center of the contact of four times the semi- 
major and semiminor axes. 

To check the accuracy of the elastic deformation results for 
d of 3, 4, and 5, the numher of equal divisions along the semimajor 
and semiminor axes was Increased by three times (d of 9, 12, and 
15), and then corresponding points were compared. The following 
equation describes the percentage accuracy of the results compared 
with the finest -mesh-size predictions: 
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R 


3 



x 100 


(3.16) 


The limiting conditions that were evaluated on a computer are 
shown in table 3.1. It was speculated that conclusions that could 
be made for these limiting conditions could also be made for any in- 
termediate conditions. The four limiting conditions shown in ta- 
ble 3.1 are two extremes of applied normal load; a light load of 
8.964 newtons (2 lbf ) , and a heavy load of 896.4 newtons (200 lbf). 
The two extremes of curvature of the solids shown in table 3.1 are 
equal spheres in contact and a ball and outer race of a ball bearing. 
The elliptical eccentricity parameter (k = a/b) is 1 for the equal 
spheres in contact and 5 for the ball and outer race. 

The equations thus far developed were programmed on the Leeds 
University International Computers Limited (ICL) model 1906A digital 
computer. 

3 . 5 Discussion of Results 

Tables 3.2 to 3.13 give the characteristics of the deformed 
shape of the contacting solids along the semimajor and semiminor 
axes when the axes are divided into three, four, and five equal di- 
visions and the conditions of table 3.1 prevail. Some observations 
can be made about these tables: 


(1) Because of the coarse grid and the elliptical pressure pro- 
file, there is not much decrease in pressure in going from the in- 
nermost to the outermost point within the contact area, 

(2) The agreement of S with S is seen to be good and is 
borne out by the ratio of the two expressed in terras of , jhe 
biggest disagreement is in table 3.7 where R 1 » 0.9870, which means 
that S is in agreement with S within 1.3 percent. Because of 
this good agreement, S will he used to define the geometric sepa- 
ration of the ellipsoidal solid near a plane. 
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(3) The ratio of the elastic deformation to the geometric 

separation of the contacting solids is seen to decrease substantially 
with increasing distance from the center of the contact zone. Fur- 
thermore, the predictions of the distance at wh J ch the elastic de- 
formation becomes insignificant compared with the geometric separa- 
tion of the solids do not change whether we have two equal spheres 

or a sphere and an outer race in contact. 

(4) The separation due to the geometry of the contacting solids 
plus the elastic deformation (S + w) is almost constant in the con- 
tact zone. The value of (S + w) at the farthest point from the 
center of the contact zone and yet still within it differs the most 
from the other values of (S 4- w) in the contact zone. 

(5) The percentage difference in elastic deformation calcula- 

tions for two mesh sizes differing by a factor of 3 was shown to be 
small. For the worst case in table 3.11, was equal to 7.494 per- 

cent. That is, the elastic deformation for d = 3 differs from the 
elastic deformation at corresponding points when u = 9 by 7.5 per- 
cent, which is extremely good, 

(6) Comparing table 3.2 with 3.5, 3.3 with 3.6, and so forth, 
which amounts to changing the normal applied load from 8.964 newtons 
(2 lbf) to 896.4 newtons (200 lbf ) , leads to the following conclu- 
sions : 

(a) R 7 does not change in the corresponding tables. That is, 
regardless of the normal applied load, the ratio of the elastic de- 
formation to the geometric separation of the solids is unchanged. 

(h) R.j does not change in the corresponding tables. This con- 
dition is undoubtedly because of the condition mentioned in (a). 

To better illustrate the results shown in the tables, figures 
3.4 to 3.9 are presented. In figures 3.4 to 3.6 t he solid curves 
represent the case of equal spheres in contact, which is represented 


2 ? 


by 0.5558 cm (0.2188 in.); and the dashed curves represent 

the ball and outer race in contact, which is represented by R = 

x 

1.286 cm (0.5055 in.), R^ = 15. Ou cm (5.406 in.). Also as n result 
of me observation made in discussing the tables that R 2 and R^ 
are not functions of the normal applied load, the results shown in 
figures 3.4 to 3.6 apply for any normal applied load. 

Figures 3.4(a) and (b) show the effects of the location along 
the semimajor and semiminor axes, respectively, on the percentage 
difference in elastic deformation for d of 3 and 9. Here an "edge 
effect" can be seen, which is a rapid rise in percentage difference 
in the film thickness when d = 3 and for corresponding points when 
^ This rapid rise is due to the pressure being either zero if 

the center of the rectangular area shown in figure 3.3 is outside the 
contact zone or of the order of 10 5 if the center of the rectangular 
area is within the contact zone. However, it is speculated that in 
lubricated contacts, where the pressure gradients are, ir general, 
more gradual than those encountered near the edge of a dry Hertzian 
contact, this edge effect is likely to he less significant. Also 
note that outside the contact zone the value of R ^ decreases. 

Figures 1.5(a) and (b) show the ettect of the location along 
the semimajor and semiminor axes, respectively, on the percentage 
difference in elastic deformation for d of 3, 4, and 5 and the more 
exact elastic deformation for d of 4, 12, :u | 15. These figures 
show a large drop In R, } from d - 4 to d = 5, which also brings 
down the edge effect considerably. There is, therefore, good case 
for letting d B 5 in any further computer evaluations. 

Figures 3.6(a) and (b) show the effect of the location along the 
semimajor and semiminor axes, respectively, on the ratio „t the elas- 
tic deformation to the geometric separation of the ellipsoidal solid 
near a plane. These figures show the distance t r<»m the semi major and 


23 


semiminor axes at which the elastic deformation becomes insignifi- 
cant, To be more specific, from the curves we see that for equal 
spheres in contact (represented by solid lines in the figures) 

< 0.05 corresponds to x > 2.6 b and y > 2.6 a. Thus, the 
elastic deformation is less than 5 percent of the geometric separa- 
tion at a distance from the center of the contact zone that is no 
less than 2.6 times the semimajor or semiminor axis. For the ball 
and outer race in contact, R£ < 0.05 corresponds to y > 1.9 a 
and x > 4.0 b. In other words, the elastic deformation is less 
than 5 percent of the geometric separation at a distance of only 1.9 
times the semimajor axis and 4.0 times the semiminor axis from the 
center of the contact zone. 

Figures 3.7(a) and (b) show the effect of the location along 
the semimajor and semiminor axes, respectively, on the geometric 
separation of the ellipsoidal solid near a plane plus the elastic 
deformation when the load is 8.964 and 896.4 newtons (2 and 200 lbf) 
and the ellipticity parameter (k) is 1 and 5. The conditions men- 
tioned in this figure correspond to those of table 3.1. These fig- 
ures show the sum of the film thickness components to be essentially 
constant within the contact zone. 

Figure 3.8 shows the effect of the number of divisions across 
the ellipse axes on the computer time for running all four conditions 
shown in table 3.1. Here we see that the computer run time quickly 
becomes exorbitant as d is increased. This is why only selected 
data for R^ were obtained. 

3.6 Concluding Remarks 

A numerical analysis of the elastic deformation of a contacting 
ellipsoidal solid and plane has been performed. The analysis as- 
sumed that the pressure in the contact zone was Hertzian. It also 
assumed that the contact zone could be divided into rectangular areas 
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with uniform pressure within each rectangular area. The resulting 
equations were programmed on a digital computer. Four limiting con- 
ditions were evaluated on the computer. They consist of two ex- 
tremes of applied normal load: a light load of 8.964 newtons 

(2 lbf), and a heavy load of 896.4 newtons (200 lbf). The two other 
extremes are of the curvature of the contacting solids: two equal 

spheres in contact, and a ball and outer race of a ball bearing. It 
was speculated that conclusions that could be made for the limiting 
conditions could also be made for any intermediate condition. 

The results indicate that division of the semimajor and semi- 
minor axes into five equal subdivisions is adequate to obtain ac- 
curate elastic deformation results. Also the elastic deformation 
becomes insignificant compared with the geometric separation of a 
sphere near a plane at a distance from the center of 2.6 times the 
semimajor axis. For a ball and outer race in contact, a similar 
observation applied at a distance from the center of 1.9 times the 
semimajor axis and 4.0 times the semiminor axis. Finally, the geo- 
metric separation plus the elastic deformation (S + w) was almost 
constant in the contact region. However, numerical values of 
(S + w) at points near the edge of the Hertzian contact show that a 
slight edge effect or error may be encountered in such regions. 

In lubricated contacts, where the pressure gradients are, in gen- 
eral, more gradual than those encountered near the edge of a dry 
Hertzian contact, this effect is likely to be less significant. 
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CHAPTER 4 

THEORETICAL FORMULATION OF THE ELASTOHYDRODYNAMIC 
LUBRICATION PROBLEM 

The elastic deformation model was developed, along with an 
appropriate nodal structure, in chapter 3. The theoretical for- 
mulation of the elastohydrodynamic lubrication problem is given in 
this chapter. The procedure here is to give the Navier-Stokes equa- 
tions of motion as well as the continuity equation. Then with the 
assumptions that are imposed on the problem, these equations are re- 
duced to the Reynolds equation. The pressure-viscosity formula of 
Roelands (1966) is used. In the numerical analysis of the Reynolds 
equation, a <j> analysis (where <f> is equal to the pressure times 
the film thickness to the 3/2 power) is used to help the relaxation 
process. A standard finite difference representation is applied to 
the various terms in the Reynolds equation. By applying a Gauss- 
Seidel relaxation method to the finite difference form of the Rey- 
nolds equation, a converged solution of <f> for the complete nodal 
structure is obtained. When is known, the pressure at the 
various nodes can be obtained. When the pressure is known, the in- 
tegrated load is calculated, and adjustments are made for the 
initially guessed film constant (H^) until the integrated load and 
the input load are in agreement. A flow chart of the computer pro- 
gram is given, as well as the complete computer FORTRAN listing. 

4. 1 Reynolds Equation 

The derivation of the Reynolds equation governing the pressure 
distribution in an elastohydrodynamic lubricated (EHL) conjunction 
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is based on the Navier-Stokes equations of motion and the continu- 
ity equation. The most general form of the Navier-Stokes equations 
of motion for a Newtonian fluid in Cartesian coordinates as obtained 
from Pai (1956) is 
Du r, Bp 


Dt 


Dv 

Dt 


= Y - + 


3Y 


Dt 
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(4.1) 


(4.2) 


(4.3) 


where 


d _ a , a , a , - a 

Dt at + u ax + v :y + w az 


(4.4) 


In equations (4.1), (4.2), and (4.3) the left side corresponds to 
the inertia terms and the right side contafnc the body force, pres- 
sure, and viscous terms, in that order. The following assumptions 
are made in solving the EHL point-contact problem: 

(1) No external forces act on the film. Thus, 

X ~ Y = Z = 0 

(2) Fluid inertia is small when compared with viscous shear. 
These inertia forces are associated with acceleration of the fluid. 
Thus, 

Du _ Dv _ Dw 
Dt “ Dt Dt 

(3) There is no variation of pressure across the fluid film. 


Thus , 


BZ 


0 


27 


(4) The viscosity (n) and density (p) are constant in the Z 
direction . 

(5) The radius of curvature of the solids bounding the oil film 
is large when compared with the thickness of the lubricant films. 
This assumption allows any effects due to curvature of the oil film 
to be neglected. 

(6) There is no slip between the fluid and bounding solids at 
common boundaries. Thus, 


Z = 0 , u=u^, v=0 


7 -h, u=u B , v = 0 

(7) Because of the geometry of the fluid film, the derivatives 
of u and v with respect to Z are large when compared with all 
other velocity gradients. 

(8) Steady-state conditions are considered. 

With these assumptions, the Navier-Stokes equations of motion 
(eqs. (4.1), (4.2), and (4.3)) reduce to the following: 


_|£ d U 


_ A 

3Y " n a >: 2 


i£ - 


c)p 

— = Q 

3Z 


(4.5) 


(4.6) 


(4.7) 


The velocity distributions (u) and (v) can be found by integrating 
equations (4.5) and (4.6) twice while using the boundary conditions 
according to assumption 6. 

The equation of continuity, representing mass conservation, is 

2ii?- u J + l(£ v) 


:>x 


3Y' 


0 


(4.8) 


Coupling the velocity distributions with the equation of continuity 
leads to the following equation, which was developed by Reynolds 
(1886) and carries his name: 
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&) + - 12u p§ ( Ph) 


3X\ n SX) + 3Y[ n 3Y 


(4.9) 


where 




Letting 


X* = £, Y* = - p = _£L 

b a M p 


n " h , p = 4 


n = — , H = 

"O R x 


(4.10) 


(4.11) 


equation (4.9) can be rewritten as 

dx*fr ax*) + ^ m) = 12u (r^^ 

where 

k ellipticity parameter, a/b 

0 dimensionless speed parameter, n u/E'R 

0 J x 

Equation (4.11) is the Reynolds equation in dimensionless form, 
nor al requirement is for dimensionless parameter (P) to be de- 
termined. Before proceeding, however, the dimensionless density (p ) , 
the dimensionless viscosity (q) , and the dimensionless film thick- 
ness (H) need to be formulated. 

4 • 2 De nsity 

At the high pressures that exist in the elastohydrodynamic 
lubricating film, the liquid can no longer be considered as an in- 
compressible medium, and the dependence of the density on the pres- 
sure must be considered. From Dowson and Higginson (1966) the dimen- 
sionless density for mineral oil can be written as 


0 = 1 + _°i009JL__ 

1 + 0.026 p 


(4.12) 


where 


P gage pressure, ton/ in. 
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Therefore, the general expression for the dimensionless density can 
be written as 


P 


1 + 


ype' 

1 + $PE' 


(4.13) 


where B and y are constants dependent on the fluid. 

4. 3 Viscosity 

As long ago as 1893, Barus proposed the following formula for 
the isothermal viscosity-pressure dependence of liquids: 

= «P (4.14) 



where 


n dynamic viscosity at gage pressure (p) 

n Q dynamic viscosity at atmospheric pressure 

a pressure-viscosity coefficient of lubricant 

The pressure-viscosity coefficient (a) in equation (4.14) character- 
izes the liquid considered and depends only on temperature, not on 
pressure. Although equation (4.14) is extensively used, it is not 
generally applicable and is valid as a reasonable approximation only 
in a moderate pressure range. 

Because of the shortcomings of equation (4.14), several iso- 
thermal viscosity-pressure formulae have been proposed that usually 
contain two or more parameters instead of Barus 's (1893) single pa- 
rameter. One of these approaches, which is used in this thesis, is 
that of Roelands (1966), who undertook a wide-ranging study of the 
effect of pressure upon the viscosity of the lubricant. For iso- 
thermal conditions, the Roelands (1966) formula (p. 95) can be 
written as 


log n + 1 200 


(log n 0 + 1.200) 


1 + -B- V 

2000 / 


(4.15) 


*log denotes the common or Briggsian logarithm, log ; 
In denotes the natural or Napierian logarithm, log*? 


30 


where 


p gage pressure, kgf/ctn 2 

Z viscosity-pressure index, a dimensionless constant 
Taking the antilog of both sides of equation (4.15) gives 

n . 10 (1 ° 8 V 1 - 2 , ( 1 + 2 mo ) Z - 1 - 2 

Rearranging this equation yields 



lor the dimensionless form of the viscosity given in equation 
(4.10), this equation becomes 



Rearranging terms in this equation gives 



where 

= 6.31xlO" 5 N s/m 2 (0 . 0631 c?) 

and 6 is a constant equal to 19 609 N/cm 2 (28 440 lbf/in. 2 ). 

In equations (4.13) and (4.16), care must be taken to ensure that 

the same dimensions are used in defining the constants. 

In the Roe lands (1966) formulation the lubricant is defined by 

the atmospheric viscosity (n ) , the viscosity-pressure index (Z) , 

and the asymptotic isoviscous pressure (p ). The equation de- 

l v , a s 

scribing the asymptotic isoviscous pressure can be written as 


P. = n n f & 

iv, as 0 y, n 

HI ok (1964) arrived at the very important 
results achieved hitherto tor an exponent i 


(4.17) 

conclusion that all KHh 
ai viscosity-pressure do- 
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pendence (eq. (4.14)) can, to a fair approximation, be generalized 

for any given nonexponential dependence simply by substituting the 

reciprocal of the asymptotic isoviscous pressure (1/p, ) for the 

iv,as 

viscosity-pressure coefficient (a) occurring in those results. This 
implies that 


„ a- 1 
a ^ 


p . 

iv, as 


(4.18) 


It might be pointed out how the values of the parameters , , 2 , 

P iv,as are obtained for a given lubricant by using Che Roelands 
(1966) formulation. From table II-l (p. 48 of Ro e lands (1966))> for 

a paraffinic mineral oil, coded 31-6, the atmospheric viscosity (y 
is found to be 41.1 centipoise. From table IV-2 (p. 106 of Roelatl 
(1966)), for the same lubricant the viscosity-pressure index (2) is 
found to be 0.67. When the atmospheric viscosity <„„) is known, the 
following expression can be evaluated: 

log(log n 0 +1.2) = 0.4493 

Making use of this equation in table XI 1-2 (p. 652 of Roelands (1966)) 
for Z = 0.67 leads to 

‘° e p lv,as * 2 ' 695 
" P iv,as ' 't® 59 N/ ™ 2 

This, then, establishes the procedure used to obtain the parameters 

V Z ’ and p iv,as for •' 8 Ivon lubricant. 

4-4 Film Thickne ss 

As was .stated in chapter 2, the separation due to the geometry 
of the two ellipsoidal solids shown in figure 2.1 can be adequately 
described by an equivalent ellipsoidal solid near a plane. The ge- 
requiremcnt is that the separation of the ellipsoidal solids 
in the initial and equivalent situations should be the sa • at equal 
values of X and Y. Therefore, from figure 4. 1 , the separation due 
to the geometry of the two ellipsoids shown in figure 2.1 can be 
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written for an ellipsoidal solid near a plane as 


S(X,Y) 


= (X - mb) 2 (Y - la ) 1 

2R 2R 

x y 


where 

m constant used to specify length of inlet region 

l constant used to specify length of side-leakage region 

For purposes of illustration the mesh described in figure 4.1 will 
still be used, where m = 4 and £ = 2. However, the equations 
developed will be written in general terms. 

Figure 4.2 gives a physical description of the film thickness 
and its components for an ellipsoidal solid near a plane. Equation 
(3.13) describes how the components are related in describing the 
film thickness. In figure 4.2 it is assumed that Y is held con- 
stant near the midplane o r the contact. 

Substituting e .on (4.19) into equation (3.13) while at the 
saint: time making this equation dimensionless as in equation (4.10) 
gives 

H-H ♦ & 2 -Q g - -;■■«■>.* + + wSSJTil „ m 

O _2 2R R R 

2R x y x 

x J 

where is a constant that is initially estimated. 

In chapter 3 the elastic deformation (w in eq. (4.20)) of an 
equivalent ellipsoidal solid near a plane in contact and subjected to 
a Hertzian stress distribution has been evaluate! numerically. 
Therefore, by using figure 4.1 and the results o r chapter 3, the 
elastic deformation can he written as 


W 'n ( X * , Y * ) 

k, 


2?,xc (nrfn)xp, 

y y i>. .d 

/ v / ; i , 1 m , n 

J 5= 1 '> ! = 1 > 

J 1 * * * * - * • • « 


where 


n constant used to determine length of outlet region 
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and 


number of equal divisions in semimajor axis (in fig. 4.1, c " 5) 
number of equal divisions in semiminor axis (in fig. 4.1, d = 5) 
m = | k - i | +1 

n = |T - j| + 1 
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(4.22) 


Equation (4.23) points out more explicitly the meaning of equation 
(4.21) while making use of figure 4.1. The elastic deformation at 
the center of the rectangular area w (fig. 4.1) caused by the 

y > J 

pressure on the various rectangular areas in in around the contact 
ellipse can be written as 


w 


9,5 tt 


P l.l °9,5 + P 2,l D 8,5 + • 


+ P 35,l °27,5 


+ P, y D q . + P ? - D,. . + . . . + P , D-_ . 
1*^ 9>4 2,2 w , 4 35,.. 27,4 


+ P 1,20 D 9,16 + P 2,20 D 8,L6 + 1 ’ ' 


P I) 

35,20 27, 16j 


(4.23) 


4 • 3 Phi (j>) Solutio n 

Having defined the density, viscosity, and film thickness, we 


can return to the problem of solving the Reynolds equation. It is 
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well known (e.8., klhomes (19,6), that the dimensionless pressure (p) 
of the Reynolds equation p lo t te( i as » function of ej , MbUs a 

uliced pressure Held, giving high values of DP/9X* and 

' P/,X * 2 ' Such a con< *f c f°n with high gradients is not welcomed when 
performing a numerical analysis by means of relaxation methods. In 

order to produce a much Rentier curve, a parameter <*) is introduced, 
where 


hH 


(4.24) 


The pressure (P) is small at large values of film thickness (H) and 

conversely. This substitution also has the advantage of eliminating 

all terms containing derivative, of products of 11 and P or H 

and f. Therefore, by „ sl ng equation (4.24), equation (4.11) can be 
written as 

_ , _ 1 

1 


3 

;vx* 


W ' 2 & - 1 ♦»>« 4 




' /..3/2 _ 3 ,,,1/2 


.lY* 




0 H \ 


i inf ~ ViCiifi) 
Mr I ix* 


Expanding this equation yields 


U' " l/2 '"VI ■ 12l/- h Yi(-»H> 

k - ,v *vi : 'Y *JJ [rJ :<X^ 


(4.25) 


a • () U*lAt. c L JH t fj^eno e Koji r t-sjLnU ,i Lion 

The finite difference method will he used to develop the various 
terms in equation (4.25). Figure 4.3 shows the mesh to be used as 
related to the dimensionless coordinates X* and Y* . Kqu.it Ion 
(4.25) must he written for tin- point (i,j) In n Rurt . 4.3 by substitu- 
ting for the derivatives, expressions that involve values of 
f1, 1 ’ and n at the surrounding points. 

.. a function of 


At the three points X* , , X * .. and X* 


‘-‘.I’ "i.P i+1,. 

,n * u ‘ represented by a parabola, where 


X* such as 
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<|> = A(X*) 2 + BX* + C (A. 26) 

The parabola and corresponding points are shown in figure 4.4. The 
expressions for ^i,j* anc * ^i+l,j Can wr *‘" ten directly 

from equation (4.26) as 

^_i 




V i+l.J 

From figure 4.4 the following equations can be written: 

Y* = X* + — 

X i,j X i-l,j + d 


■ A < x !-i.j ,2 + BX i-i.j + c 

(4.27) 

- Atxjj) 2 + BX; J + C 

(4.28) 

- A < x l + i,j )Z + Bx m,j + c 

(4.29) 


X* . = X* 
1 + 1,1 


1-1.1 + d 


Substituting these expressions into equations (4.28) and (4.29) gives 

. 2 


1>. 


,j ■ A ( x * i - i . j + a ) + + 1 ) + c 

■i.j ■ + i ) 2 + 'K-i.i + d) + c 


(4.30) 


(4.31) 


Therefore, given equations (4.27), (4.30), and (4.31), it is possible 
to solve for A, B, and C to give 
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The following derivatives can he written by using equations (4.32), 
(4.33), and (4.34): 
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Having developed those basic quations with the dummy variable 
(i^), we can now proceed to develop the various terms in equation 
(4.25) by using the finite difference format developed in this sec- 
tion. The following equation is written for the point i,j: 
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From equations (4.38) and (4.35) the following equations can be 
wr i t ten : 
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1 / ; 1 /? .(1 \ ;> 
,X*(n H i X* ) - ,K2S '* 


i41J v»- s 


(4.41) 


(4.42) 


itl, 


1+1. i 1+1. i 


411. . + 11. , .) 

i . ) 1-1,1 


+ 1_1,j %/n. 

•i, , , V i-l,) 


i-l . 


Ill,,, . - 411, . + 311, . .) 
1 + 1,1 i . 1 i-l . 1 


( 4 . 4 3 ) 


Tims, tin* following equations nan in- ilirectlv written: 
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_1 Jg lt\ 

3X*\r] ax*/ 


= 0.25 d 


ffi+1 rl 

K+l.j 


(34*.,, . - 44. . + <J> . .) 

i+l»J i,J i-l,j 




(4.44) 


9 (pH) = 0.5 d( PjJ1 .H. 


ax * v... - ... u„, 1+1> .n i+1> . - P^.H.^.) (4.45) 

The derivatives with respect to Y* can be obtained directly by sub- 
stituting c for d, the subscript i,j~l and i-l,j, and the sub- 
script i, j+1 for i+1 , j . 


J+l (3H i,j + l- 4H i,j +H i,j-l) 


-J-fe h 1/2 JH\ _ 2f p i,j+l 

3Y*(n H 9Y*j 0,25 c l/f. , +l Vi^ 

+ 77 1,J 1 V** 7 TT (H -i i+i - 4H - • + 3 H . . ,)1 

n i,j-l 1,J 1 1,J+1 *.J i.J-1 

- : l/£ ja.\ = o 25 c 2piu±i fu 

8Y*^n dY*J 5 c n (3 ^i ■ 

' L i.J+1 ’* 


(4.46) 


. , “ 4(j) * . 4* (f> . . ) 

.J+i i.j i,j-r 


+ ^ ( *i^ + i" 4 *i.J + 3 *i.J- 1 ) 


(4.47) 


Substituting equations (4.43) to (4.47) into equation (4.25) while 
collecting terms gives 

h.jh+i.j + Vjh.j-i + c i,jh-i,j + 5 i,j*i,j+i -h.jh.j -"i.j ■ 0 

(4.48) 


where 


„ • £ 


A i.j “ ^i+i.j + Vi.j 


“i.j ■ (w ( "i,jti + 3 “i,j-i> 


C i,j “i+l,J + -‘"i-l.j 


(4.49) 

(4.50) 

(4.51) 

(4.52) 
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D i,j (dk) (3ll i,jxl + Vj-l } 


JtI i> j-1 


4(l 'i+l,j + Vl.j) + 4 (dk) (U i,j+i + Vj-i) 


Ml \ M 1+1 i ^ (3H . - 4H. . + H 

2H J 1 l 1+1 ’ 3 v 1+ 1>J i+l, J i,j i-1 , j ' 

i, j 


4H . . + H. , .) 

-t i i . 1 4 7 


■i-i.j VWo - "tj + + (£) 




- 4H. . + H. . 

i -j A 4 i * 


j +i i,j i,j-i 


J i,j-1 V H i,j-l (n i,j+l " 4H i,j + 3H i,j-l>]] (4 - 54 > 


24Ub 

dR H 3/2 P i+l.j H i+l,j " P i~l,j H i-l,j ) 
x i,i 


(4.55) 


4.7 Nodal Structure 


The nodal structure used in obtaining most of the results is 


shown in figure 4.5. This nodal structure was arrived at after much 


exploration in which the number of nodes in the semimajor and semi- 


minor axes, as well as the distance from the center of the contact 


to the edges of the computing zone, was varied. The . odal structure 


shown in figure 4.5 was considered to be suitable when the minimum 


film thickness did not change when either additional nodes were placed 


in the semimajor and semiminor axes or the distance from the center 


of the contact to the edges of the comp ting zone was extended. 


From figure 4.5 the following can be written: 


m = 4 


n = 1.15 


l = 1.6 


(4.56) 


c = 5 


d = 13 



/ 




W-^7 
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These values used to define the nodal structure were used for most of 
the evaluations. One exception was for high-speed cases when the 
constant used to determine inlet length (m) had to be extended from 
4 to 6, with the other values being held constant. 

4 . 8 Bo undary Conditions 

The boundary conditions are the following: 

(1) At the edges of the rectangular computation zone (fig. 4.5) 
the pressure is zero, which therefore implies that <f> is also zero. 
Specifically, this means that along the bottom of figure 4.5, 

4*1,1 = along the left side, <j>^ ^ = 0; along the top, ^ ^ = 0; 

and along the right side, <)>,_ . = 0. 

6/ > J 

(2) At the cavitation boundary. 


P 


3P = JfP_ _ 
3X* 3Y* 


This condition is commonly known as the Reynolds condition and will 
be satisfied by simply resetting <j>. . to zero whenever it occurs 

1 9 J 

as a negative value. 

4 . 9 Initial Conditions 


Outside the Hertzian contact ellipse the pressure is initially 
assumed to be zero, and therefore 4> is also equal to zero. That 
is , ip = 0 when 

(X* - m) 2 + (Y* - 5.) 2 > 1 

Inside the Hertzian contact ellipse the pressure is initially as- 
sumed to be Hertzian; that is, 

P = 2^bP V* “ (X * " m)2 " (Y * " 1)2 (4 . 57) 

or 

3FH 3 ^ 2 f ~ 2 ' ~~ 2 

1 = 2 '^bF~‘ V 1 “ (X ~ m) " (Y * - (4.58) 


(X* - m) 2 + (Y* - J?.) 2 < 1 


when 
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4.10 Relaxation Method 

If the subscript n is the iterant and <J>. . is the particular 

1 > J 

solution to be found, the relaxation method known as the Gauss- 
Seidel method can be expressed as 


<|> . . . - = 4 > . . 

i,j,n+l i,j,n 


- x(- ^)A +1 i l li n ~ C i, j^j-1, j ,n+l D i..i^i.i+l,n 


- C. 


L i,i 


where X is an overrelaxation factor that is initially set to 1.6. 
Therefore, equation (4.59) is used, starting with node (2,2) and then 
continuing with (3.2), . . ., until (M,2) followed by 
(2,3), (3,3), . . ., (M, 3) , and ending with (2,N), (3,N) . . ., (M,N), 


where 


M = (n + m)d - 1 
N = 2£c - 1 


(4.60) 


The relaxation procedure described by equation (4.59) is con* 
tinued until 


N 

z 

= 9 9 


i 


^■J.j >n+l ^i,,1 ,n j 
‘h, j ,n+l 


< 0.1 


4.11 P ressure Loo p 

The relaxation method provides values of A . for every point 

i » J 

within the mesh. Having determined <J> , , we can write the dimen- 

i > J 

sionless pressure as 


P. . 





-3/2 


(4.61) 


With these new values of the dimensionless pressure, new values of 
the dimensionless viscosity, density, and film thickness can be 
evaluated. Thus, the coefficients of equation (4.59) 
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w 


(A. . , B. . , . . L. .) should also change. Accordingly, it is 

i * J i> J i * J 

necessary to reenter the relaxation loop. This process is con- 
tinued until the following inequality is satisfied: 

N M 


z z 


j=2 , 3, . . . i=2 , 3 , . . . 
4.12 Normal Applied Load 


p - P 

K Qml 

P i, j ,n+l 


The constant in equation (4.20) was initially estimated. 


The next task is then to find the correct value for H . In order to 

o 

do this, the integrated normal applied load must be evaluated, where 


/ nrl-n r 21 

J PdY dX 


(4.62) 


'0 Sq 

By applying Simpson’s rule, this double integral can be rewritten as 


4 K a b 
Vrd 


J(',T 
► CJP 


j 




3, i 


+ :n\. 
f zr 


+ 2P , 




3,4 


i, 3 


3,z--z 


+ 


J ( Z P , , + V , + ZP . 

m+n -J,. nr+n- . . i nrhi- 


n+u- J , r 


+ •’ 1’ , + i’ , , ) 

nrhi-Z,. n+n- . ! 


f (ZP 


f !> 


■M- 1 , Z m+n- 1 , J ‘ nrhi-1,4 ‘ m+n-1,3 


,’P + i> 

irr+n- 1 . • - J m+n.- 1 , Z - - 1 


This equation can be written simply as 

m+n-1 

2 

i=2","3 , . . 

where 


F = 4El§b 
9cd 


m+n-1 

z 


a / 2P-1 

t -» / \ — v q 0 

2 "P . 


‘ E 


(4.63) 


q l = 0 

if 

i 

is 

odd 

q l = 1 

if 

i 

is 

even 

q 2 = 0 

if 

j 

is 

odd 

% = 1 

if 

j 

is 

oven 


By using the trapezoidal rule, the normal applied load ran also he 
expressed as 


iw _ 

wtmamm 


UKPROnUCIBIU TY OP Tlltf 
‘ A'. r.\t;r; .y poor 
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F 


itH-n-1 



1=2,3,... 


21-1 

S 


i = ? 


,3,.. 


P. . 
ij 


(4.64) 


Both equations (4.63) and (4.64) were programmed to serve as a check 
on one another. 

When the current values of H q and F (eq. (4.63)), as well as 
the initial normal applied load (F) are known, a new estimate for 
(H^) can be expressed as 


i = H ($) 
o o\F/ 


0.13 


(4.65) 


With a new value for H q the film thickness (eq. (4.20)) is recal- 
culated, and reentry into the relaxation process is required. This 
process is continued until the following convergence is satisfied: 


^< 0 . 


0005 


Once this convergence criterion has been satisfied, the pressure and 
film thickness in and around a point contact are established. 

4.13 Flow Chart s 

Figures 4.6 and 4.7 are flow charts for the numerical solution 

on the digital computer of the equations developed in the analysis. 

Figure 4.6 is the flow chart of the main program. There are essen- 

tially three loops within the main program: In the relaxation loop, 

^i j n+3 * S 8 enera ted. In the pressure loop, the new values of 

4>i j n+ j of the relaxation loop result in new values of pressure 

P . * * which in turn result in new values of film thickness H 

bj,n+l i,j* 

new values of viscosity n . . , and new values of density p. . The 
final loop is the normal-load loop, which ensures that the integrated 
normal applied load agrees with the initially specified value. 

Figure 4,7 is the flow chart of the subroutine SUB6. Here it 
can be seen that a number of calculations occur only once and need 
not be repeated on reentering this subroutine. With a new pressure 


43 


distribution the elastic deformation is recalculated and this sub- 
routine is left with a new film thickness and, therefore, a new 

*i,j,n+i' The F0RTRAN listings of the main and SUB6 computer pro- 
grams are given in the appendix. The FORTRAN listings of the EHL 
point-contact analysis were programmed on the NASA Lewis Research 
Center's Sperry Univac 1100/42 computer, which has 131 K primary 
memory and 262 K extended memory. 

4 • 14 Concluding Remarks 

In this chapter a procedure for the numerical solution of the 

complete, isothermal, elastohydrodynamic lubrication problem for 

point contacts is given. This procedure calls for the simultaneous 

solution of the elasticity and Reynolds equations. In the elasticity 

analysis the contact zone was divided into equal rectangular areas. 

It was assumed that a uniform pressure was applied over each area. 

In the numerical analysis of the Reynolds equation the parameter 
3/2 

(() - PH , where P is the dimensionless pressure and H the di- 
mensionless film thickness, was introduced in order to help the re- 
laxation process. The nodal structure, boundary conditions, and 
initial conditions were given and a Gauss-Seidel relaxation method 
was used. The computer program has three major loops: the -elaxation 

loop, the pressure loop, and the normal-load loop. The last loop re- 
quires the integrated load to be in agreement with the input load 
within some tolerance. Upon the convergence of all three loops the 
pressure and film thickness in and around a point contact are estab- 


lished. 






CHAPTER 5 

FULLY FLOODED RESULTS 

The most important practical aspect of the elas tohydrodynamic 
lubrication (EHL) point-contact theory developed in chapter 4 is 
the determination of the minimum film thickness within the contact. 
That is, maintaining a fluid film of adequate magnitude is extremely 
important to the operation of some machine elements. In the present 
chapter the influence of contact geometry as expressed in the el- 
lipticity parameter and the dimensionless speed, load, and material 
parameters on minimum film thickness is investigated for a conjunc- 
tion fully immersed in lubricant (i.e., fully flooded). In the nu- 
merical work the ellipticity parameter is varied from 1 (a ball-on- 
plate configuration) to 8 (a configuration approaching a line 
contact). The dimensionless speed and load parameters are varied 
over ranges of about two and one orders of magnitude, respectively. 
Conditions equivalent to using solid materials of bronze, steel, and 
silicon nitride and lubricants of paraffinic and naphthenic mineral 
oils are considered in obtaining the exponent on the dimensionless 
material parameter. Thirty-four different cases are used in obtain- 
ing the fully flooded minimum-f ilra- thickness formula. A central- 
film- thickness formula is also developed. 

In this chapter, contour plots are shown that indicate in detail 
the pressure spike and the two side lobes in which the minimum film 
thickness occurs. These theoretical solutions for film thickness 
have all the essential features of previously reported experimental 
observations based on optical interferometry. 


I: 
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5« 1 Dimensionless Groupin g 

The variables resulting from the isothermal EHL point-contact 
theory developed in chapter 4 are 


X 

R 

y 

h 

E r 


effective radius in x-direction, mm 
effective radius in y-direction, mm 
film thickness, mm 
effective elastic modulus, N/mm 2 


P iv,as «y»Pto t ic isoviscous pressure obtained from Roelanda (1966) 
N/mm 2 

U surface velocity in x-direction, mm/sec 

°0 atmospheric viscosity, N • sec/mm 2 

viscosity pressure index, a dimensionless constant 
F normal applied force, N 

°’ B C0 " Sta " tS used “ density of flufd, ™2 /N 

It has been found by Dowson and Higginson (1966) that density has 
little effect on minimum film thickness for line-contact situations; 
therefore, one may assume the same is true for point-contact situa- 
tions. Even though the compressibility effect is still considered 
in the EHL theory developed in chapter 4, the constants used to de- 
fine the fluid in the density equation wiJJ. not be used in the 
minimum-film-thickness formulation. Therefore, the 11 variables men- 
tioned above were reduced to nine, a and 6 being eliminated. From 

the nine variables the following five dimensionless groupings can be 
written: 


(1) Dimensionless film thickness 

h 


H = 


R 


where 


-i- + ~L_ 

r Ax r Bx 


(5.1) 


(5.2) 


The radii of curvature in equation (5.2) 


are shown in figure 2.1. 
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(2) Dimensionless load parameter 

F 


W = 


E'R 2 

X 


where 


E' = 


- v 


A B 

(3) Dimensionless speed parameter 


U = 


v 

E'R 


where 


(U A + U B> 


(4) Dimensionless material parameter 

G = E ’ 


(5.3) 


(5.4) 


(5.5) 


(5.6) 


where p . 

iv, as 


P iv, as 

is the asymptotic isoviscous pressure obtained from 
Roelands (1966) . The asymptotic isoviscous pressure can be approx- 
imated by the inverse of the pressure-viscosity coefficient 

(p . % 1/a) . 

r iv,as 

(5) Ellipticity parameter 


(5.7) 


k = ~ (5.8) 

where 

a semimajor axis of contact ellipse 

b semiminor axis of contact ellipse 

The ellipticity parameter is determined entirely from the definition 
of the radii of curvature (r^, r fix> r^, and r fiy ) , and the deriva- 
tion can be found in chapter 2. 

The dimensionless film thickness can be written as 


H = f (k , U » W, G) 


(5.9) 
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The most important practical aspect of the EHL point-contact theory 
developed in chapter 4 is the determination of the minimum film 
thickness within the conjunction. Therefore, in the fully flooded 
results to be presented in this chapter the dimensionless parameters 

(k, U, W, and G) will be varied and the effect upon minimum film 
thickness will be studied. 

^ Effect of Ellipticity Parameter 

The ellipticity parameter (k) is a function of the radii of 

curvature of the solids only (r r r and r ’i tw j<4 

X Ax’ r Bx’ Ay* 3nd r By )> The radit 

of curvature in the x-direction for both solids A and B are used 
in defining the dimensionless speed and load parameters. Therefore, 
only the radius of curvature of solid B in the y-direction will 
be changed in varying the ellipticity parameter from 1 (a ball-on- 
plate configuration) to 8 (a configuration approaching a line con- 
tact). In doing this the dimensionless speed (U) , load (W) , and 
material (G) parameters were held constant at the following values: 


-11 


(5.10) 


U - 0.168jxl0 
w = 0.1106X10 -6 
0 = 4522 J 

Care was taken to ensure that the highest ellipticity parameter 
(k - 8) was in the elastic region. The approach used was to relate 
the opt.ating conditions for this near- 1 ine-contact situation to the 
clastic region defined in figure 5.i, which was obtained from Dowson 
and Whitaker (1965). From this figure the elastic, intermediate, 
and rigid regions are defined for G = 5000 (which is very close to 
the value of G in eq. (5.10)) and for various values of dimension- 
less speed (1!) and dimensionless load (U {) ) . The dimensionless speed 
(U) in Dowson and Whitaker (1965) is exactly that used in this thesis. 
The dimensionless load (W () ) used in figure 5.1 differs from that used 
in this thesis, and a tie between the two must be developed. The 
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dimensionless load parameter as defined by Dowson and Whitaker (1965) 
is 


W 


D 



(5.11) 


where 


p - Force 

Unit length 

Therefore, in relating this to an elliptical contact in which the el- 
lipticity parameter (k) is large compared with unity, a dimensionless 
load parameter can he written as 


W D 2a E' R 


( 5 . 12 ) 


X 

For this equation and k = 8 the location of U and Wp is shown 
by a circular symbol in figure 5-1. As can be seen, it lies within 
the elastic region. Thus, for the conditions given in equation 
(5.10) and for a ellioticity parameter (k) of 8, the conjunction is 
truly elastohydrodynamic. For k less than 8 the results move 
further into the elastic region. 

Table 5.1 gives 10 values of k and the corresponding minimum 
film thickness (H^ ) as obtained from the EHL point-contact theory 
developed in chapter 4. Having these 10 pairs of data, the object 
is to determine an equation that describes how the ellipticity pa- 
rameter affects the minimum film thickness. The general form of this 
equation can be written as 


H 


nun 


H 


Ae 


Bk 


(5.11) 


min, Ly 

A least-square exponential curve fit to the 10 pairs of data points 


was used 
Besides a 



in obtaining values for A and B 
least -square fit a coefficient of 


. 10 

1 n equat ion (5.1 1) 
> 

determination (r*“) 


is 
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obtained. The value of r 2 reflects the fit of M,e data to the re- 
sulting equation: ] being a perfect fit, and zero the worst pos- 
sible fit. The minimum film thickness for a line contact (H ) 

min t 1 / 

used in equation (5.13) was determined by finding the H that 

min , L 

gives a coefficient of determination closest to 1. The line-contact 
minimum film thickness was thus deduced from the present set of 
results for the limiting case in which the ellipticity parameter (k) 

approached infinity. This value of H turned out to be 

*“ 6 * J 

7.082x10 with a corresponding coefficient of determination of 

0.9990, which is an excellent fit. Furthermore, the values of A 
and B in equation (5.13) as obtained from the least-square fit are 

A ~ 0.9966 % 1.00 (5.14) 

B = -0.6752 % -0.68 (5.15) 

Fro„ Katies (5.11), ( 5 . 14 ), ..,„<) (5., 5 , tlw (ollowln,. proportion- 
ality ran bo written, which shows the effect of ellipticity parameter 
on minimum film thickness: 


ii . 
min 


/ 1 “ 0.68 k . 

U - L» ) 


(5.16) 


WiK>re H min = h min /R x’ predictt>l1 b Y the relationship that gives the 
best least-square fit to the numerical solutions. It is most sig- 
nificant that A turns out to be 0.9966, or approximately I .00, 
■shue as k • 0, • (). Therefore, even though tile smallest 

value of k used in obtaining equation (5.16) was unitv, it would 
seem that 'equation (5.16) could be applied to smaller values since 
in the limiting case (k * 0), equation (5.16) satisfies the physical 
situation. For the other extreme ot large k, a line-contact situa- 
tion is approached and the agreement with existing results is again 
good. From Dowson and Higginson (l'>66) the line-contact minimum 
Him thickness for the .1 imensiunless parameters given in equation 

r, *" )) ^ 7,7 -’ l)xln Compare this with 7.()H2xl(f b t rum the present 
results. I he dltlerence ot '» percent could Well he the result of 
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Dows on and Higginson (1966) using an exponential pressure-viscosity 

relationship instead of the Koelands (1966) formulation used in the 
present work. 

Substituting equations (5.14) ami (5.15) into equation (5.15) 
gives H mln , the dimensionless minimum film thickness obtained from 
the least-square formulation. The for the 10 elliptic, 'ty 

parameters are given in table 5.1. The percentage difference be- 
tween the minimum film thickness obtained from F.HI. point-contact 
theory (H^) and the minimum film thickness from the lea.sr -square- 

fit equation (H . ) is express, ’ as 
mm 1 


n _ / m i n min. . 

°i -| -‘h i 10, » 

min 


(5.17) 


bote that in table 5.1 the magnitude of b, never exceeds • 3 percent. 

figures 5.2 to 5.10 give contour plots of pressure and film thick- 
ness for ellipticity parameters of H, 6. 4, 5, 2.5, 2, 1.75, 1.5, and 
1.25, respectively. In these figures the other dimensionless param- 
eters (U, 0, W) were held fixed, as described in equation (5.10). 

1,1 partS (a) of LlK,Sl ' I * gures , rumour plots of dimensions pressure 
are given; in puts (b), eontour pints of d i mens , ou I ess film thi.k- 
tiess are given; and in parts (c), t firee-d i mens i ona I r. T r„ sen, at ions 
pressure are ,i..- . The 4 symbol in parts (a) and (b) indicates 
t,H ‘ ‘•“ nt4,r th< - "•rtxi.u -art. Note that, because of the di- 

mensionless represent.it i,>n ot the coordinates, the actual He tzfan 
nuitaet ellipse heroines a circle regardless ot the value of the el- 
lipticity parameter. The Her, .inn eon, act circle is shown in each n, 
parts (a) and (b) by asterisks. At the top ot ea.-h par, the con- 
tour label and the corresponding value are given. Ihe inlet region 
to tl 'e left and the exit region is to the right 


< 7 %, 
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For an ellipticity parameter (k) of 8, the maximum pressure is 
near the center of the contact (fig. 5.2(a)); and even though the 
conditions are in the elastic region, no pressure spike occurs. As 
the ellipticity parameter is decreased (going from fig. 5.2(a) to 
5.5(a)) the maximum pressure moves to the right (i.e., downstream 
of the center of the Hertzian conjunction) and increases in value. 

In decreasing the ellipticity parameter even more (going from fig. 

5.6(a) to 5.10(a)) the label of the pressure spike goes from A 
to C although the numerical values of the contours are different 
in each case. The pressure gradient is much larger at the exit end 
of the conjunction than at the inlet region. The pressure distribu- 
tion is seen even more dramatically in figures 5.2(c), 5.3(c), . . ., 

5.10(c), corresponding to ellipticity parameters of 8, 6, . . ., 1.25. 

The three-dimensional representation of pressure shows the details 
of the pressure spike more clearly. 

Figures 5.2(b), 5.3(b), . - ., 5.10(b) show contour plots of 
film thickness when the ellipticity parameter (k) is 8, 6, 4, 3, 2.5, 

2, 1.75, 1.5. and 1.25, respectively. For k = 8 the minimum film 
thickness is in a small area directly behind the axial center of the 
contact As k is decreased the minimum film area moves away from 
the axial center of the contact. For k = 1.25, two mini-mum-film- 
thickness areas occur at the sides nearer the edge of the Hertzian 
circle. These results, showing the "side lobes" in which minimum- 
film- thickness areas occur, produce all the essential features of 
previously reported experimental observations based upon optical in- 
terferometry. 

Figures 5.11(a) and (b) show the variation of pressure and film 
thickness, respectively, in the X-direction close to the midplane 


HKPK0D!K;;oi: 


of the contact for three values of the ellipticity parameter. As 
been ttue for all the ellipticity parameter results presented, 
the values of the dimensionless speed, load, and material parameters 
»ere held fixed as per equation (5.10). In f lgure 5 . ll(a) „ e flnd 

that for k - 6 no spike occurs, but this well-known feature of 
theoretical solutions to the elastohydrodynamic problem is evident 
for k = 2.5 and k = 1.25. 

In figure 5.11(b) for k - 1.25 the central region is not 
parallel with the X-axis. The reason is probably that compressi- 
bility effects are considered in the theory developed in chapter 4. 
That is, when compressibility is considered, the film thickness in 

the center is reduced by the amount that the fluid volume decreases 
at high pressure. 

5. 3 Influence of Speed 

By changing only the surface velocity in the x-direction (u) 
the dimensionless speed parameter (U) (eq. (5.5)) changes, but the 
other dimensionless parameters (k, W, and G) remain constant. The 

values at which these dimensionless parameters were held constant in 
the calculations are 


(5.1 


k = 6 

W = 0. 7371xl0~ 6 
C = 4522 J 

Table 5.2 gives the dimensionless speed parameter (U) and the 
corresponding minimum film thickness (H^) as obtained from the EH] 
point-contact theory developed in chapter 4. There are 15 different 
values of the dimensionless speed parameter covering nearly two 
orders of magnitude. Having these 15 pairs of data, the objective 
is to determine an equation that describes how the dimensionless 
speed affects the minimum film thickness. The general form of this 
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equation can be written as 

H min = IljK (5.19) 

By applying a least-square power fit to the 15 pairs of data 

^ U i’ H min,i^’ 1 = 1 » * • •» 15 >» the values of I and K were 
found to be 


I = 560.18 


(5.20) 


K 0.67542 0.68 ^ 2j^ 

The coefficient of determination (r 2 ) for these results was excellent 
at 0.9998. Substituting equations (5.20) and (5.21) into equation 
(5.19) gives the values of shown in table 5.2. The percentage 

difference (D^ between the minimum film thickness obtained from the 
EHL point-contact theory (H^) and the minimum film ob- 

tained from the least-squares fit (H^) is expressed in equation 

(5.17) and given in table 5.2. Note that the variation of D.. is 
less than ±2 percent. 

From equations (5.21) and (5.19) the effect of dimensionless 


speed on dimensionless minimum film thickness can be written as 

, 0.68 


h . * u 1 - 

rain 


(5.22) 

Figures 5.12 to 5.26 give contour plots of pressure and film 

thickness for dimensionless speed parameters ranging from 5.050xlO~ U 
-13 

to 8.416x10 . The other dimensionless parameters (k, W, and G) 

were held fixed as described in equation (5.18). As was true for 
the ellipticity results, contour plots of pressure and film thickness 
are designated a and b, respectively, and three-dimensional repre- 
sentations of pressure are designated c. 


In figure 5.12(a) the maximum-pressure area is near the center 
of the contact. As the dimensionless speed is decreased to and in- 
cluding U - 5.892x10 (fig. 5.21(a)), the maximum-pressure area 
moves downstream to the right. For dimensionless speeds greater than 
U = 5.892x10 1 the maximum-pressure area moves back tovard the cen- 


w 
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ter of the contact. Note that the "spike" pressure exceeds the 

Hertzian maximum pressure at U = 5.892xl0~ 12 , but at smaller speeds 

"".12 

(e.g., U = 4.208x10 ) the Hertzian pressure is dominant. Also 

note in these figures that the pressure in the inlet region is 
higher at high speeds than at low speeds. 

In figure 5.12(b), a high-speed situation, two minimum-film- 
thickness areas appear midway between the center of the contact and 
the Hertzian circle. As the speed decreases, the two miniraum-f ilm- 
thickness areas reduce to one, which is located at the axial center 
of the contact. With further reductions of speed the minimum-film- 
thickness area still remains in the axial center but moves closer 
to the Hertzian circle. 

In figure 5.12(c) a three-dimensional representation of the 
pressure for the highest speed case is shown. Here no pressure spike 
occurs, and the pressure rises to a peak that was shown in figure 
5.12(a) to be near the center of the contact. As the speed is re- 
duced the pressure spike emerges, and the top of the pressure pro- 
file becomes flatter. Once the pressure spike occurs, it moves 
toward the exit of the conjunction as the speed decreases. 

Figures 5.27(a) and (b) show the variation of pressure and film 
thickness, respectively, on the X-axis at the midplane of the con- 
junction for three values of dimensionless speed parameter. In fig- 
ure 5.27(a) the dashed line corresponds to the Hertzian pressure dis- 
tribution. Figure 5.27(a) shows that the pressure in the inlet region 
is higher for the high-speed (U = 5.050xl0 _U ) profile. For 
b - 0.8416x10 and U = 0.08416x10 the pressure spike orig- 
inates very near to the Hertzian pressure, and as the speed in- 
creases the pressure spike moves upstream. 

The typical elastohydrodynamic film shape with an essentially 
paralleJ section in the central region is shown in figure 5.27(b), 
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Also there is a considerable change in film thickness as the di- 


mensionless speed is changed, as indicated by equation (5.22). This 
illustrates most clearly the dominant effect of the dimensionless 
speed parameter (U) upon the minimum film thickness in elastohydro- 
dynamic contacts. 

5. 4 Influence of Load 

By changing only the normal applied load (F) in equation (5.3) 
the dimensionless load parameter (W) changes while the remaining di- 
mensionless parameters (k, U, and G) remain constant. The values 
at which these parameters were held constant are 

k = 6 "h 

U = 0. 1683xl0 -11 > (523 

C = 4522 J 


Table 5.3 gives the dimensionless load parameter (W) and the 
corresponding minimum film thickness (H^) as obtained from the EHL 
point-contact theory developed in chapter 4. There are eight differ- 
ent values of the dimensionless load parameter, covering over an 
order of magnitude. Having these eight pairs of data, the objective 
is to determine an equation that describes how the dimensionless 
load affects the minimum film thickness. The general form of this 
equation can be written as 

H min = QwL (5.24) 

By applying a least-square power fit to the eight pairs of data 

^i* ^min,i^* i ~ 1 * • • 8 }, the values of Q and 1 . were 

found to be 


Q = 2.1592x10"° ( 5>2 

L = -0.072924 % -0.073 ( 5>2 

The coefficient of determination (r 2 ) for these results was 0.9260. 
which was good, but was the lowest obtained in deriving the minimum- 
film- thickness equation ( 5 . 33 ). 


Substituting equations (5.25) and 
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(5.26) into equation (5.24) gives the values of H . shown in 

min 

table 5.3. The percentage difference (D^) between the minimum 

film thickness obtained from the EHL point-contact theory (H . ) 

min 

and the minimum film thickness from the least-square-fit equation 
^min^ ex P resse( i in equation (5.17) and given in table 5.3. In 
table 5.3 the variation of is within +3 percent at all times. 

From equations ( r .24) and (5.26) the effect of load on minimum 
film thickness can be written as 


H . * W 

mm 


-0.073 


(5.27) 

Figures 5.28 to 5.35 give contour plots of pressure and film 

thickness for dimensionless loads ranging from 0.1106xl0~ 6 to 

—6 

1.290x10 . The other dimensionless parameters (k, U, and G) were 

held fixed as described in the relationship (5.23). Once again, in 
parts (a) of these figures, contour plots of pressure are given; in 
parts (b), film thickness contours are presented; and in parts (c) , 
three-dimensional representations of pressure are given. 

In figure 5.28(a) the maximum pressure occurs directly behind 
the center of the contact, with no pressure spike occurring. As the 
load is increased the pressure spike emerges. With further increases 
of load the pressure spike moves toward the exit of the conjunction. 


These results can be more clearly seen in the three-dimensional 
representation given in parts (c) of figures 5.28 to 5.35. 

In figure 5.28(b), the low-load case, the minimum film thick- 
ness occurs directly behind the center of the contact. As the load 
is increased the minimum-film- thickness area still remains in the 
axial center of the contact but closer tc the Hertzian circle. With 
further increases of load, two minimum-film- thickness areas appear 
equidistance from the axial center and closer to the Hertzian circle. 
In figure 5.35(b), the highest load case considered, the minimum 
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film thickness is off to the sides in two areas close to the Hertzian 
circle. 

The variation of pressure and film thickness in the in- 
direction along a line close to the midplane of the conjunction is 
shown in figure 5.36 for three values of the dimensionless load pa- 
rameter. The values of the dimensionless speed , material, and el- 
lipticity parameters were held fixed as described by equation (5.23) 
for all computations at various loads. In figure 5.36(a) note that 
as the dimensionless load is increased the inlet pressure becomes 
smaller. For the highest load case shown in figure 5.36(b), film 
thickness rises between the central region and the outlet restriction 
in the same manner as seen in figure 5.11(b). Again this is attrib- 
uted to compressibility effects of the fluid. Also note that at a 
load of W = 0.5528x10 ^ the film thickness is slightly smaller than 
at a load of W = 1.106x10 The reason is that at the lower load 
the minimum film thickness is closer to the axial center of the con- 
tact than at the higher load. As was pointed out in discussing fig- 
ures 5.18(b) to 5.35(b), the location of the minimum-film-thickness 
region changes as the dimensionless load is changed. 

5. 5 Effect of Material Properties 

Contrary to what was found in the previous three sections, the 
effect of the dimensionless material parameter on minimum film thick- 
ness is not a simple matter. As can be seen from equations (5.3), 
(5.5), and (5.7), when either the material of the solids (as ex- 
pressed in E') or the lubricant (as expressed in n- and p ) 

0 r iv,as 

is varied, not only does the material parameter (G) change, but so 
do the dimensionless speed (U) and load (W) parameters. Only the 
ellipticity parameter can be held fixed. For all the results pre- 
sented in this section the ellipticity parameter is held fixed at a 


value of 6. 


Table 5.4 gives the four material- 


parameter results. The gen- 
eral form showing how the minimum film thickness is a function of the 
dimensionless material parameter is given as 


C = TG 


(5.28) 


where 


C = 


min 


d _ e “ 0 - 68 k )LJ 0.68 w -0.073 


(5.29) 


Note m equation (5.29) that the exponents are rounded off to two 
significant figures so that any error could be absorbed in T given 
in equation (5.28). By applying a least-square power fit to the 
four pairs of data, the values of T and V were found to be 

T = 3.6891 

(5.30) 

V = 0.48669 % 0.49 (5.31) 

The coefficient of determination for these results was 0.9980, which 
is excellent. Substituting equations (5.30) and (5.31) into equa- 
tion (5.28) gives the values of shown in table 5.4. The per- 

centage difference (D^ shown in table 5.4 varies by only 2 percent 
at all times. Therefore, from equations (5.28) and (5.31) the effect 
of the dimensionless material parameter on the dimensionless film 
thickness can be written as 

H a r 0 * 49 

min G (5.32) 

5 - 6 Minimum-Film-Thickness Formula 

The proportionality expressions (5.16), (5.22), (5.27), and 
(5.32) establish how the minimum film thickness varies with the el- 
lipticity, speed, load, and material parameters, respectively. This 
enables a composite minimum-f ilm-thickness formula for a fully flood- 
ed, isothermal, elastohydrodynamic point contact to be written as 


H 


= 3.63 U 0 , 68 G 0 , 49 W ' 0,073 (1 - e -0 ' 68 k , 


min " J * DJ u u w U - e ", ( 5 . 33 ) 

In equation (5.33) the constant 3.63 is different from that in equa- 


tion (5.30) to account for rounding off the material-parameter ex- 
ponent . 

T *bl* 5 ’ 5 8lVes 34 Afferent cases used in obtaining agna- 
tion (5.33). In this table. corresponds to the minimum film 

thickness obtained from the EHL point-contact theory developed in 
pter 4, and H mln is the minimum film thickness obtained from 
equation (5.33). The percentage difference between these two values 
expressed by Dj, which is defined in equation (5.17). i n ta- 
ble 5.5 the values of D 1 are within ±5 percent. 

It is sometimes more convenient to express the side-leakage 

factor in equation (5.33) in terms of the radius of curvature ratio 

W lnStMd ° f the eUiPticity parameter (k) through the follow- 
ing relationship: 


/r ) 0 ' 64 
k - 1.03 If 


(5.34) 


where 


J- = JL + JL 

R y r Ay r By 

-L = -i- + JL 

R x r Ax r Bx 


(5.35 


(5.3e 


A least square power fit was used in obtaining equation (5.34). This 

equation is valid for R y /R x between 1 and 40. The coefficient of 

determination associated with this power fit was 0.9997, which is 
excellent. 

Jsing equation (5.34) avoids the need to evaluate elliptic inte- 
grals of the first and second kind in the determination of k. The 
minimum film thickness can thus be derived directly from a knowledge 

of the radii of curvature of the contacting bodies (r r 

, % Ax* r Bx* r Ay* 

and V- 

It is interesting to 


compare the new point-contact, minimum- 
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* thickness formula (eq. (5.33)) with the corresponding equation 
generated by Dowson (1968) for line contacts 


H =2.65 "0*70 g 0.54 w 0-13 

min,L 


(5.37) 


The powers of U, G, and W in equations (5.33) and (5.37) are quite 
similar considering the different numerical procedures upon which 
they are based. It is also worth noting that the power (W) in 
equation (5.33) is extremely close to the value of -0.074 proposed 
by Archard and Cowking (1966) in their study of point contacts. 

5. 7 Central-Film-Thickness Formula 

There is interest in knowing the central film thickness, in 
addition to the minimum film thickness, in elastohydrodynamic con- 
tacts. In this section a central film thickness will be formulated. 
The procedure used in obtaining the central film thickness is the 
same as that used in obtaining the minimum film thickness. 

Table 5.6 gives 10 values of the ellipticity parameter (k) and 

the corresponding central film thickness (H ) as obtained from the 

c 

EHL point-contact theory developed in chapter 4. The other dimen- 
sionless parameters (0, W, and G) were held constant as defined by 
the relationship (5.10). Having these 10 pairs of data, the objective 
is to determine an equation that describes how the ellipticity param- 
eter affects the central film thickness. The general form of this 
equation can be written as 


(5.38) 


A least-square exponential curve fit to the 10 pairs of data points 


1 ■ 1 , • . . , 10 


k i’ I 1 “ H 


was used in obtaining values of I and J in equation (5.38). The 
value of H in equation (5.38) was found to be 8.69xl0 -6 and was 
determined by finding the value of H that gives a coefficient of 

C , L 


REPRODUCIBILITY OF THE 
ORIGIN AL PAG j is poor 


■ 
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determination (r 2 ) closest to 1. The values of t and J in equa- 
tion (5.38) and the coefficient of determination (r 2 ) as obtained 
from the least-square fit are 


I = 0.61178 % 0.61 
J = -0.73377 % -0.73 
r 2 = 0.96356 


(5.39) 

(5.40) 

(5.41) 


The value of the coefficient of determination (r 2 ) in equation (5.41) 
is not as good as that obtained for the ellipticlty portion of the 
minimum-film- thickness formula, which was 0.9990. From equations 
(5.38), (5.39), and (5.40) the following proportionality shows the 
effect of the ellipticlty parameter on central film thickness: 


H c “ (1 - 0.61 e 0,73 k ) 


(5.42) 


Substituting equations (5.39) and (5.40) into equation (5.38) gives 
H c , the dimensionless central film thickness obtained from the least- 
square formulation. The values of for the 10 values of ellip- 

ticity parameter (k) are given in table 5.6. The percentage differ- 
ence between the central film thickness obtained from the EHL 
point-contact theory (H c ) and the central film thickness obtained 
from the least-square fit (H c ) is expressed as 


H - H 

c c 


(5.43) 


Note in table 5.6 that is within the range -11% < ^ • 

This is a larger range than that found for the minimum-film-thickness 
formula, which was f 3 percent. 

Table 5.7 gives the dimensionless speed parameter (U) and the 
corresponding central film thickness as obtained from the F.HL point- 
contact theory developed in chapter 4. The values of the other di- 
mensionless parameters (k, W, and C.) were held constant as defined 
in equation (5.18). The general form of the equation that describes 
how the dimensionless speed affects the central film thickness can 
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be written as 


H = KU L 
c 


(5.44) 


By applying a least-square power fit to the 15 pairs of data 

{(U i' 1 " !» • • •, 15), the values of K and L were found 

to be 

i '* 537 '" (5.45) 

L - 0.66722 % 0.67 (5.46) 

The coefficient of determination (r 2 ) f or these resnits was 0.9990. 

Substituting equations (5.45) and (5.46) into equation (5.44) glues 

the values of H c shown in table 5.7. The percentage difference 

(0 2 > between the central film thickness obtained from the EHL point- 

contact theory (H.) and the central film thickness obtained from the 

least-square fit (iy is expressed in equation (5.43) and given in 

table 5.7. In table 5.7, D, is In the range -4% < D., < 61. This 

is to be compared to the range of Dj for the minimum film thick- 
ness, which was ±2 percent. 

From equations (5.44) and (5.46) the effect of dimensionless 
speed on dimensionless central film thickness can be written as 

H cr Jj0 . 67 

c U (5.47) 

Table 5.8 gives the dimensionless load parameter (U) and the 

corresponding central film thickness Oy a. obtained from the EHL 

point-contact theory developed in chapter 4. The values of the other 

dimensionless parameters (k, u. and a) were held constant as defined 

equation (5.23). The general form of the equation that describes 

how the dimensionless load affects the central film thickness can 
be written as 


H c - MW 


,N 


(5.48) 


By applying a least-square power fit to the eight pairs of data 
((Wf, H^j), i * 1, . . . t g} > the va j ues Q f ^ ;ln< j 


were found 
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to be 


M = 2.8508xl0" 6 

(5.49) 

N - -0.067248 % -0.067 (5 5Q 

The coefficient of determination (r*, for these results was 0.7303 
which is „„ t as good a fit as the 0.0260 obtained for the load por- 
tion of the minimum- film- thickness formula. Substituting equations 
(5.6,) and (5.50) into e, cation (5.48) gJv e s the voices of „ c shoOT1 

in table 5.8. The percentage difference (D ^ . 

'^2' between the central film 

thickness obtained from the EHL point-contact theory fly a „d the 

central film thickness obtained fro, the least-square fit („ c , , s 

pressed in elation (5.43) and given in table 5.8. I„ tablets, D 
is ±5 percent. 

From equations (5.48, and (5.50) the effect of dimensionless 
loed on dimensionless central film thickness can be written as 

H or w ~0.067 

c (5.51) 

Table 5.9 gives the four material parameters and the correspond- 
ing central film thickness m i , 

thickness (h c > as obtained tnm the polnt . cont „ ct 

theory developed in chapter 4. The ellipticity parameter (k, was 

held fixed at 6 for the results shown in table 5.9. The general 

rra of the equation that describes how the dimensionless material 

parameter (c, affects the central film thickness can he written as 

rmV 


Q = Tr; 


where 


(5.52) 


H 


(1 - 0.61 e' U,?3 k ;U 0 -67 w -0.067 

By applying a least-square power fit to the four pairs of data, the 
values of T and V were found to he 

T "2.6595 

(5. 53) 

V « 0.52776 % 0.53 

(5.54) 

The coefficient of determination (r 2 ) for n 

' r ) tor these results was 0.4804. 
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Substituting equations (5.VJ) and (5.54) into equation (5.52) gives 

the values of H shown in table 5.8. The values of as ex- 

^ 2 

pressed in equation (5.43) are shown in table 5.9. In table 5.9, the 
values of are within the range of -4% < - 5%. 

From equations (5.52) and (5.54) the effect of dimensionless 
material parameter on the dimensionless central film thickness can 


be written as 

H c - C 0 * 53 (5.55) 

The proportionality expressions (5.42), (5.47), (5.51), and 
(5.55) establish how the central film thickness varies with the el- 
lipticity, speed, load, and material parameters, respectively. This 
enables a composite central-f iim-thickness formula for a fully flooded 
isothermal, elastohydrodynamic point contact to be written as 

H c = 2.69 u 0 * 67 g °- 53 W _0 ' 067 (1 - O.ol e~°’ 73 k ) (5.56) 

If it is desired, the side- leakage factor in equation (5.56) can be 

expressed in terms of the radius-of-curvature ratio (R /R ) instead 

y x 

of the elliptic ity parameter (k) by using equation (5.34). 

Comparing the central-film-thickness formula (5.56) with the 
mi nimum- f ilm- thickness formula (5.33) reveals a slight difference. 

In equation (5.56) the load exponent is small but negative, as it. was 
for the minimum- f i lm- thickness formula. This is in contrast witt; 
the recent numerical study of Ranger et al. (1975), who found a small 
but positive exponent on the dimensionless load (W) in their formula- 
tion of a central film thickness. 

Table 5.10 gives the different cases used to obtain equa- 
tion (5.56). In this table, correspond to the central film 

thickness obtained from the I'iiL point-contact theory ieve loped in 
chapter 4, and corresponds to the central film thickness ob- 

tained from equation (5.5e>). The percentage difference between these 
two values is expressed by !>•>, which is defined in equation (S.4i). 




In table r .10 the values of are within _ L 10 percent. 

5. 8 Concluding Remarks 

By using the procedures outlined in an earlier chapter the in- 
fluence of the ellipticity parameter and the dimensionless speed, 
load, and material parameters on minimum film thickness have been 
investigated. The ellipticity parameter was varied from 1 (a ball- 
on-plate configuration) to 8 (a configuration approaching a line 
contact^ . The dimensionless speed parameter was varied over a range 
of nearly two orders of magnitude. The dimensionless load parameter 
was varied over a range of about one order of magnitude. Situations 
equivalent to using solid materials of bronze, steel, and silicon 
nitride, and lubricants of paraffinic and napthenic mineral oils were 
considered in an investigation of the role of the dimensionless ma- 
terial parameter. Thirty-four different cases were sed to generate 
the following minimum- film- thickness and central-film-thickness 
relationships: 


n . - 3.63 u 0 - b 8 c. 0 - 4 V [> - 0n (i - «-°' 68 k ) 

mm 

5 - 2.69 „0.67 c .0.53 w -U. 067 (1 . Q . 61 ,-0.73 k, 
c 


The ellipticity parameter (k) car be written as 


a \°- M 

k - 1.031 


Contour plots have been presented that indicate in detail the 
pressure distribution and the film thickness. In some solutions, 
pressure spikes were in evidence. The theoretical solutions of film 
thickness have all the essential features of previously reported ex- 
perimental observations based upon optical interferometry. 

The importance of the present chapter lies in the fact that it 
presents f c r the first time a satisfactory theoretical film- 
hickness equation for elas tchydrodynamic point contacts operating 
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under fully flooded conditions. The exponents on the various dimen- 
sionless parameters governing minimum film thickness in such conjunc- 
tions are quite similar to those developed by Dowson (1968) for 
line contacts. The most dominant exponent c-curs In associ-.tlon 

” Uh th " SP ' ed para " eter ' “Mm the exponent on the load parameter 
Is very small and negative. The material parameter also carries a 
Significant exponent, although the range of this parameter in en- 
gineering situations is limited. Ranger, et al. (1975) have de- 
veloped a central-film-thickness formula for the contact geometry 
of a ball on a plate from which an estimate can be made of the 
minimum film thickness. However, the work presented in this chap- 
ter is valid for any contact geometry and proceeds directly to the 
evaluation of the minimum film thickness. 

Perhaps the most significant feature of the prised minimum- 
fi.lm- thickness formula is that it can be applied to any contacting 
solids that present an elliptical Hertaian contact region. Many ma- 
chine elements, particularly rolling-element bearings, possess such 
geometry, and it is expected that the new mi„i»„„-f ilm-thickness 
equation will find application in such fields. 
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. CHAPTER 6 
STARVATION RESULTS 

It was not until the late 1960's and early 1970's that the in- 
fluence of lubricant starvation upon elastohydrodynamic behavior 
received serious consideration. Prior to this time it was assumed 
that the inlets were fully flooded. This assumption seemed to be 
entirely reasonable in view of the minute quantities of lubricant 
required to provide an adequate film. However, in due course it was 
recognised that some machine elements suffered from lubricant star- 
vation. 

How partial filling of the inlet to an elastohydrodynamic con- 
junction influences pressure and film thickness can readily be ex- 
plored theoretically by adopting different starting points for the 
inlet pressure boundary. Orcutt and Cheng (1966) appear to have 
been the first to proceed in this way for a specific case corre- 
sponding to a particular experimental situation. Their results 
showed that lubricant starvation lessened the film thickness. 
Wolveridge, et al. (1971) used a Grubin (1949) approach in an 
analysis of starved elastohydrodynamic lubricated line contacts. 
Wedeven, et al. (1971) analyzed a starved condition in a ball-on- 
plate geometry, and Castle and Dowson (1972) presented a range of 
numerical solutions for the starved line-contact elastohydrodynamic 
situation. In these references the analysis yielded values of the 
proportional reduction in centerline film thickness from the fully 
flooded condition in terms of a dimensionless inlet boundary pa- 
rameter. 

In the present chapter, 15 cases in addition to Litre e pro- 
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sented in chapter 5 were used to obtain the starvation results. 

From the results a simple dimensionless inlet boundary distance was 
written. This inlet boundary distance defines whether a fully flooded 
or a starved condition exists in the contact. Furthermore, it was 
found that the film thickness for a starved condition could be 
written in dimensionless terms as a function of the film thickness 
for a fully flooded condition and the inlet distance parameter. Con- 
tour plots of pressure and film thickness in and around the contact 
are shown for fully flooded and starved conditions. The theoretical 
findings are compared directly with previously reported results. 

6.1 Boundary between Fully Flooded and Starved Conditions 

Figure 6.1 shows the computing area in and arou*'^ the Hertzian 
contact. In this figure, as defined in chapter 4, the coordinate X 
is made dimensionless with respect to the semiminor axis (b) of the 
contact ellipse and the coordinate Y is made dimensionless with re- 
spect to the semimajor axis (a) of the contact ellipse. The elliptic- 
ity parameter (k) is defined as the semimajor axis divided by the 
semiminor axis of the contact ellipse (k = ./b) . Because of the di- 

mensionless form of the coordinates X and Y, the Hertzian contact 
ellipse becomes a Hertzian circle regardless of the ellipticity pa- 
rameter. This Hertzian contact circle is snown in figure 6.1 with a 
radius of unity. The edges of the computing area, where the pressure 
is assumed to be ambient, are also denoted. In this figure the di- 
mensionless inlet distance (m) , which is equal to the dimensionless 
distance from the center of contact to the inlet edge of the com- 
puting area, is shown. 

Lubricant starvation can be studied simply by reducing the di- 
mensionless inlet distance (m) . A fully flooded condition is said to 
exist when the dimensionless inlet distance ceases to influence in 
any significant way the minimum film thickness. When starting from 
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a fully flooded condition and decreasing m, the value at which the 
minimum film thickness first starts to change is called the fully 
flooded - starved boundary and is denoted by m*. Therefore, lubri- 
cant starvation was studied by using the basic elastohydrodynamic 
lubrication point-contact theory developed in chapter 4 and observ- 
ing the effect of reducing the dimensionless inlet distance. 

Table 6.1 shows the effect of changing the dimensionless inlet 
distance upon the dimensionless minimum film thickness for three 
groups of dimensionless load and speed parameters. For all the re- 
sults presented in this chapter the material parameter (G) is fixed 
at 4522 and the ellipticity parameter is fixed at 6. In this table 


it is seen that as the dimensionless inlet distance (m) decreases the 


dimensionless minimum film thickness (H , ) decreases. 

' min 

Table 6.2 shows how the three groups of dimensionless speed 
and load parameters affect the location of the dimensionless inlet 
boundary distance (m*) . Also given in this table are the correspond- 
ing values of dimensionless central and minimum film thickness for 
the fully flooded condition as obtained by interpolation of the nu- 
merical values. The value of the dimensionless inlet boundary (m*) 
shown in table 6.2 was obtained by using the data from table 6.1 when 
the following equation was satisfied: 


^minjF ( H rain) 


m=m* 


H . ^ 

min, F 


= 0.03 


( 6 . 1 ) 


The value of 0.03 is used in equation (6.1) since it was ascertained 
that the data in table 6.1 could only be obtained to an accuracy of 
±3 percent. 

The general form of the equation that describes how the dimen- 


sionless inlet distance at the full^ flooded - starved boundary (m*) 
varies is given as 
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m* - 



( 6 . 2 ) 


The right side of equation (6.2) is similar to the forms of the 
equation given by Wolveridge, et al. (1971) and Wedeven, et al. 
(1971). By applying a least-square power fit to the data obtained 
from table 6.1, the following can be written: 


m* = 1 + 3.06 



-.0.58 



(6.3) 


A fully flooded condition exists when m _> m*. and a starved condi- 
tion exists when m < m*. The coefficient of determination (r^) for 


these results is 0.9902, which is excellent. 


If in equation (6.2) the dimensionless minimum film thickness 
is used instead of the central film thickness, the following is 
obtained : 


m* = 1 + 3.34 



10.56 


H 


min,F 


(6.4) 


The coefficient of determination for these results is 0.9869, which 


again is excellent. 

From Wedeven, et al. (1971), using the symbols of this thesis, 
the dimensionless inlet distance at the fully flooded — starved 
boundary can be written as: 


“w 



12/3 



(6.5) 


Comparing equation (6.3) with equation (6.5) indicates close agree 
ment with Wedeven, et al. (1971). The latter, however, predicts a 
slightly higher value of the fully flooded - starved boundary than 
predicted from the present results. 
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6.2 Starvation Film Thickness Formulas 

Having clearly established the limiting location of the inlet 
boundary for the fully flooded conditions (eqs. (6.3) and (6.4)), 
an equation defining the dimensionless film thickness for lubricant 
starvation conditions will be developed. The relationship between 
the dimensionless central film thickness in starved and fully 
flooded conditions can be expressed in general form as 


H 

H 


c , S = 
c,F 



( 6 . 6 ) 


Table 6.3 shows how the ratio of the dimensionless inlet distance 

parameter to the fully flooded - starved boundary [ (m - l)/(m* - 1) ] 

affects the ratio of central film thickness in the starved and fully 

flooded conditions (H /H ). A least-square power curve fit to 

c , b c , r 

the 16 pairs of data points 


7 H c,s\ 

(m - 1 

'll 

Ik,)/ 

\m* ~ 1 

4J 


was used in obtaining values for 
For these values of C* and D* 
thickness for a starved condition 


i = 1,2, . . ., 16 

C* and D* in equation (6.6). 
the dimensionless central film 
can be written as 


H 


c,S 



(6.7) 


By using a similar approach while making use of the data in table 6.3 % 
the dimensionless minimum film thickness for a starved condition can 
be written as 


H . c = H , t 

mm,S min,F\m* 


(JLr-l ) 

, F\m* - 1/ 


0.25 


( 6 . 8 ) 


Therefore, whenever m < m*, where m* Is defined by either equa- 
tion (6.3) or (6.4), a lubricant starvation condition exists. When 


this is true, the dimensionless central film thickness is expressed 

by equation (6.6) and the dimensionless minimum film thickness is ex 
pressed by equation (6-8). if m > m*. where m* ls defined by 

either equation (6.3) or (6.4), a fully flooded condition exists. 

The dimensionless central and minimum film thicknesses for a fully 
flooded condition (H^ and H^) were developed in chapter 5 
and are expressed in equations (5.56) and (5.33), respectively. 

(That 1S ’ H min in e <l* (5.33) is equivalent to p in e q. (6.8), 

and H c -in equation (5.56) is equivalent to in eq . (6 . 7)<) 

The ratio of dimensionless inlet distance to the fully flooded - 
starved boundary as obtained from Wedeven, et al. (1971), expressed 
On - D/Objj - D , is also given in table 6.3. By comparing 
these results with the results obtained from the present thesis 

[(m - l)/(m>. ~ U], it can be seen that for group 1 the agreement 
is excellent. However, the agreement in groups 2 and 3 is not as 
good. A possible explanation for this difference can be than an 
approximate expression for the Hertzian deformation is used in the 
Wedeven, et al. (1971) analysis. They indicate their equation 
(eq. (6.5)) is only valid for small m* or more specifically 
m * < 3 ’ group 2, m* = 3.71 and in group 3 m* = 5.57. Since 
no such assumption is required in deriving equations (6.3) and (6.4), 
they would seem to be more appealing. 

Figure 6.2 shove the influence of Infer boundary parameter upon 
central film thickness for the Wedeven. et al. (1971) results and 
those obtained from the present thesis. Prom this fi c „re ,t is ob- 
served that the wedeven, et al. (1,71, results s lve sl, s htl, higher 
values of the eentral film thickness ratio of starved to fuily flooded 
condition than Lhos obtained from the present results. 

6 . 3 Contour Plot of Rosuj t s 

To explain more fully what happens in 


going from a fully flooded 
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to a lubricant starvation condition, figures 6.3 to 6,19 are pre- 
sented. As in chapter 5, in parts (a) of these figures, contour 
plots of dimensionless pressure are given; and in parts (b) of these 
figures, contour plots of dimensionless film thickness are given. 

In parts (a) and (b) the + symbol indicates the center of the 
Hertzian contact. The Hertzian contact circle is shown in each of 
parts (a) and (b) by asterisks. At the top of each of parts (a) 
and (b) , the contour labels and the corresponding values are given. 

In figures 6.3(a), 6.4(a), . . 6.7(a), contour plots of di- 

mensionless pressure (P = p/E 1 ) are given for group 1 of tabie 6.2 
and for dimensionless inlet distances (m) of 4, 3, 2, 1.5, and 
1.25, respectively. The contour labels and values are kept constant 
in going from figure 6.3(a) to 6.7(a). In figure 6.3(a) a fully 
flooded condition exists. Once starvation occurs the severity of the 
situation increases as m is decreased, thus implying that the most 
severe starvation case is shown in figure 6.7, where m - J.25. In 
figures 6.3(a), 6.4(a), and 6.5(a) a pressure spike is clearly vis- 
ible, whereas in figures 6.6(a) a.'.d 6.7(a) no pressure spike is 
present. Note in figure 6.7(a), the most severe starvation case, 
that the contour labeled H does not extend as far to the left as 
it did for the fully flooded pressure results shown in figure 6.3(a). 

In figures 6.3(b), 6.4(b), . . ., 6.7(b), the contour plots of 
dimensionless film thickness (H = h/R^) are given for group 1 of ta- 
ble 6.2 and m of 4, 3, 2, 1.5, and 1.25, respectively. The film 
thickness results shown In figure 6.3(b), . . 6.7(b) correspond 

to the pressure results shown In figures 6.3(a), . . ., 6.7(a). The 
central portion of the film thickness contours has become more paral- 
lel as starvation has increased (n; decreasing) with the minimum film 
thickness decreasing. Note ais< that the film thickness contour 
values for the most severely starved condition (fig. 6.7(b)) are much 
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lower than the film thickness contour values for the fully flooded 
conditions (fig. 6.3(b)). 

In figures 6.8(a), 6.9(a), . . and 6.12(a), contour plots of 
dimensionless pressure (P = p/E 1 ) are given for group 2 of table 6.2 
and for m of 6, 4, 3, 2.5, and 2, respectively. The contour values 
are the same for each of these figures. In figure 6.8(a) (the fully 
flooded conditions) the contour for the largest pressure (contour A) 
is present, but as starvation occurs in figures 6.9(a), 6.10(a), 
and 6.11(a) this contour is absent. Furthermore, for the severely 
starved condition shown in figure 6.12(a) both the A and B contours 
are absent. This implies that or group 2 the pressure peak present 
in a fully flooded condition gets flatter as starvation progresses - 

Figures 6.8(b), 6.9(b), . . 6.12(b) give contour plots of film 

thickness for group 2 of table 6.2 and m of 6, 4, 3, 2.5, and 2, 
respectively. The minimum film thickness areas in these figures occur 
on the axial center of the contact and move slightly to the right as 
starvation becomes severe. Note from the contour values that film 
thickness decreases substantially in going from a fully flooded con- 
dition (fig. 6.8(b)) to a severely starved condition (fig. 6.12(b)). 

As was found for the low-speed results, the central portions of the 
film thickness contours become parallel as starvation is increased. 

In figures 6.13(a), 6.14(a), ...» 6.i9(a), contour plots of 
dimensionless pressure (P p/E f ) are given for group 3 of table 6.2 
and for m of 6, 4, 3, 2.5, 2, 1.75, and 1.5, respectively. The 
contour values are the same for each of these figures. Figure 6.13(a) 
gives the pressure contour for a fully flooded condition; figure 
6.19(a) gives the pressure contour for a severe i,: starved condition. 

As starvation occurs the pressure profile flattens, since in figure 
6.17(a) contour A is absent, in figure 6.18(a) contours A and B are 
absent, and in figure 6.19(a) (the severely starved condition) eon- 
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tours A, B, and C are absent. Also from these figures it is found 
that the distance from the center of the contact to the upstream 
location of the largest contour value (labeled H) decreases as the 
severity of lubricant starvation increases. 

Figures 6.13(b), 6.14(b), . . ., 6.19(b) give contour plots of 
film thickness for group 3 of table 6.2 and m of 6, 4, 3, 2.5, 2, 
1.75, and 1.5, respectively. In figure 6.13(b) (the fully flooded 
condition) the minimum film thickness occurs to the sides of the con- 
junction in two areas that are midway between the center of the con- 
tact and the Hertzian circle. As m is decreased or the severity of 
starvation increases (going from figs. 6.14(b) to 6.19(b)) the mini- 
mum film thickness area remains in the axial center of the conjunction 
but moves to the right, nearer the Hertzian circle. Note the simi- 
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larity among the film thickness contours of figure 6.7(b) (group 1 
case), figure 6.12(b) (group 2 case), and figure 6.19(b) (group 3 
case) . 

The dimensionless pressure (P = p/E*) on the X-axis is shown 
for three values of dimensionless inlet distance and for groups 
1 and 3 of table 6.2, respectively, in figures 6.20(a) and (b) . 

The value of Y is held constant near the axis of symmetry of the 
conjunction. in these figures as a conjunction becomes starved (as 
m is decreased) the pressure spike diminishes. 

Figures 6.21(a) and (b) show the dimensionless film thickness 


(H * h/R^) cm the X-axis for three values of dimensionless inlet 


distance and for groups 1 and 3 of table 6.2, respectively. The 
value of Y is held fixed close to the axis of symmetry of the con- 
tact. In these f-gures, particularly figure 6.21(b). the central 
region becomes flatter as starvation occurs. Also, in going from a 
fully flooded condition to a starved condition the film thickness 
decreases substantially. 
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6.4 Concluding Remarks 

By using the theory and numerical procedure outlined in chap- 
ter 4, the influence of lubricant starvation upon minimum film thick- 
ness in starved elliptical elas tohydrodynamic conjunctions has been 
investigated* This study of lubricant starvation was performed by 
moving the inlet boundary closer to the center of the conjunction. 
From the results it was found that the location of the dimensionless 
inlet boundary (m*) between fully flooded and starved conditions 
could be expressed simply as 


or 



n0.58 



m* = 1 + 3.34 



10.56 


H 


min,F 


That is, for a dimensionless inlet distance (m) less than m* , star- 
vation occurs and for m > m*, a fully flooded condition exists. 
Furthermore, it has been possible to express the central and minimum 
film thicknesses for a starved condition as 


H 


c , S 


H 

c 


(js-z 

, F\m* - 


t) 


0.29 


H 


min, S 


H 


min , F 


\m* - 1/ 


where 

H fully flooded dimensionless central film thickness 

c , r 

H . „ fully flooded dimensionless minimum film thickness 

min, F 

m dimensionless inl^f distance 

m* dimensionless inlet distance at the fully flooded - starved 

boundary 

Contour diagrams of t ho pressure and film thickness in and 
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around the contact have been presented for both fully flooded and 
starved conditions. It Is evident from the contour diagrams that 
the pressure spike becomes suppressed and the film thickness de- 
creases substantially as the severity if starvation increases. 

The results presented in this chapter, when combined with the 
findings of the previous chapter, enable the esse:.cial features of 
starved, elliptical, elastohydrodynamic conjunctions to be ascer- 
tained. 
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CHAPTER 7 

SUMMARY OF CONCLUSIONS 

A procedure for the numerical solution of the complete, iso- 
thermal, elastohydrodynamic lubrication problem for point contacts was 
given. This procedure calls for the simultaneous solution of the 
elasticity and Reynolds equations. In the elasticity analysis the 
contact zone was divided into equal rectangular areas. It was as- 
sumed that a uniform pressure was applied over each area. In the nu- 

3/2 

merical analysis of the Reynolds equation the parameter <p - PH , 
where P is dimensionless pressure and H is dimensionless film 
thickness, was introduced to help the relaxation process. The 
pressure-viscosity analysis of Roelands (1966) was used. The numer- 
ical co ling of the elasticity and Reynolds equations results In a 
converged solution for the pressure profile. This pressure profile is 
then integrated over the computing zone to give the value of the cor- 
responding normal applied load. This load is then compared with the 
input lc id and corrections are made to the film thickness until these 
two loads are in agreement. 

The most important practical aspect of the elastohydrodynamic 
lubricated point -contact theory developed is the determination it the 
minimum film thickness within the contact. That is, the maintenance 
of a fluid film of adequate magnitude is extremely important to the 
operation of some machine elements. The minimum film thickness for a 
fully flooded conjunction was frond to be a function of the ellipticity 
parameter and the dimensionless speed, load, and material parameters. 

In the results the ellipticity parameter was varied from I (a ball -on- 
plate configuration) to 8 (a cont igurat ion approaching a line contact). 
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The dimensionless sneed parameter was varied over a range of nearly 
two orders of magnitude. The dimensionless load parameter was varied 
over a range of one order of magnitude. Situations equivalent to the 
use of solid materials of bronze, steel, and silicon nitride and lubri- 
cants of paraffinic and naphthenic mineral oils were considered in an 
investigation of the role of the dimensionless material parameter. 
Thirty-four different cases were used to generate the minimum film 
thickness and central film thickness formulas given below as 

= 3.63 U°* 68 G 0 - 4 V°' 073 (i - e' 0 * 68 k ) 


H 


min,F 

" V— W— (i - Q.61 * 
where the dimensionless speed parameter 


f 0. 67 0. 53..-0. 067 , 


-0.73 k 


) 


(7.1) 

(7.2) 


U = 


V 

E'R 


the dimensionless load parameter 


W = 


E'R 2 

x 


(7.3) 


(7.4) 


the dimensionless material parameter 

g = -*1_ 

p 

iv,as 

and the dimensionless ellipticity parameter 


k - * 
b 


(7.5) 


(7.6) 


Equations (7.5) and (7.6) can also be written in more convenient form 


as 


G = E'a 


.0.64 


k - 1.03 U 

V x> 


(7.7) 


(7.8) 


In equation (7.8) the ellipticity parameter is expressed strictly in 

terms of the radii of curvature and thereby eliminates the common 

practice of evaluating the elliptical integrals of the first and 
second kind* 



The Importance of equation (7.1) lies in the fact that it pre- 
sents for the first time a satisfactory theoretical film-thickness 
equation for elastohydrodynamic point contacts operating under fully 
flooded conditions. The exponents on the various dimensionless pa- 
rameters governing minimum film thickness in such conjunctions are 
quite similar to those developed by Dowson (1968) for line contacts. 
The most dominant exponent occurs in association with the speed pa- 
rameter, while the exponent on the load parameter is very small and 
negative. The material parameter also carries a significant expo- 
nent, although the range of this parameter in engineering situations 
is limited. Ranger, et al. (1975) have developed a central-film- 
thickness formula for the contact geometry of a ball on a plate from 
which an estimate can be made of the minimum film thickness. How- 
ever, the work presented in this chapter is valid for any contact 

geometry and proceeds directly to the evaluation of the minimum film 
thickness . 

Perhaps the most significant feature of the proposed minimum- 
film- thickness formula is that it can be applied to any contacting 
solids that present an elliptical Hertzian contact region. f if .y ma- 
chine elements, particularly rolling-element bearings, possess such 
geometry, and it is expected that the new minimum-film-thickness 
equation will find application in such fields. 

Contour plots of the fully flooded results have been presented 
that indicate in detail the pressure distribution and the film thick- 
ness. In some solutions, pressure spikes were in evidence. The 
theoretical solutions of film thickness have all the essential fea- 
tures of the previously reported experimental observations based upon 
optical interferometry. 

In addition to the fully flooded studies, the influence of lubri- 
cant starvation upon minimum film thickness in starved elliptical 
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elastohydrodynamic conjunctions has been Investigated. This study of 
lubricant starvation was performed by moving the inlet boundary closer 
to the center of the conjunction. From the results it was found that 
the location of the dimensionless inlet boundary (m*) between fully 
flooded and starved conditions could be expressed simply as 

lO. 58 

(7.9) 


m* = 1 + 3.06 


«c, P 


or 


m* = 1 + 3.34 



-|0.56 


H 


min,F 


(7.10) 


That is, for a dimensionless inlet distance (m) less than m*, starva 
tion occurs and for m > m*, a fully flooded condition exists. Fur- 
thermore, it has been possible to express the dimensionless central 
and minimum film thicknesses for a starved condition as 




v0.29 


H „ 
c , S 

- H / m ■ i 
c,F\m* - 1, 

) 

(7.11) 

min,S 

= H ( m - 

min,F\m* - 

k)°- 25 

(7.12) 


Contour diagrams of the pressure and film thickness in and around 
the conjunction have been presented for starved conditions. It is 
evident from the contour diagrams that the pressure spike becomes sup- 
pressed and the film thickness decreases substantially as the severity 
of starvation increases. 

The starvation results when combined with the fully flooded re- 
sults enable the essential features of starved, elliptical, elasto- 
hydrodynamic conjunctions to be ascertained. 


82 


REFERENCES 

Archard, J. F., and Cowking, E. W. (1965-66) Elastohydrodynamic lubri- 
cation of point contacts , Proc. Inst. Mech. Engrs. London, vol. 180. 
pt. 38, p. 47. 

Archard, J. F. , and Kirk, M. T. (1964) Film thickness for a range of 
lubricants under severe stress, J. Mech. Eng. Sci .. vol. 6, p. 101. 

Barus, C. (1893) Isotherma, isopiestics, and isometrics relative to 
viscosity, Amer. J. Sci .. vol. 45, p. 87. 

Blok, H. (1965) Proceedings of the International Symposium on Lubri ca- 
tion and Wear. D. Muster and B. Sternlicht, eds . , McCutchan, Berkeley, 

p. 1. 


Cameron, A., and Gohar, R. (1966) Theoretical and experimental studies 
of the oil film in lubricated point contact, Proc. Roy. Soc. (London), 
vol. 291A, p. 520. L 

Castle, P., and Dowson, D. (1972) A theoretical analysis of the starved 
elastohydrodynamic lubrication problem for cylinders in line contact, 
Elastohydrodynamic Lubrication: 1972 Symposium, Inst. Mech. Engrs 

(London), p. 131. 

Cheng, H. S. (1970) A numerical solution of the elastohydrodynamic film 
thickness in an elliptical contact, J, Lubric, Tech ., vol. 92, p. 155. 

Christensen, H. (1964) The variation of film thickness in highly loaded 
contacts, ASLE Trans . . vol. 7, p. 219. 

Dowson, D. (1965) Elastohydrodynamic lubrication - an introduction and 
a review of theoretical studies. Paper Rl, Inst. Mech. Engrs., 

(London), May 1965. 

Dowson, D. (1968) Elastohydrodynamic, Proc. Inst. Mech. Engrs. (London), 
vol. 182, pt. 3A, p. 151. L 

Dowson, D., and Higginson, G. R. (1959) Numerical solution to the elas- 
tohydrodynamic problem, J. Mech. Eng. Sci .. vol. 1, p. 6. 

Dowson, D., and Higginson, G. R. (1961) New roller-bearing lubrication 
formula, Eng, (London) , vol. 192, p. 158. 

Dowson, D., and Higginson, G. R. (1966) Elas tohydrodynamic Lubrication. 
Pergamon Press, New York. ~ " ~ * 

Dowson, D. , and Whitaker, A. V. (1965-66) A numerical procedure for the 
solution of the elastohydrodynamic problems of rolling and sliding con- 
tacts lubricated by a Newtonian fluid, Proc. Inst. Mech. Engrs. (London} 
vol. 180, pt. 3B, p. 57. U 

Dyson, A., Naylor, H., and Wilson, A. R. (1965-66) The measurements of 
oil-film thickness in elastohydrodynamic contacts, Proc. Inst. Mech. 
Engrs. (London) , vol. 180, pt. 3B, p. 119. ’ 

Grubin, A. N., Vinogradaova, I. E., and Ketnva, F., eds. (1949) Investi- 
£ation of the Contact Machine Components . Central Sci. Res. inst. Tech. 
Mech. Eng., Book 30, (D.S.I.R. translation 337), Moscow. 


83 


Harris, T. A. (1966) Rolling Bearing Analysis, Wiley, New York. 

Hamrock, Bernard J., and Anderson, William J. (1973) Analysis of an 
arched-outer-race ball-bearing considering centrifugal forces, 

J. Lubric. Tech., vol. 95, p. 265. 

Hertz, H. (1882) The contact of rigid elastic bodies, J. fur die Reiner 
und Angew andte Mathematik . vol. 92, p. 156. — 

Martin, H. M. (1916) The lubrication of gear-teeth, Eng. (London), 
vol. 102, p. 119. — 

Orcutt, F, K. , and Cheng, H. S. (1966) Lubrication of rolling-contact 
instrument bearings, Gyro-Spin Axis Hydrodynamic Bearing Symposium , 
vol. 2, tab 5, Mass . Inst. Tech. Instrum. Lab., Cambridge. 

Pai, S. I. (1956) Viscous Flow Theory , vol. I, Van Nostrand, Princeton. 

Ranger, A. P., Ettles, C. M, M., and Cameron, A. (1975) The solution of 
the point contact elasto-hydrodynamic problem, Proc. Royal Soc. 

(London) , vol. 346, p. 227. 

Reynolds, 0. (1886) On the theory of lubrication and its application 
to Mr. Beauchamp Tower's experiments, including an experimental deter- 
mination of the viscosity of olive oil, Phil. Trans. Roy. Soc. (London), 
vol. 177, p. 157. 1 1 L 

Roelands, C. J. A. (1966) Correlational Aspects of the Viscosity- 
Tempera ture-Pressure Relationship of Lubricating Oils . Druk. U.R.B.. 
Groningen. 


Sibley, L. B., and Orcutt, F. K. (1961) Elastohydrodynamic lubrication 
of rolling-contact surfaces, ASLE Trans. . vol. 4, p. 235. 

Timoshenko, S«, and Goodier, J, H. (1951) Theory of Elasticity, 2nd ed., 
McGraw-Hill, New York. 

Wedeven, L. E., Evans, D., and Cameron A. (1971) Optical analysis of 
ball bearing starvation, J. Lubric. Tech ., vol. 93, p. 349. 

Wholmes, T. L. (1966) The Effect of Surface Quality on Lubricating Film 
Performance , Ph.D. Dissertation, University of Leeds. 

Wolveridge, P. E., Baglin, K. P., and Archard, J. G., (1971) The 
Starved Lubrication of Cylinders in Line Contact, Proc. Inst. Mech. 

Engrs . (London) , vol. 185, p. 1159. 






APPENDIX - LISTING OF COMPUTER PROGRAMS 


l* 
24 
34 
44 
5 * 
64 
74 
8 * 
94 
1 04 
114 
1?4 
1 34 
] 4* 

1 54 
164 
174 
1P4 
194 
20 4 

2 14 
??4 
234 


IMPLICIT DOUBLE PRECISION <A-H,0-Z> 

DIMENSION P?<95 t 20,3> 

DIME NS ION DENS<7 60 ),VlS(7t0 1,XMU(7fe0 1 f ZPRI7 60) 

DI HFNSION A ( 760 1 , R t 760 1 t C l 760 1 ,DLZl7t01*XL (7 60 1 ,XM ( 7601 
COMMON NSAVE , MB J f NX f NY 

CO MM ON A 1 # B 1 # P1 ,BHO,RX, RY ,XK # PI ,SA ,Sb,T 3, T? 

COMMON U ( 7foD I ,D< 1520 ) , $ ( 760 ) t H ( 1 620 » ,PR< 1520 > ,PHI < 76D > ,PRSV ( 760 ) 
COMMON EP ,HMI N,ZC,ZD, XN ,YM 
C 

C EL ASTOHYDRODYNAMIC LUBRICATION OF POINT CONTACT D R OCR AM 

C 

C ALL INPUT DATA IS IN THE UNITS OF NF U T ON t C E N T I M E T E R , Sf C ON 0 

C 

C 

C INPUT 

c 

DLROSMz 100.000 
PH IMX r 100.00-10 
MR jrO 
M B JO - 0 
MB J09 I 1 
MB JQ6 - 1 
MB JO 6 - 1 


? 4 $ 
754 
?64 

2 74 
r 84 
29 4 
304 
314 
7 74 

3 3 * 

3 44 
364 
3C4 
374 
384 
3 B 4 
4(14 

4 1 4 
4 ?4 
4 3 4 
4 44 
4 64 
4 6 4 
4 74 
4 P 4 
490 
5P4 
1 1 4 
674 
6 34 
*>4 4 
* 6 4 
564 
574 
684 
604 
6D4 
614 
674 
t 34 
684 
6 64 
f 64 

6 74 
6«4 
694 
7 n 4 
714 
774 

7 34 
744 
7‘4 
764 
774 
784 
794 
804 


M AUPrO 

jr NN :0 

R P I AN * i . uno 

772 - - 1 3 DO 

NPQ-n 

MH 9 I 0 

M H 0 i 1 

OEL: . 000000 inn 
» T: 3 . 14 1 S 9265359 DD 
ORE i 1 .600 
DEL 4 = 1 .0600 
DEL 1 : • 5 0 0 
ROSMi 1 DD .000 
C 

C CONVEX SURFACES are POSITIVE 

C CONTAVE SURFACES ARE NEGATIVE 

C 

RA V =1 • 1 1 I 7500 
R A X * 1 . 1 1 1 7 S 0 C 
R P X : 1 .0 1 7 
R R Y : -1 * 1861 IS 00 
X I A - • 300 
G A : 1 .9994 70 7 
P 1 VAS^M 868 .700 
VI SO : .000004 1 IDO 
Z : .6700 
U A : 50 .OOG 
UB :tJA 

E L PH A : 5 .82 7440-6 
RE TA: 1 . 683480-5 

V T SE : .000000006 3 1 00 

V A : 0 • 

V 0 : VA 

p: 20 .non 

HO:-. HO 00364600 

c 

C NODAL STRUCTURE INPUT 

C 

Zn n 3.0 00 
XN : 4 .000 
NX 167 
7 C:S.nDP 
NY : 4 

YMI 1 .600 
NX YINX 4 NY 
NX 1 : NX - 1 
NY 1 1 2 • N Y 
N R r N X 

Rr APIS ,991 ( (P?( I , J • 1 ).I:i,NX),j:l,NYll 

99 f ORMAT ( 8010.6 I 
00 98 j: 1 t NY 1 
00 9 1 I : 1 ,N X 
H z I ♦ ( J - 1 1 4NX 
97 PR 1 NJ-P 2 U ,J, 1 > 

Op CONTINUE 
975 CONTINUE 
X I P : X I A 


85 


p j* 
p?* 
03* 
R4* 
85* 
06* 
P7* 

00 * 

09* 
90* 
91* 
9?* 
93* 
94* 
9 S* 
96* 

07* 

90* 

99* 
100 * 
101 * 
1C?* 
10 3* 

1 C 4 * 
105* 
106* 

I C 7* 
10R* 
109* 
110* 
111 * 
11 ?* 
113* 
114* 
116* 

1 16* 

1 1 7* 

1 1 0* 
119* 
l?n* 
i?i* 
I??* 

1 ? 3* 

1 ? 4 * 
1?5* 

1 ? 6 * 

1 ? 7* 
1?P* 

1 ?o* 

1 30* 
131* 

I 3?* 

1 3 3* 

1 34* 
135* 
136* 
137* 


1 IP* 
139* 
140* 
14 1* 
147* 
143* 
144* 
146* 
146* 
14 7* 

1 4 P * 
149* 
16C* 
1* 1* 
1* ?• 
1* 3* 

) * 4* 

1 ' 5* 
l« 6* 
\ l 7* 
1* P* 

1 * 'J o 

1# n* 

I 6 |« 

1 6 ?* 
16 3* 
16 4* 


6 R -G A 

JPHOVrO 

JRMQ7I0 

JRHrQ 

MA XS?r 1 

Df l?;. 000500 

Df L 3- . ] Dl) 

Mil 

HNjr? 

C 

c CURVATURE SUM ANO DIFFERENT 

RX = I 1 . /RAX 1 ♦ ( 1 ./Rflx > 

PY: < 1 • /RA V 1 ♦ I 1 . /PB Y » 

RHOrRx *P Y 
GAMMArf PX-RV J /PHO 
C 

c ELLIPTIC INTEGRALS E AND 1 


XK-SORT ( ? • or 0 1 

** is THf flLIPTKITY PARAMETER 

1 AKrSQPT * 1 .-( 1 ./ ( XK**? I ») 

CALL CCl IF .£ «AK V T£P) 

XJKrsCRT((2.*r-r*( 1 .♦GAMMA >)/([*<! .-GAMMA) 1 I 
A 9 z X K - X JK 
XK 2X JK 

15 < ARS ( A9) .L T .DEL 1 GO TO ? 

GO TO 1 

StMlMAJOR ANOSfHlMJNOR AXLS OF CONTACT ELLIPSE 


C 

f 

c 


? S? = 1 1 .-XI A**? I/GA* I 1 ,-XIB**? )/GP 

Ti-( 3 «*(XK**ri*r*p)/(p]*pH 0 i 

A 1 : HI *S? )**i I ,/3. 1 
P 1 :A l /XK 


2 1 =-5*70 
??-.5*?C/XK 
? 3z?] **? 

?4r?^*i»p 

A A i r$ ? / p i 


T 3r,‘, * rx*(B1**?| 

T? = ,S*RV* I A 1*0? 1 
WP1T1 <6,31 RHO .GAMMA t XK *L f F 
3 FORMAT I4MPH0- ,016*S,5X .6HGAMMA - 
1 ?Hf = f 0l6.5 t 5X f 2HF: t Dl6.S) 


.016.5.5x f ?HKi t 0l6.5.5X, 


OlMf NSlONLfSS PAPAHF TFP GROUPING 


C 

C 

C 


EP=?./S? 

SMUT. S* f UA^UB » 

S m V - • 5 * 1 V A « VP ) 

V - SORT (Smli**^«smv**? 1 
THf TATA T AN ( SMV/SMUI 
Or l p/p i VAS 
U - V I S 0 *R X * V / f P 


‘ • r / j K 

Wl ?AR:p*px/ < ? ,*f P * A 1 1 
APCHAP; I • / ( l««(?.*RY)/( 3 « *R X | | 

<MlNAT? # 04*(ARCHAP*G*UI**.74/iyL7***0741 
DOy^ 0 ‘ji? .6 5*G *♦ . S 4 *u**. 7 /yl 7 AP **. 1 3 
MM INCH IG**.54I*IU**.651/«UL?**.073) 
R p IUIfc,l901 HMJNCL 
19C rORMATUH .PHMMINCAL: .016.6 » 

WPlTE<6.irM 00 u SON 



WPITl <6.41 Al.Rl t Pl # fp 


INITIAL PRr^SURf AND FILM THJCKNf 5 S 

7R~1./(U«|?.*PY)/I3.*RXM 
up zp« ( pmo**? » n r 
Ml <6.61 UP ,M0 

6 FORMAT UN ypt % D16. 5, 5X.3HHOI, 016.51 
GUVISf/viSO 
Q ? - f P / I 9600 . S?6« 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE JLS POOR 
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its* 

166* 

167* 

ice* 

1*9* 
170* 
171* 
17?* 
17 3* 

1 74* 
17S* 
176* 
177* 

1 7 fl* 
179* 

ian* 
ip i* 

ip?* 

I R 3* 

1 P4* 
IPS* 

1 P 6* 

1 P 7* 

1 PR* 
IRQ* 

J 90* 
19 1* 
19?* 

1 9 3* 

I 9 4* 
19 5* 
19** 

19 7* 

J 9R * 

1 9 9* 

?nn* 
? ri * 
?T?* 
?0 3* 
? r m * 
?r« • 
?r** 

?V 7* 
?^P* 
?r>9< 

?io* 
?i l* 
? i ?* 
71 3* 
714* 
715* 
716* 
7 1 7* 
?18* 
7 1 9* 
??n* 
?? 1 * 
???* 
??3* 
?? 4 « 

??s* 

77 6 * 
??7* 

??«* 

??9* 
7 3P * 

? 3 1 * 
? 3 ?• 
7 3 3* 
734* 
7 3 S* 
7 36* 
7 3 I* 
?3fl* 
7 39* 
740* 
?*» 1* 
?47* 
74 3* 
?4 4* 
74 S* 
746* 
74 ?• 
74 B • 


QSrf IPHA*EP 
06 rRt T A *[ P 
14 CONTINUE 

CALI SUP 6 I HO * M t MN 1 1 
IF IM.EQ.P) bO TO 333 

Y :0 
HNlrl 

4 P 1 TE C 6 ,SPll ; PHI MX 
S*n F 0 PMAT( 1 h f 6 HPHIMX: t D H .51 
333 CONTINUE 
C 

C VISCOSITY ANr DENSITY 

C 

DO 7 Nr I ,NX Y 

Df NS (N 1 r 1 .♦ <0 5 *PR (N >1 / U .♦L' 6 *PB <N) I 
VI S(N ) rQ j ** ( l . -( 1 . ♦ 0 ?*PP (N M **? ) 

XMU(N) rDENSt N 1 /VI S IN) 

7 CONTINUE 
C 

C RELAXATION COfFFJENTS A,B,C*D ( L*AND M 

C 

3 P 0 CONT 1 NUT 

SUHl rn .ODD 

0 0 9 J : ? , f Y 
DO ID I r 7 , N X 1 
N - I ♦ I J - 1 ) * N X 
N 1 r N ♦ 1 

N 7 r n - 1 
N 3 r N ♦ N X 

IF IJ.FC.NY) N 3 - N 
1,4 r N-N x 
V 1 r XHtM N 1 I 

Y ?r YMU I *7 ) 

Y 3 : VMU ( N l ) 

Y 4 : X M I ! I N 4 I 

Y 5 - M ( N 1 ) 

Y 6 : H ( N ? I 

Y Tz H ( N ) 

YPiHlM) 

Y 9 r H ( N 4 ) 

Y ID "Y 1 *SQR T C Y ' ) 
v 1 \z Y?*SC,RT I Yf ) 

Y 1 ?ty 3 * S 0 P T I Y P ) 

Y 15 -Y 4 *S 0 PTIY 9 ) 

A INU 7 3*1 ^ • * Y 1 4 Y 71 

H < NI = 74 * I Y 3 « 3 • * Y 4 I 
C I N 1 r 7 J * | Y 1 ♦ 3 . * Y 7 ) 

DL 71 N)r 74 *M.*Y 3 *Y 4 ) 

XI 1 = 4 .*<? 3 *(Y 1 *Y?I'* 74 *IY 3 4 Y 4 ) ) 

*17 = 1 . 5 / 1 Y 7 ** 5 . 5 ) 

XL 3 = 73 * IY 10 *( ?.*Y 5 - 4 .*Y 7 *V 6 ) • V 1 l*(YS- 4 . * V 7 ♦ 3 • * Y 6 ) ) 
XL 4 = 74 *IY 17 *I 3 .*Y»- 4 .*Y 7 *Y 9 I*V 13 *IY 8 - 4 .*Y 7 # 3 .*Y 9 M 
XL I N 1 r X t 14 XL?*IXL 3 *XL 4 ) 

XN 1 “ 1 ?.*U*B 1 *PX/I V 7 ** 1 . 5 ) 

X«?= I Y 5 *Uf NS 1 N 1 )-Y 6 *DENS(N?))*C 0 SI T Ht T A ) 

XH ir ( YR^Df NS I M 3 ) -Y 9 * 0 FNS IN 4 I ) *SIN( 7 HL 1 A ) 

XM(N):XMl*(71*XM?«7?*XM31 

10 CONT I Nur 
9 continue 
3 n 7 sumi m .non 

HSUH9:p 

C 

C PTLAX at ION TOPHU; A 

c 

DO 3 DR j:?,NY 
00 309 1 = 7 • NX ) 

M N = I ♦ I J - 1 I *N X 
HNA IHN ♦ 1 
MNP =MN~ 1 
MNC ’MN 4 N X 

I r I J.r 0 .NY ) MNC r H N 
MNP : MN -N X 

?PR<MN)tPH] (HN)-0RE*(PHI tMN)*IXM|MN)-A|MN)*PMl |MNA)-BtMN)*PMj I MND ) - 

1 -C l*N )* PH I (MNP )-Ul 7 IHN loPMUNNC * 1 / » L l*N I 1 
jf (?rR (MNI.L 1 .PHIHX I CO 10 55 P 

MSUH 9 :*SU* 9* 1 
7 PR I HN » : PMl M* 

ssn c oNTiNtir 

If l?PRIHN).Lr.n.) 00 TO 34 9 
Y I ft = < ?PP (PN I - PHI l«N I I //PIMPN I 
SUM 1 :SllM HAPS C V I S I 

no to 3 a b 

34 9 ?PR I HN ) 0 .0 00 


7*9* 

750* 
751* 
757* 
753* 
75*1* 
?ss* 
75 b* 
?5 7# 
758* 
?S9* 

?bn* 

?H* 

?b?* 
?bi* 
7bA* 
265* 
766* 
?fc 7* 
768* 
769* 
770* 
771* 
777* 
773# 
778* 
775* 

7 7 b* 
777* 

7 7 fl * 
779* 
780* 
7«1* 
787* 
78 3* 
?«9* 

78 5* 

7 6b* 

7 8 7* 
?P 8* 
789* 

79 0* 
791* 
797* 
?9 3* 
?9*»* 
795* 
79b* 
797* 
?*>8* 
799* 
300* 
in* 
107* 
301* 
in*** 
305* 
Iflb* 

in?* 

sr-M 

10 9* 

nn* 
n i* 
117* 

31 1* 

31 *»* 

3 1 S* 
3)6* 
317* 
31** 
319* 
3?n* 
371* 
177* 

37 3* 
37** 
375* 
37b* 

37 7* 
3?** 
379* 

3 3T • 
331* 

3 37* 


38 8 PM I I HN 1 r7PP (HM 
309 CONTINUE 
3np CONTINUE 
3 1 D CONTINUE 

H AUPrHAUP ♦ ! 

IT ISUH1 .LT.DfL 3) GO TO 311 
GO TO 307 
C 

c VISCOSITY ANP DENSITY I Tt RATION 

31 1 $UH?rn.nnn 

WPITf (6,680) HAUP t HS(fH9 
9*0 rOPHATUH tSHHAUPl, 18 f SX f bMHSUM9r # IA| 

HAUPrn 

psu^in.oor 

00 315 J; 7 ,N y 
DO 31b I r ? , NX ] 

N = I ♦ ( J- I )*NX 

PP(N) = (PMTINI/(M(N)**1.S)«(P0SH-1 # ) *PP$V IN ) | /PQSN 

PPSV(N)*PP(N1 

PSUNrPSUM*PPIN) 

Of NSN - 1 •♦(05*PPfN) ) / i ).«Ot*PP(N) | 

VISN*0 !**<]•-( 1 •♦Q7*PPIN) )**7) 

XHUNrnENSN/VJSN 

Y99t(*MUN-XMU<N | l/XHUN 

SUH?r< UH?*ARS I Y99 ) 

OE NS I N ) -Of NSN 
V I SCN 1 rV I SN 
3 1 b X*U t N 1 1 XHUN 
315 CONTINUE 
HNlrl 

ie (HH9.ro. n go to **59 

C X L L SUFbCHO,H,MNI ) 

*< c 9 CONTINUE 

JE NNi Jf NN ♦! 

WP] TE (6 ,B10 ) < UM? 

810 FpPHATClH , 5MSUH? ) ,01 b. 5 1 
IE ISUH7.LT.Dft 1 I GO TO 13 
IE CHAXS7.GT.3PO) GO TO 95U 
MAX<.?:HAXS7* 1 
GO TO 3C0 
C 

C APLIEO NOPHAt LOAD 

c 

1 3 CONTINur 

M IXS?r 1 

IF<HH9.r0.1) ro TO **01 
WPITf <b,**Pl > Jf NN 
8P1 rOPPATCIh ,9HJENNIf f Pr f lb) 
jr NN; P 

0U9T3 .0 no 

DO 16 1 - 7 «N X 1 
OU 1-7.000 
0U3I0.000 
in; I/? 

! wnrl h«i h 

IE ( T6H.Nl .1 1 ou 1 ; 1 .OOP 
DO 17 j:?,NY 
Nr I • | j-i )«NX 
I 7 OU 3rQU 3«PP|N) 
lb OUH - 0U9 « OU J*OU 1 

PpAPr*» »*1P*A 1*R 1*0U«*/ f 3.*^C*20 I 
PCHTCAr? .*EP*A1*R 1 *P S UH / < ^ C *70 I 
WPITC (b.SOOl PPAP # PC Hf C 8 

Sno ^ OPHA T MH t 5HP6AP;,01b.5 t SX # ?HPCMICn: f Dlb.Sl 

c Nrw cfNTPAL riiH thickness 
c 

IE CJtNN.tr .51 X Jt NNi . 500 
IE { Jf NN .GT.5 1 X Jf NNi .600 
I E I JE NN .GT . 10 1 Y JT NN r .800 
IE ( Jf NN «G I. IS I X Jf NN r .900 
IT IIBMHCI .t T. 1 .0-8) HOr-HO 

IE t HO. I T.p.l HOBAP;HO*(XJLNN«< l.-XJENNI*(P/P6AP)«*?77| 
IE I HO. GT. 0.1 H0RAp:H0*tlJlNN«tl.>YjENN)*(PBAP/P|«*777) 

7 S UH 3 : CP-PR Ap >/p 
SUNfrapsi 7SUH31 

VPMf U.iM P , PR A P ,HO,HCbAP «SUHJ 

18 EOPmAT < 1 M .7MPr,01b.8 t SX,SMPBAPr # Dlb.8,5X t 3MHr:,OU.8 # 

I b a f b MHPB A P : v 0lb»8 *SX , ShSUH 3 . ,016.81 

IT CSUH3.tT.nfL7 1 GO TO *9 
HP: HOP A R 
191 CONTINUE 

DO ANP J: 1 ,NY 


3 3 3* 
3 34* 
1 »6* 
336* 

3 ?«* 

T 39 * 
34n* 
^4 j * 
34 ?? 
34 3* 
344* 
34*.* 
346* 
34 7* 

3 4 8 * 
34 4* 
3*0* 
Tr 1 * 
35?* 
353* 

3 64 * 
3*5* 
'^6* 

7 5 7* 
358* 
359* 
3fc 0* 
HI* 

3< ?♦ 

36 3* 
364* 
365* 
3/6* 

36 7* 
364* 
364* 

* 70 * 
371* 
37?* 

3 7 3* 
374* 
376* 

3 76* 
37?* 

3 7 8 * 

3 74* 

38n* 

3 ' 1 * 

3fc ? * 

38 3* 
384* 

36 ** 

36 6 * 

38 7* 

36 6* 

38 9* 
340* 

34 1 • 
34?* 

34 3* 

34 14* 

365* 

’ 66 * 

3 f 7 * 

T( >m* 

364 * 

urr* 

4 f I* 

4 0 ?• 

4 0 3* 

40 4* 

4 ( 6 * 

40 7.* 

417* 

4? 8* 

4 l''C* 

4 10* 

4 )| * 

4 1?* 

4 1 1* 

4 I 4* 

4 1 « • 

416* 


C 


no 44 1 7 : 1 ,NX 

N 3 : I * ( J-] )*My 
HlN 3 |:M 0 *PX*( 6 (N 3 | 4 W | N3n 
441 PHIlNM'-PP (N3) * |HfNnM 

44 n coMiNur 

no to 3 oo 

14 C ON T 1 NUF 

k ( jom. r q . i > re to 70 r 

J 6 H : 1 

CT L 1 : . 5 00 
GO To 7 6 1 

7P ■) i r *J 8 H 0 ?.rc.n r,o ?c 7 ^o 

J n H 0 7 ~ ] 

nr l l : • ? p o 

(> n TO 7 5 6 

7*>U iriJMMCjY.rO.il on TO 756 
J M H C Y : 1 

OTL 1 r . ion 
7 f 5 CONTlNuf 

HP 6 4 0 J : l t o y 
no 64 ? J 1 , N> 

N 7 ! « NX « < J- 1 ) 

P ? < T • J , 1 > : p p < «4 | 

64 ? P?{ I ,J, 7 ) :h!N J 

640 CONTp.-ur 

nn 694 j: 1 ,Nv 
no 696 T - I * N X 
J1 ‘NY* J 


69 6 
694 


6*4 


7P? 

7*1 

4 P ] 
4 H 4 

4 1 f j 


4H 3 
404 

4 n 5 

31 7 
316 


30 1 
3P ? 
30 3 
-JO 4 

J n 5 
30 6 


J ? - N Y « J -j 

p?t I , J1 , 1 i:p?f I f J ?f | , 

V?i 1 * J 1 *? »:P? I I ,j? ,?j 
CONTlNUf 

r opm a r i n io 1 s i P ? * 1 

HOrMOPI p 

no 73 1 J - ! v N Y 
no 70 ? r: 1 9 N x 
N3-7*(J-J|*Nx 
H IN 3 | -MC*Px*< MN 3 I *UIN 3 I ) 

PHI(S 5 |:PI>(V 3 l*(H(NU*«l.b) 

CONTI^jf 

go to jnn 
CONTINUF 
CONT INUf 

WPITf 16 . 409 ) IPfi|N| # Nr| # 7 bCl 
rOPMMtlH ,?HPP t / t 7 m 1 H , 10013 . 5 . / 1 1 
“ Pln 16,4101 1 PMI CM ,Nri , 7 P 01 
0 pH 8 t I I H ,/HPH ( / *70 1 1 M , 10013.5 /!J 
WPTTt ( 6 , 405 ) lMINI.N:i, 70 ul * 

W*MTM 6 , 4 P 41 * ^ 1 N ) ,N : j t 7 0l ) 
kOUi r 6,405 ) ISINI ,N: | , 7 CL> 

r o wm* r mm .lOOU.s/ll 

“ 1 , "* iln inrNMNi ,n:| 7 nni 
rt)BM*T(iH , 7 moi nwiv,/, 70 « ih ,»nou b />. 

rOPH AT( |„ , 9 mvI 5 C 05 ITv ./,70 1 J h 

up I M 1 6 , 30 1 ) I A IN I ,N: 1 f 7 U0 1 

W 9 I M I 6 , 30 > ) •Ml.I.CiilTOb? 

yP1I( 1 6,30 3 ! iriN) #N : 1 # 7 LC) 

«^U 16 , 3 P 4 . IOt ? I N I ,N: J , ;on I 
W p !n 16,305 ) « X t CN| ,N“, , 7 u 01 

UPITM 6 ,JP 61 CXMlNI,N:i, 700 l 

f OPh A T ( in , IMA ,✓, 701 1M 


>ipdi j.s./n 


6 OP*A T | Ih 

1 OPm*t | , H 

10 PMUIH 
r 0 P M * F 1 |H 
f n W “A T | in 


. IMP , / , 701 1 M 
•INC,/, 701IM 

• IMP,/, 701 1 M 

* 1 Ml , / , 7 1 ) I IH 

tlM M,/,70,1M 


• 1 LO 1 3. 5 , / ) | 

• 10 0 13 . 5 ,/)) 

• lilCt3.S,/l| 
tlUOl J. 5 ./I 1 
* 100 13 . 5 , /|| 
« lu 013 . 5 , /II 


6 * 1 

f ij 


69 4 
69 8 


nO 650 j j fNy 

no #51 7 - 1 ,NX 

5 I *NX* < j 1 | 
p ’l T , J , 1 I P P | »| | 
n ? I T , J , ? ) m ( n ) 
ros r iNur 
no 6 9 * J 1 , N Y 
no 699 I | , NX 
J 1 N Y ♦ j 
J ? N Y ♦ 1 - j 

P’lTfJltll: P? f | , J? , 1 j 
*’ ’ * ’• J i ..’» *■ . i i , j.>, 
n OK T i vuf 


I • ?> 


89 




4 ] 7* 
4 ] 64 
4 I 94 
4204 
4714 
4??4 
4?3* 

4?44 

4?S4 

4?* 4 
4 ? 74 
H ?64 
4294 
4304 
4 X ) 4 
4374 
4 3 34 
4 344 
4 354 
4 364 
4 3 74 
43P4 
4 794 
4404 
44 ]4 
4424 
44 74 
4444 
44 54 
4464 

44 74 
4444 
4494 
4 5 04 
45)« 
4524 
4534 
4544 
4554 
4 5 64 

45 74 
4 f 44 
4594 
47 04 

46 J 4 
4724 

47 74 
4744 
4 7 54 
4764 
46 74 
4744 
4694 
4704 
4 7 14 
4 17 4 
4 7 74 
4 744 
4 754 
4 774 
4 7 74 
474# 

4 794 
4 6 04 
44 14 
4024 
4* 74 
4644 
4 6*4 
4 664 
46 74 
4004 
4694 
4 90 4 
49 14 
4924 
4934 
4944 
4954 
4964 
49?t 
4944 
4994 
5004 
50|4 


PUNCH 65? t l I f P? f 1 , j,n ) ,!=1 V NY) ,J= | ,NVf I ,N = 1 t ?| 

65? FORMAT I 6010 .5 > 

C 

c flow rah anal vsis 
c 

on 776 J - 7 * NY 
DO 777 Ir?,Nxi 
N - I « I ,J - J )4NX 
N1 

N 2 z N - ] 

N 7 - N t N X 

If ( J . F Q . N Y I N3=N 
N4 rN-NX 

PGPX: # 54?p4(PP INI) -PR CN? 1 1 
PORYi .54/C4 fPP<N3| -PR ( N 4 ) ) 

OYlrv4DrN51N|4H(Nl4C05ITHtTAl 

0Yltu4DFN5(Nl4H(N»45lN(THLTAl 

0K5:0rN5tN)4CH(Nl443l/u?.4RX4VISCNM 

OX ?r<jx 54PGRX/P1 

0Y?rgx54prRY/Al 

CX tcx 1 -ox? 

0 Y ~ 0 Y 1 -0 Y? 

QToTzOX 40Y 

07V-S0RT (0X44?«QY64?) 

P?(I*J*31:QTV 

0 A NC ZA T AMQV/QX 1 

UP nr 16 t 776 I N t Q T V ,OANG 

77fl rORNATCIH , ?HNz , H 0 * 5 X , 4HU 7 V : 1 0 1 3. 5 » 5* * 5H0 AN 6r t 01 3 *6 ) 

WRITt (6,7791 0X1 ,0X?,0Vl,yY2 

779 fOB«*TtlH .<.H0X|:,013.S,SX,i|H0X?:,DI 3. 5 , SX , 9 HOY 1 - , 01 3 .5 ,»X . 

I 9H0T2:,ni3.M 

UPIU (6,760) OTOT , QV,OX # PGRX ,P6RY 

7»n roBMXTiiM .5HCIOT:.ni3.5,SX,3HQY:,DlS.i,SX.3HOX: f DlJ.5.5x, 

1 bHPr,BX-,013.5.!.X,bHBCBYr t DI3.5l 
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™„,wsts„, H1 „ „o»u „„ 

MU1MIMIII AXtS Wilts Tilts, AS, , l,,V,„„, „„ , JIVlsloNS 


AM) CONDI I ION 1 OK TAPI. I- 3. 1 PHKVAII.S 


O •« , rdiiwlc»s| Pressure n 


1 la 

13a 

15a 

17a 

19a 

21a 

23a 

25a 

27a 


N i nr 


io’ 

0832 
. 0754 
. 0020 
0372 

0 


P«i 


0. 1250 
. 1200 
. 1093 
. 0899 
. 0539 

0 


29a 

31a 

33a 

355 

37a 


13 b 


41a 

43a 

45a 

47a' 

49a 

11a 


I 


"l 


n, t , r 

U- 

9>i inat inn, \\ 

“1 — — 

halm 

Total separation, S . w 

Patio 

s 

nil 

ill. 

H, ; " 

1 

cm 

in. 

R a 

1.0 

. 9? 

)0 

99 

0.0004 1 0”* 
. 0H2H 
. 05 59 
. 0055 
. 0521 
. 03 M 
. 03 10 

0.9340 10‘ 
. 0320 

• 0299 
. 02 5 o 
. 0205 
.0151 
. 0122 

99. 05 
19. 01 

0. 09 H 
3. OflH 

1 . 459 
7195 

. 4170 

0.0074 10’ 
. 0074 
. 0,374 
. OH74 
. 0079 
. 0914 
. 1052 

0. 0344- 10’ 3 
. 0344 
. 0344 
. 0344 
. 0340 
. 0300 
0414 

0. 0391 
. 0378 
. 0757 
. 1921 
. 7276 



.1*2- 2 

. 0103 

. 20 5 i. 

IM, 

1 

.0491 



i I 

. 0229 

. 01190 

j 1 804 

1494 

. 05HH 




. 0203 

. oorto 

• 1 2a2 

. 1 7 03 

. 0702 


. 999H 

.0103 

. 0072 

- 0944 

.2111 

. 083 1 


. 999H 

. 0105 

. 0005 

.0710 

. 2479 

. 0970 


. 9990 

.0152 

. 00150 

. 0550 

. 2003 

. 1 1 35 


. 9997 

. 0M0 

. 0055 

. 0440 

. 3325 

. 1309 
. 1 49 H 
. 1701 


. 9997 

. 0130 

. 005 I 

. 0355 

. 3 HO 5 


. 9990 

. 0122 

. 004 ,: 

. 0290 

. 432 1 


. 9990 

. 0114 

. 0045 

. 0240 

. 4 HO 9 

.1917 
. 214H 


. 9995 

.0107 

. 0042 

. 0201 

. 5450 


. 9995 

.0102 | 

• 0040 

. Ol 70 

. 00 HI 

. 2394 
. 2053 


. 9994 

. 0097 1 

. 003 H 

• 0145 

. 0739 
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1A, " K 3 - ■ aiAriAt TK,,,ST,L ' S OF ™> ; COMPONENTS OK FILM THICKNESS ALONG THK .SKMIMA.IOR AM, 
SKMIMINOR AXKS WHEN THESE AXES- ARK OIVIDKI) INTO TIIRKK EQUAL DIVISIONS 
AND CONDITION 2 OF TABLE 3. 1 PREVAILS 


Coordinates 

X Y 

7 b 7a 
9a 
11a 
13a 
15a 
17a 
19 a 
21a 
23a 
25a 
27a 

1 29a ' 

9b 7a 


0.3953 - 10 6 0. 5734 1 0 6 1. 
•3457 .5014 

•2144 .3109 


1. 844 10' 
1 . 039 
1.243 
. 7801 
. 5824 
. 407(1 
. 3917 
. 3376 
. 2967 
. 2047 
. 2390 
.2179 
1. 639 
1.243 
. 7861 
. 5824 
. 4676 
. 3917 
. 3317 
. 2967 
. 2049 
. 2393 
.2182 


mat it mi. w 

Rat i.. 

lotal separation, .s , w 

Halit. 

in. 

R w 

2 S 

cm 

in 

’ *3 

.720] io' 3 

35.33 

1. 897 10~ 3 

0.7407- 10" 3 

' 0.8130 

- 045 1 

0.275 

1 . 900 

. 7479 

1. 088 

. 4892 

1 . 830 

1 . 922 

. 7565 

3. 102 

.3095 

.6021 

2.092 

. 8235 

4. 897 

.2293 

.2720 

2.723 

1 . 072 

4. 170 

. 1841 
1542 

. 1468 
. 0882 

3.653 

1 . 438 

3.924 



4 . 83 i 

i . 902 


. 1329 

. 0572 

0. 238 

2.456 


. 1108 
1 042 

. 0392 

7. 809 

3. 098 



. M«_ «Mf 

i 8 . 1 1 8 

3. «20 


.0941 

. 0207 

11. 78 

4. 638 


, 0858 

.0158 

14. 00 

5. 535 



1.832 1.922 

■ 0027 2.092 

.2723 2.723 

■ 1409 3.653 


: — L! 1 - ---- I | ■ 01 50 | 14. W) j 5 .531 [. 

- CHARACTERISTICS OF THE COMPONENTS OK FILM THICKNESS AI.ONt. THE SKMIMA.IOR 
SKMIMINOR AXES WHEN THESE AXES ARE DIVIDED INTO KOKH EQUAL DIVISIONS 
AND CONDITION 2 OF TABI F 3 . 1 PREVAILS 


Kalin I K la si ic dr farina tin 


X V N ni/ 


0.4004-10 0.5807 10 

•3730 .5419 


1 . 803 10 ' 
I. 746 
1.514 
1. 177 
■ 8143 
. 0350 
. 5258 
. 4503 
. 3942 
. 3510 
. 3105 
. 2883 
. 2647 
. 2440 
. 2270 
.2120 
1.740 
1.514 
1. 177 
. 8143 
. 0360 
. 525 H 
. 4503 
.3942 
.3510 
. 3 1 05 
- 2883 
. 2047 
.2440 


0.7333 10" 
. 0874 
. 5902 
. 4033 
. 3200 
. 2500 
.2070 
. 1773 
. 1552 
. 1382 
. 1246 
. 1135 
. 1042 
. 01103 
. 0890 
. 0837 
. 0874 
. 5902 
. 4634 
. 3207 
. 2500 
. 2071 
. 1773 
. 1553 
. 1382 
. 1240 
. 1 135 
. 1042 
. 0903 


03. 43 
1 1 . 89 
3 . 905 
1 . 603 
. 6702 
3 543 
. 2 1 00 
. 13 50 
. 092o 
. 0000 
. 0487 
.0370 
. 0288 
. 0228 
. 0184 
. 0150 
11.90 
3. 908 
1 . 004 
. 0708 
. 3547 


FR EQUAL DIVISIONS 
ILS 

’ j Total separation. S . w 

Katin 

cm 

in. 

H 3 

1. 892 10"' 

0. 7448 10‘ : 

0. 7391 

1 . 893 

. 7452 

. 8450 

1 . 896 

. 7405 

1. 180 

I. 911 

. 7525 

2. 748 

2 2.019 

. 7948 


1 2. 427 

. 9555 


6 3.023 

1 . 1 90 

— 

j 3. 769 

1 . 484 


4.053 

1 . 832 



) 5.009 

2. 232 


0.810 

2. 681 


7 8. 072 

3. 1 78 


9.459 

3. 724 


10.97 

4. 318 


12. 00 

4. 959 


14. 34 

5. 647 


l . 89 3 

. 7452 

. 8450 

1 . 890 

. 7464 

1. 178 

1.911 

. 7523 

2. 747 

2. 018 

. 7940 


2. 426 

. 9550 


3 . 020 

I 189 


3. 707 

1 483 


4.051 

1 831 


5.064 

2. 230 


0. 805 

2. 079 


8. or, 7 

3. 170 


9. 454 

3. 722 





in. no 

4. 315 


12.59 | 

4 955 


14.33 

5. 043 





w 
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I-ABLK 3. 7. - CHARACTERISTICS OF THE COMPONENTS OF FILM THICKNESS A I. UNO THE SEMIMA-IOR AM) 
SEMIMINOR AXES WHEN THESE AXES ARE DIVIDED INTO FIVE EQEAE DIVISIONS 
AND CONDITION 2 OF TABLE 3. 1 PREVAILS 


Coordinates I Pressure, p 


Y 

N cm^ 

psi 

k 1- 

pm 




a 


11a 

0. 4027 1 0 e 

* 0. 5840 10 6 

1.000 

1. H62- 10“ 

13a 

.3859 

. 5597 

. 9999 

1. 787 

15a 

.3499 

.5075 

. 9998 

1.G37 

17a 

.2876 

. 4172 

. 9996 

1 414 

19a 

. 1726 

. 2503 

. 9993 

1. 124 


Katicj Elastic deformation, w Ratio Total separation. S • w 

*i 4 ~ I I H„ * i ' 



10' J 99. 05 

i in 

0. 098 

3. 008 
1.459 
. 7195 
. 4170 
. 2 65 H 
. 1 804 
. 1282 
. 0944 
.0716 
. 0556 
. 0440 
. 0255 
. 0290 
.0240 
.0201 
. 0170 
.0145 
19. 02 
6. 702 
3. 010 

1 . 460 
. 7202 
. 4174 
. 2660 
. 1805 
. 1283 
. 0945 
.0717 
. 0556 
. 0441 
. 0355 
. 0290 
.0240 
. 0201 
. 0170 
.0145 


1 . 880 1 O' 
1. 00! 

1. 881 
1 . 884 
1 . 895 

i 1 . 972 
I 2.264 
I 2.690 
3.218 

3. 840 

4. 547 

5. 339 
6.213 

7. 165 

8. 197 

9. 307 
10. 49 
11.76 
13. 10 
14.52 

1 . 881 
1. 881 
1 . 884 
1 . 895 
1. 971 

2. 263 
2. 687 
3.216 

3. 838 
4.544 
5.337 
6.208 
7. 160 
H. 192 
9. 299 

in. 49 
11. 75 
13. 09 
14.51 


0. 7403 10' 
. 7404 
.7407 

. 7417 
.746! 

. 7764 
. 8915 
1.059 
1.267 

1. 512 

1 . 790 

2 . 102 

2. 446 
2. 821 
3.227 

3. 664 

4. 133 
4.629 

5. 157 
5.715 

. 7404 
. 7407 
. 7416 
. 7460 
. 7761 
. 8911 
1. 058 
1.266 
1.511 
1. 789 
2 . 101 

2. 444 
1 2. 819 

3. 225 
3.661 

4. 128 

4. 625 

5. 153 
5 . 711 



■am E E ,1. CHAHACT ERISTICS OF THF COM PONKNTS OF KII.M THICKNESS AI.ONO T 
SKMIMA.IOH AM. SKMIMINOU AXKS WMKN THESE AXES AUK DIVIDED INTO THREE 
l Ql A I DIVISIONS ANJ) CONDITION 3 OF TAHLK 3. 1 PREVAILS 


Kl.islic <ir|urni;H icn. u 


0. 0248 lO 1 ’ 0.0360 10*' ft. 0399 10' 


0. 0157 10" 
. 013:* 

. 0087 
. 0040 
. 0028 
.0022 
. 0018 
. 0015 
. 0018 
• 0012 
. 0011 
. 0010 
. 0140 
.0124 
. 0098 
. 0083 
. 0072 
. 0064 
. 0058 
. 0053 
. 0049 
. 0045 
■ 0042 


Total s t» j 1 a 

■at inn. s . y.' 

cm 

in. 

0.0409 10" 3 

0. 0161 10* 3 

. 0409 

. 0161 

.0414 

- 0163 

) . 0480 

. 0189 

. 069 1 

. 0272 

. 0983 

. 0387 

. 1339 

. 0527 

1760 

. 0693 

.2249 

. 0883 

.2789 

. 1098 

. 3396 

. 1337 

. 4067 

.1601 

. 0411 

.0162 

.0414 

0163 

.0434 

.0171 : 

.0511 

.0201 1 

. 0630 

.0248 J 

.0782 

. 0308 

.0970 

.0382 

.1189 

. 046 H 

. 1438 

.0566 

. 1720 
j .2029 

.0677 

.0799 


TAIII.E 3.11. - CHARACTERISTIC!: (IF THE COMI-ONENTS Ol- HI M THICKNESS All.NI, THE 
SKMIMA.IOH AM) SKMIMINOU AXKS WHEN THESE AXES ARE DIVIDED INTO rol l; 
KQl'AL DIVISIONS AND CONDITION 3 OK IAHI.K 3 1 IMDAAII.S 


I’ri'ssiirc, p 
N ■nC mm 


f- Li.sl ir dt’l' -i ni.i i j- 1 ; i 


I 0. 0252 in* 1 0.03(15 in (i 0 04()1 1 o' 3 n. n| :, K jo :< [r;i 53 


0.040!) in' 
. 040!) 

040!) 

0411 
. U45 5 
ori'io 
nTiio 
Io;?i 
I2nh 
1 3:)3 
1 030 


ni 4.‘. .230(1 

.0113 2715 



.00581 4153 



7 


T 




\ 





TABLE 3 “ ' CHARACTKB ‘ ST ' C8 OF THK fOMTONENra OF FILM THICKNESS ALONG THE SEMIMAJOR AM) 
SEMIMINOR AXES WHEN THESE AXES ARE DIVIDED INTO THREE EQUAL DIVISIONS 


AND CONDITION 4 OF TAHI.E 3 . 1 PREVAILS 



TABLE 3. 12. - CHARACTERISTICS OF THE COMPONENTS OF FILM THICKNESS ALONG THE SEMIMA.JOR AND 
SEMIMINOR AXES WHEN THESE AXES ARE DIVIDED INTO FOUR EQUAL DIVISIONS 


AND CONDITION 4 OE TABLE 3. 1 PREVAILS 


iordi nates 


Pressure, p 


N cm 


0. ! IBB 10 


psi 

0. 1694 10 


1 la 
1 3 it 
15a 


. 1090 
.0914 
. 0555 


15H1 
. 1320 
. 0H05 


17U 0 0 

19:1 I I 

21a ; 


23a 

25a 


27a 

29a 

31a 

33a 

35a 


:n., 

1 39a 

lh 9a 
Hi I 


>1) 


ri> 


o. | 

lh i 


1090 
. 09 1 4 
. 0555 


'»* I 


i) ! 



* 


. 15H1 
. 1320 
. 0805 

0 



Elastic def 

nrmat inn, \ V 

Ratio 

T >ta! scpai 

■ill lull, S • VV 

Hath* 


S 

rm 

in. 

H., A 
*■ S 

cm 

III. 

l! :i 

1 . c 

00 

0. 8004 • 10' 3 

0.3411 10' 3 

03. 53 

0, 8801 10' 3 

0.3405 10' 3 

0. 8604 



.7922 

. 3119 

9. 039 

. 8744 

. 3404 

. 9549 



.0444 

. 2537 

2. 735 

8 80 ’ 

. 3 405 

1. 290 



.4310 

. 1097 

. 9417 

. 8849 

. 3 448 

3. 869 



. 2286 

. 0900 

. 3032 

. 9 82 2 

. 3 807 




. 1001 

. 0054 

. 1478 

1 . 290 

. 4077 




. 1339 

. 0527 

. 0854 

I. 702 

. 0,099 


. 9999 

. 1130 

. 0445 

. 0542 

2. 144 

. 8057 




. 0983 

. 03 H7 

. 0307 

2. 770 

1 . 09.3 




. 0809 

0342 

. 0200 

.3. 4.32 

1-351 




. 0780 

. 0307 

. 0191 

4 D-l 

1 . 0.38 


. 999K 

. 0709 

. 0279 

. 0145 

4. 40k 

1 . 990 


. 999H 

. 0050 

. 0250 

. 01 12 

5. 850 

2.303 


. 9998 

. 0599 

. 0230 

0084 

0. 807 

2. Oho 


. 9997 

. 0550 

. 0219 

. 0071 

7, 838 

’ 3. OHO 


. 9997 

. 05 1 8 

. 0204 

. 005 8 

8. 440 

3.522 


1 . 000 

■ 8319 

. 3275 

17 00 

8800 

. 3407 

. 5433 

1 . 000 

. 7028 

. 304.3 

0. 410 

8 810 

. 3471 

1 . 1 62 

. 9999 

. 0012 

. 2003 

2 948 

. 8854 

.3 48*, 

1 855 

. 9999 

. 5470 

. 2150 

t :»<n 

- 4124 

3542 


. 9998 

4780 

. 1882 

. 8844 

. 014 

4010 


. 9997 

, 4280 

1*85 

. 5700 

. 174 

. 4042 


. 9990 

. 3880 

. 15.3 0 

.3 844 

. ,3 81 

. 4455 


9995 

. 351,4 

• J 4<);i j 

. 27rt4 

. *214 

. 0433 


. 999-1 : 

. 3292 

1290 | 

2005 

423 

. 757* 


4993 

3 05 8 

1204 

'-4 

2 5 1 

881.1 


.9991 : 

. 2d 5 5 

t 12 4 


*■10 

1 *1.30 


. 9989 1 

2077 

. 10:5 4 

097.3 ! 

020 

1 184 


94 t,n ' 

. 25 4) 

0442 

07 8", | 

4* 2 

1 3*. 3 



•9. ! 

2 3 h 1 1 

‘>'*.<7 

. 008 1 I : 

3 440 

1. ...51 


9984 j 

. 2253 1 

0»h7 

0533 1 - 

4 45.3 

1 . 75.3 



1 


15 33 4 45.3 
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TABLE 5. 1. - EFFECT OF ELLIPTICITY PARAMETER ON 


MINIMUM FILM THICKNESS 


Elliptic it'y 
parameter, 
k 

Minimum film thickness 

Difference between 

"min “ d "min’ 
D l> 

percent 

Obtained from 
EHL point- 
contact theory, 

"min 

Obtained from 

least- square 

fit, 

8 min 

min 

1 

3. 367x10"® 

3.464x10"® 

+2. 88 

1.25 

4. 105 

4.031 

-1.80 

1.5 

4. 585 

4. 509 

-1. 22 

1. 75 

4.907 

4.913 

+. 11 

2 

5.255 

5. 252 

-.05 

2.5 

5.755 

5.781 

+. 45 

3 

6.091 

6. 156 

+1.08 

4 

6.636 

6.613 

-.34 

6 

6. 969 

6.961 

-. 12 

8 

7.048 

7.050 

+. 02 


TAELE 5. 2. - EFFECT OF DIMENSIONLESS SPEED 
PARAMETER ON MINIMUM FILM THICKNESS 


Dimensionless 
speed param- 
eter, 

U 

Minimum film thickness 

Difference between 
"min and "min’ 

D,. 

percent 

Obtained from 
EHL point - 
contact theory, 
"min 

Obtained from 
least -square 
hi, 

**min 

0.084 16x10' 11 

3. 926x10*® 

3.915*10"® 

-0.275 

. 1683 

6. 156 

6.252 

+ 1. 564 

.2525 

8. 372 

8.223 

-1.780 

.3367 

9. 995 

9.987 

-.078 

.4208 

11.61 

11.61 

-.004 

.5892 

14. 39 

14.57 

+1.280 

.8416 

18. 34 

18.54 

+ 1. 104 

1.263 

24.47 

24. 39 

-.320 

1.683 

29.75 

29.61 

-.467 

2. 104 

34. 58 

34.43 

-.432 

2. 525 

39.73 

38. 95 

-1.977 

2.946 

43.47 

43.22 

-.576 

3.367 

47.32 

47. 30 

-.042 

4.208 

54. 57 

I 54.99 

+.765 

5.050 

61.32 

62.20 

+ 1.430 


IT TABLE 5. 3. - EFFECT OF DIMENSIONLESS LOAD PARAMETER 

f ON MINIMUM FILM THICKNESS 

« ______ 


Dimensionless 
load param- 
eter, 

W 

Minimum film thickness 

Difference between 

"min *" d "min’ 
D l< 

percent 

Obtained from 
EHL point- 
contact theory, 

"min 

Obtained from 
least -square 
fit, 

"min 

o. uo6*ur 6 

6. 969x10"® 

6. 941x10 ® 

-0.41 

.3211 

6.492 

6. 599 

+ 1.65 

.3686 

6.317 

6.358 

+.64 

.5528 

6.268 

6. 172 

-1.52 

.7371 

6. 156 

6.044 

-1.81 

.9214 

6.085 

5.947 

-2.27 

1. 106 

5. 811 

5.868 

+.98 

1.290 

5.657 

5. 803 

♦2 58 



PEECTO1NG PAGE BLANK NOT FILLEU 





102 


TABLE 5. 4. * EFFECT OF SOLD MATERIAL AND LUBRICANT AS REPRESENTED IN DIMENSIONLESS 


MATERI/L PARAMETER ON MINIMUM FILM THICKNESS 


Solid ma- 

Lubricant 

Dtmenalonleaa 
material pa- 
rameter, 

G 

I Hme nil cm leas 
ipeed param- 
eter, 

U 

DlmenuionleM 
load param- 
eter, 

W 

Minimum film thickneaa 

Difference between 

H mln ““ “mln« 
D l< 

pcrctnt 



Obtained from 
EHL point- 
contact theory, 

»mln 

Obtained from 
least-aquare 
Jit, 

**mln 

Bronze 

Wu«iLitic 

2310 

0.3296M0" 11 

0.7216x10'® 

6.931x10'® 

6.073x10’® 

-0. 84 

Bronze 

Naphthenic 

3391 

.9422 

.7216 

17. 19 

17.404 

+1. 25 

Steel 

Paraffinic 

4322 

. 1683 

. 3606 

6. 317 

6.336 

+.31 

Silicon 

nitride 

Paraffinic 

6705 

. 112? 

.2456 

6.080 

6.038 

-.70 


TABLE 5. 5. - DATA SHOWING EFFECT OF ELLIPTICITY, LOAD, SPEED, AND MATERIAL ON MINIMUM FILM THICKNESS 
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TABLE 5.9. . EFFECT OF ELLIPT1CITY PARAMETER 


ON CENTRAL FILM THICKNESS 


Elllpticlty 

parameter 

k 

Central film thickness 

Difference between 

Obtained from 
EHL point - 
contact theory, 

H c 

Obtained from 
least-square 
fit, 

fi c 

"c “ d H c . 
°2> 

percent 

1 

6. 860xl(T 6 

6. 13P>T0‘ 6 

*10. 52 

1.25 

6.964 

6.565 

-5. 73 

1. 5 

7.001 

6.921 

*1. 14 

1.75 

7.015 

7.218 

+2. 80 

2 

7.402 

7.465 

+. 85 

2. 5 

7.653 

7. 841 

+2. 46 

3 

7. 845 

8. 102 

♦3. 28 

4 

8.292 

8.408 

+ 1. 40 

1 6 

8.657 

8.625 

37 

I 8 

8. 672 

8.675 

+. 03 


TABLE 5. 7. - EFFECT OF DIMENSIONLESS SPEED 


PARAMETER ON CENTRAL FILM THICKNESS 


Dimensionless 
speed param- 
eter, 

U 

Central film thickness 

Difference between 

Obtained from 
EHL point - 
contact theory, 

H c 

Obtained from 
least -square 
fit, 

»c 

H c and H 
» 2 , 

percent 

0. 08416*10" 

4. 917-10" 6 

4. 720-10" 6 

-4. 00 

. 1683 

7.517 

7. 495 

29 

.2525 

9. 999 

9. 825 

*1. 74 

. 3367 

11.40 

11.91 

>4. 47 

.4208 

13.07 

I 13. 81 

+5.66 

.5892 

17. 13 

17. 29 

+. 93 

. 8416 

21.35 

21. 94 

+2. 76 

1.263 

2j. 62 

28. 76 

-2. 90 

1.683 

35.50 

34. 83 

-1. 89 

2. 104 

41.05 

40. 43 

-I. 51 

2. 525 

46.64 

45.66 

-2. 10 

2.946 

51.08 

50,61 

92 

3.367 

55.56 

55. 33 

41 

4.208 

63.81 

64.20 

+.61 

j 5.050 

71.25 

72. 51 

+ 1. 77 


TABLE 3.8, EFFECT OF DIMENSIONLESS LOAD PARAMETERS 


ON CENTRAL FILM THICKNESS 


Dimensionless 

r 

Central film thickness 

load param- 

i 


eter, 

Obtained fro it 

! Obtained (n, ii 

W 

EHL point 

1 least -square 


j contact theory, 

| fit, 


H 

c 1 

H 

c 




o. looti to'k 

, H.6,)7 I0* t; ! 

8, 424 10 

2211 

7. 796 

7. 990 

. 3686 

7 . >05 

7. 720 

. 5528 

7 . 309 1 

7. >12 

.7371 i 

7. >17 

7, 368 

9214 

7 611 ^ 

7, 258 

1. 106 j 

7 416 j 

7. 170 

1.290 | 

6 762 

7 096 


Difference between ] 
1!,. and II , | 

0 2 . ‘ j 

percent 


2 H;* 
> 2 . 49 

>2. m 

*2. 78 
1 98 
-4. 63 
J. 32 
•4 94 
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% 



TABLE 5. 9. - EFFECT OF SOLID MATERIAL AND LUBRICANT AS REPRESENTED IN DIMENSIONLESS 


MATERIAL PARAMETER ON MINIMUM FILM THICKNESS 


Solid ma- 
terial 

Lubricant 

Dimensionless 
material pa- 
rameter, 

G 

Dimensionless 
speed param- 
eter, 

U 

Dimensionless 
load param- 

W 

Central film '^ickness 

Difference between 

H and H , 
c c 

D 2- 

percent 

Obtained from 
EHL point- 
contact theory, 

H c 

Obtained from 
least square 
fit, 

»c 

Bronze 

Bronze 

Steel 

Silicon 

nitride 

Paraffinic 

Naphthenic 

Paraffinic 

Paraffinic 

2310 

3591 

4522 

6785 

0. 3296*10" U 
.9422 
. 1683 
. 1122 

0. 7216 * 10" 6 
.7216 
. 3686 
. 2456 

8.422 t0“ 6 
20.70 
7. 505 
7.825 

8. 226*10"® 
20. 99 
7.819 
7. 585 

-2.33 
+ 1.40 
+4. 18 
-3.07 


TABLE 5. 10. - DATA SHOWING EFFECT OF ELLIPTICITY, LOAD. SPEED. AND MATERIAL ON CENTRAL FILM THICKNESS 



| 

■j i 

* a 


4 




E 


TABU: 6.1.- EFFECT OF STARVATION ON MINIMUM 
FILM THICKNESS FOR THREE GROUPS OF THE 
DIMENSIONLESS SPEED AND LOAD PARAMETERS 


Dimension- 

Group 

less inlet 
distance, 

1 2 

3 

m 

Dimensionless load parameter, \\ j 


0.36H610" 6 0,7371* 10*** 

0.7371 ■ K)- <5 




Dimensionless speed parameter, u | 


0.16H;M0' U 1.6H3-10" 11 

5.050- 10' 1 1 


Minimum film thicknes* 

C ft 

mm 

6 

j 

01.32 ID"* 5 

4 

** . 317* 10 20.27 

57 . 50 

3 

« - 26 1 27. H4 

51.70 

2.5 


It* . H9 

2 

5.097 23.40 

39.91 

1.75 


34.01 

1 . 5 

5.2.10 

27. >10 [ 

1.25 

3.945 




TABLE 0.2 . EFFECT OK DIMENSIONLESS SPEED AND LOAD PARAMETERS ON 
DIMENSIONLESS INLET DISTANCE AT FULLY FI.OODH) - STARVED BOUNDARY 

“P | Dimrnsi< ’"- Dimenninn- rD,„, t .n.si„n-| f ully Bonded I Lully (loo, le, I [ I), mens,,,,, 
lews speed lea., load leasrud.ua « «-..U-ul film nun,,,,,,,,, lea. ,,„e, 
P “ ri *7“*‘ r ' parameter, parameter, thuknrsa 1 , 1 ,,, ,|, sl4 „ a , 

1 W “* b 'V.K »*■«*. fully (Raided 

"nun, K starved 

boundary , 


0. InH.'i- 10"* 0. itiWi • 10' 


I 7 . 4M0 I (>’ 
| :»:i . f>j 
7o.o7 


0.221 10 ' 
20 . 20 
oi). yj 


TABLE 8.3. - EFFECT OF DIMENSIONLESS INLET -'STANCE UN DIMENSIONLESS 
CENTRAL FILM THICKNESS RATIOS 

upToinienaion- I Rat,,, „| central I Kut.o „l I Inlr, boundary 

lee. inlet film thicknesses film thu-knesses parunu-ter , parameter.,! 
distance, (or starved and l„ r starved and („, . |, . „ Wedeven, et al 

m flixKied condi- flooded comlUiuns, ny7 j ) 


H . « H r 

c , S C , t 


**n:in,S H mui, F 


l m - 1) (m w . l 














<a) Distortion of the element. 



(b) Local elastic deformation. 
Figure 3. 1. - Types of elastic deformation. 


I 








z 


ri 



Figure 3.2. * Elastic deformation of a semi-infinite body 
subjected to a uniform pressure over a rectanqular 
area. 







percentage difference in elastic deformation 









Figure 4 . 2 . Components ol film thickness for ellipsoidal solid near a pla 





Calculate k. a, b 
R 


| Calculate n and p by using eqs. 14.16) and 14.13) 

Calculate Aij.Bjj.Cjj.Dj j.Ljj, and My 
by using eqs. (4. 49) to (4. 55) 

Calculate ^j j n+ i by using eq. 44.59) 


/ l5 \ 

'pressure^ 

loop 

ranverg^ 


fun+1 n+l H i.j 


1721 Yes 


/ ls \ 

r \ r elaxa-\ 

tton loop 
converg- 
\ed? / 


SSU (4 63) ^ CalCUlate * byUSln ^ e<1 ' *3 


X ls \ 

normal - 
load loop 
converg- 
\ed? / 


Write results 


Figure 4 6. * Flow chart of main program. 


the first time \ Yes Calculate P, j initially 
by using eg. (4 57) 


Calculate n 
by using eq. (4 


Calculate H: : by using eq. (4 20) 


$. - p .u.3i7 

l.j.n r i,j ri i t j 


f jgure 4 7. Flow chart of subroutine SUB6, 










Dlmtflsbntoss 

prtsuirt, 

P - p IV 
A 0.7<W0' : 

B .66 

C .60 

D .50 

t .40 

F .30 

G .20 

H .10 



C 

B 



IV# VVMWUf 




'S* rj ? ”i£ n ,! l ° n ' M u * mur » ‘ f « l nlm •"* thr.rdim.nslon.1 r.pr.untl'ion o I pr.ssur. tor 

•lllptlclty ».m.tor *> of *. Th. M mmMhs p.r.n.1*, U. «. *nd G «. h,ld wn.Un. * hm In mJE&M 




Dimensionless 

pressure, 

P • p/E a 

A 0.83bcl0’ 3 
B .75 

C .65 

D .55 

£ .45 

F .35 

G .25 

H .15 



/ / 


\ \ \ \ \ \ \ 
'-A. \ \ \ \ \ \ \ 

i ' \ , \ i i 


\ \ \ \ ’ / / / / / / / 

\ ' V / / / / / / / 



(a) Contour plot of dimensionless pressure. 

Figure 5. 3. - Contour plots of dimensionless pressure and film thickness and three-dimensional representation 
of pressure for eiliptlclty parameter (k) of 6. The dimensionless parameters u, W, and G are held constantas 
defined in equation (5. 10K 
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(c) Three-d (mens tonal representation of pressure. 
Figure 5.3. - Concluded. 
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Dimensionless 
pressure, 
P-p /£' 

A 1 . OOxlO - 3 

B .95 

C .90 

D .80 

E .70 

F .50 

G .30 

H .10 


H 



E 

D 


C 

ft 


A 




Contour plot of dimensionless pressure. 

FI ? Ur * L 4 * n 1 C « n ,l? ur pk),$ of dfmeni *onless pressure and film thickness and three-dimensional representation of ores- 
.Ton (5 101 P N P,ram * Ihedimenstonless «. W. and G ,re heM «. % 


Dimensionless 
Him thickness, 
H - h/R x 

A 6.8xl0‘ 6 

B 7,0 

C 7.3 

D 7.7 

E 8.2 

F 8.8 

G 9.6 

H U.O 


H 



(bl Contour plot of dimensionless film thickness. 
Figure 5.4. * Continued. 




Dimensionless 

pressure, 

P ■ p/E’ 

A 1.05x10’ 3 
B 1.00 

C .95 

D .90 

E .80 

F .70 

G .50 

H .20 



... « M .vui k'v i Vi p»e>»ure. 


F>9ur * ' Conlour PtoholdlrrunslonMss praswra and dim thickness and thraa-dimanslonal repre- 

*r i ,Mlp,IC !lL P * ram * ter *’ of 1 ^ mramatarj u. W. endTare 

new constant as daflnad In aquation (5. 10). 




Dimensionless 
pressure, 
P-P IV 

A 1.10(10* 3 

B 1.05 

C 1.00 

D .90 

i .80 

f .70 

G .50 

H .20 



8 A 


A 



^ leisure. 

^nlitton of pre'sju re* tor 1 **! fiptic f»y 5 pa ^ a mVf er^k t of * ^ ai'nf.nXl*™' ,hr "~ dim,nslcnl1 

ir« held constant as defined In equation (5 101. ? ' 5 ' ^ dlm,niionl * s$ (•»«"*<»'» IJ. w. and c 
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I. J'H» jf'«) 






Dimensionless 
film thickness. 
H-h/R x 

A 5. 4x10“ 6 
B 5.7 

C 6.0 

D 6.5 

E 7.0 

F 8.0 

G 9.0 

H 10.0 



(bl Contour plot of dimensionless film thickness. 
Figure 5.7. - Continued. 







Dimensionless 

pressure, 

P - p/E* 

A 1.4010* 3 
B 1.35 

C 1.30 

D 1.20 

E UO 

F .90 

G .60 

H .20 



(a) Contour plot of dimensionless pressure. 

Figure 5. 9. - Contour plots of dimensionless pressure and film thickness and three-dimensional repre- 
sentation of pressure for ellipticity parameter (k) of 1. 5. The dimensionless parameters U, W, and G 
are held constant as defined in equation (5. 10). 



X 


Jk 



Dimensionless 

pressure, 

P-p/E' 

A 1.7x10’ 3 

B 1.6 

C 1.5 

0 1.4 

E 1.2 

F 1.0 

G .7 

H .3 



U) Contour plot o( dimensionless pressure. 

1 ° - ; Co " l<HJI ' PWJofdlmtnslontej prnsura and film tfikknajs and ttiracdlmansional 

W a^ctShlm’T, bralllpflc^p.ramifarfltlofl.ZS. Th« dirrwnnonless paramatars u. 
w. and G are held constant as deflnad in muation (5. 10). 


Dimensionless 
mm thickness, 
H - h/R x 

A 4. 3klO" 6 

B 4.6 

C 5.0 

D 5.5 

E 6.0 

F 6.6 

G 7.4 

H 1.2 



(b) Contour plot of dimensionless film thickness. 
Figure 5. 10. * Continued. 
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Dimension lew 
pressure, 
P-p/E’ 

* UWO ' 3 

B 1.6 

C 1.4 

D 1.3 

E 1.1 

F .9 

G .6 

H , 3 




it) Contajr plot of dimensionless pressure 


Figure 5. 12. - Contour plots of dimensionless pressure and 
sure tor dimensionless speed parameter (U> of 5. 05(kl0' " 
stint as defined in oration G. 111. 


film thickness and three-dimensional representation of pres 
. The dimensionless parameters k. W, art G are heto con- 
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Dimensionless 

pressure, 

P-p/E' 

) fiYin'3 




(a) Contour plot of dimension's pressure. 


Figure 5. 1 j. - Contour plots of dimensionless pressure and film thickness and three 
pressure for dimensionless speed parameter OJ) of 4. 208x10““ The dimensionless 
held constant as defined in equation (5.18). 


thickness and three-dimensional representation of 
, The dimensionless parameters k, W, and a are 


£ 






Dimensionless 

pressure, 

P-p/E' 

A 1.8x10* 3 
B 1.6 

C 1.4 

D 1.3 

E 1.1 

F .9 

G .6 

H .3 



i: 


\ 

\w 



\ J />:!; : / 



(a) Contour plot of dimensionless pressure. 

Figure 5. 14 - Contour plots ot dimensionless pressure and film thickness and three-dimensional representation of 
pressure br dlmens ontoss speed parameter (U) of 3. 367X10 11 . The dimensionless parameters k. W. and C are 
held constant as defined In equation (5. 18). 
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rW 
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| 

i i 





Dimensionless 

pressure, 

P - p/E* 

A 1.8x10“ 3 
B 1.6 

C 1.4 

D 1.3 

£ 1.1 

F .9 

G .6 

H .3 



(a) Contour plot of dimensionless pressure. 

Figure 5. 15. - Contour plots of dimensionless pressure and film thickness and three-dimensional representation 
of pressure for dimensionless speed parameter <U) of 1 946xl0 -11 . The dimensionless parameters k. W. and G 
are held constant as defined in equation (5. IB). 
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Dimensionless 
pressure, 
P-p IV 

A 1,8x10' 3 
8 1.6 

C 1.4 

D 1.3 

E 1.1 

f .9 

G .6 

H .3 





A 



(a) Contour plot of dimensionless pressure. 

^pressure ter ^uTef^SSiiVo' ^ ^^l?dlm3OTSoflh^* ,,,en *2 ,, ?* rep, ' 5 * nt, " on 

•re held constant as defined in equation 15.18). ' Th ' dlmenslonl « s constant* k, W. endG 





Jrnmm 


t 


* 



Oi-nemlcnlws 
pressure, 
P'p IV 

* 1.8jc10‘ 3 

8 1.6 

C 1.4 

0 1.3 

f 1.1 

P .9 

G .6 

H .3 




U> Contour plot c/dlmom ton loss prossure. 


r^urt >. u. * contour plots erf dlmtnstonloss prossura and fl 
or P^Msgro tor dimension lass spaod paramo tor (Ui of 2. lotxj 
sr* hold constant as doflrvd In equation 15. ID. 


Ihkknass and thraa-dlmans tonal raprasantatton 
. Tho diman sion Iasi constants k. w, ant G 


Dim«nsionlrts 
film thicknus, 
H-h/R x 

A M.fclQ' 6 
B 5.4 

C 36.0 

D 37.0 

E 39.0 

F 42.0 

G 46.0 

H 51 0 


H 




Ibl Contour plot of dimen lion Itu film thlc*n«». 
Fiqurf 5.17. ■ Continue. 
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Dimensionless 

pressure, 

P-p/E' 

A 1.8x10' 3 
B 1.6 

C 1.4 

^ 1.3 

E 1.1 

f .9 

G .6 

H .3 



/ KJ ! i I 

\ / A ///// / / 



(al Contour plot of dimensionless pressure, 

teflon of pressurefer j ^< 83 !^* T t ' f !|^' d,m(>ns,onal "Pr«en- 
W and P. %ra hoiH rnnrfinl .. 1 _ * * * The dimensionless parameters k. 


IAi r .... u . 1™* '•IMBICI IUM 

W, and G are hold constant as defined In equation (5, 18). 


Dimensionless 
film thickness, 
H«h/R x 

A 30. (WO" 6 

B 30.5 

C 31.0 

0 32.0 

E 34.0 

F 37.0 

G 41.0 

H 46.0 



(b) Contour plot of dimensionless film thickness. 
Figure 5,18. - Continued. 




Dimensionless 

pressure, 

P-p/p 

A 1.8x10' 3 
B 1.6 

C 1.4 

D 1.3 

£ 1.1 

F .9 

G .6 

H .3 
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Dimensionless 
Him thickness 
H-h/R x 

A 24. 8x1 O’ 6 

B 25.3 

C 26.0 

0 27.5 

E 29.0 

F 32.0 

G 36.0 

H 41.0 



(b) Contour olot of dimensionless film thickness. 
Figure 5.19. - Continued. 
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Dimensionless 

pressure, 

P -p/E' 

A 1.34x10' 3 
B 1.45 

C 1.40 

D 1.30 

E 1.15 

F .90 

G .60 

H .30 



C 


A 



(a) Contour plot of dimensionless pressure. 


Figure 5.20. - Contour plots of dimensionless pressure and ^hkknessandtjy 
representation of pressure for dimensionless speed parameter tU) o Ml WO ■ 
less parameters k, W. and G are held constant as defined by equation a 181. 


ee-dimensional 
The dimension- 


Dimensionless 
film thickness, 
H-h/R x 

A 18.5X10' 6 

B 19.0 

C 20.0 

0 21.0 

E 22.0 

F 24.0 

G 27.0 

H 32.0 


H 



(bi Contour plot of dimensionless film thickness. 
Figure 5.20. - Continued. 




174 


Dimensionless 

pressure, 

P - p /£' 

A l.Mxl O' 3 
B 1.48 

C 1.40 

0 1.30 

E 1.15 

F .90 

G .60 

H .10 


H 



v * 

B ■ 

P A 



r ‘"'"'muh pie»ure, 

!ure 'Ofdlm«niionl«s s^mI ,hrM ' dlrnenslonal feprewntttton of prw 

stint is defined in equation <5. 18). he d,m<ns,on «s parameters V, W, and G are held con- 







OlmensionlMs 

pressure. 

P-p/E' 

* I.44xI0'3 

8 1.40 

C 1.3? 

D 1.22 

E 1. 05 

E .80 

C .50 

H .10 




**•"*"• ter dim^ton^ n?»cknen and mrte-dimtniional 

4ff held wm tint **, dtfintd by •quition (5.1J). d } m«nsionl#4i paranufer j k. 


Dimensionless 
film thickness, 
H • h/R x 

A 11. 7xl0" 6 

B 12.0 

C 12.4 

0 110 

£ 14.0 

r i6.o 

G 19.0 

H 24.0 


H 



(b) Contour plot of dimensionless film thickness. 
Figure 5. 22. - Continued. 


(c) Three-dimensional representation of pressure. 
Figure 5 22, - Concluded. 


Dimensionless 

pressure, 

P-p/E' 

A 1.4*xl0' 3 
B 1.45 

C 1.39 

D 1.31 

E 1.20 

F 1.00 

G .70 

H .10 



(a) Contour plot of dimensionless pressure. 

Figure 5.23. - Contour plots of dimensionless pressure and film thickness and three-dimensional representation 
of pressure for dimensionless speed parameter OI> of 3. 367xl0’ lz . The dimensionless constants k, W, and G 
are held constant as defined In equation (5. 18). 


Dimensionless 
film thickness, 
H-h/R x 

A 10.1X10' 6 

B 10.3 

C 10.6 

D 11.0 

E 12.0 

F 14.0 

G 17.0 

H 22.0 


H 


G 



(b) Contour plot of dimensionless film thickness. 
Figure 5.23. - Continued. 



Dimensionless 

pressure, 

P-p/E* 

A 1.50(10' 3 
B 1.44 

C 1.34 

D 1.20 

E 1.00 

F .70 

G .20 


G 




B 


A 



(a) Contour plot of dlmens ton less pressure. 

Figure 5. 24. - Contour plots of dimensionless pressure and film 
resentatton of pressure for dimensionless speed (U) of 2.525x10 
k, W, and G are held constant as defined In equation (5.18). 


thickness and three-dimensional rep- 
1 . The dimensionless parameters 


Dime jsionless 
film thickness, 
H-h/R x 

A S.felO' 6 
B 9.0 

C 9.6 

D 10.6 

E 12.0 

F 14.0 

G 18.0 


fl- 



it)) Contour plot of dimensionless film thickness. 
Figure 5.24. - Continued. 




Dimensionless 

pressure, 

P.p/E' 

A 1.5*10'* 

B 1.45 

C 1.35 

D 1.20 

E 1.00 

F .70 

C .20 


G 



C 

r* 



Flgur* 5. 25. - Contour plots o ( dlmensk>nl*s> prtssurt tnd film thkkn.ss jnd thrM-dimtnslonil 
rjprmntitlon ofprwjurt for dlm«nsk>nl«ss spood pirinwttr lUtof |,tt*lo'“ Thodlmon- 
stontess parameter s k, W, and G are hefcl constant as defined in equation (5. 181. 
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Dimension l«s 
film thickness, 
H * h/R x 

6 . 25 X 10' 6 

6.40 

6.60 

7.00 

7.50 

8.70 

iaoo 



plot of dimensionless film thickness. 
: lguf • 5. 25. - Continued. 



Mi Three 4immio«Ml fomentation of pressure 
Figure V 2V - Concluded, 




















Dimensionless 

pressure, 

P-p/E’ 

A L 17xlO' 3 
B 1.14 

C 1.10 

D 1.00 

E .90 

F .70 

G .40 

H .10 


H 



(a! Contour plot of dimensionless pressure. 

Figure 5. 30. * Contour plots dt dimensionless pressure and film thickness and three-dimensional representation 
of pressure for dimensionless load parameter (W) o' 0. 3tt6xl0 . The dimension less parameters k. U, and G 
are held constant as defined in equation (5. 231. 








U* TTirw dt r«pre**ntotKm of prmuft 

FbureVU ConikjdaJ. 
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Dimensionless 

pressure, 

P-p/E' 

A 1.53(10" 3 
B 1.45 

C 1.35 

D 1.20 

E 1.00 

F .70 

G .20 



(a) Contour plot of dimensionless pressure. 

Figure 5. 32. Contour plots of dimensionless pressure and film thickness and three-dimensional 
representation of pressure for dimensionless load parameter (W) of 0. 7371xlO' 6 . The dimension- 
less parameters k, U, and G are held constant as defined by equation (5.23k 
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10 Ttiree-dtmens tonal representation of pressure 
Figure 5 . £. - Concluded. 


Dimensionless 

pressure, 

P-p/E' 

A 1.67x10" 3 
B 1.60 

C 1.50 

D 1.35 

E 1.10 

F .70 

G .20 


G 
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B 

A 



(a) Contour plot of dimensionless pres: ire. 

Figure 5. 33. - Contour plots of dimensionless pressure and film thickness and three-dimensional 
representation of pressure for dimensionless load parameter (W) of 0. 9214x10 . The dimension- 
less parameters k, U. and G are held constant as defined in equation (5. 23). 
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Dimensionless 

pressure, 

P - p/E’ 

A 1.75x10' 3 
B 1.70 

C 1.60 

0 1.40 

E 1.10 

F .70 

G .20 



.■J r> 
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(a) Contour plot of dimensionless pressure. 

Figure 5. 34. - Contour plots of dimensionless pressure end film thickness and three-dimensional 
representation of pressure for dimensionless load parameter (W) of 1. 106xl0' 6 . The dimension- 
less parameters k, U, and G are held constant as defined by equation (5. 23). 
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I. 


Dime ns ion less 
film thickness, 
H ■ h/R x 

A 6.1x10" 6 
B 6.4 

C 7.0 

D 7.7 

E 8.7 

F 10.0 

G 12.0 



(b) Contour plot of dimensionless film thickness. 
Figure 5. 34. * Continued. 
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(bl Dlmen»lonless film thkkness 

Figure 5. ». Vorletlon of dimensionless pressure end film thldmess on X-uls 
for three v«.ues of dimensionless loeri perimeter. The velue of Y Is held fixed 
n«r uUI ctn»tr of contact. 










Edges of computing arec, where pressure is assumed ambient 0 * 
Figure 6. 1 . - Computing area in and around the Hertzian contact circle. 



Inlet boundary parameter, tm - D/|m' 1) 

Figure 6. 2 . Influence of inlet boundary parameter upon central film thickness ratio 








220 


Dimensionless 
film thickness, 
H - h/R x 

A 6.5x10 
B 7.0 

C 7.5 

D 8.0 

E 8.5 

F 9.0 

G 10.0 

H 11.5 


... H 



(b) Dimensionless film thickness. 


Figure 6. 3. - Concluded. 




Dimensionless 
film thickness, 
H-h/R x 

A 6. 4x10" 6 
B 6.7 

C 7.2 

D 7.8 

E 8.5 

F 9.3 

G 10.2 

H 11.2 



Dimensionless 

pressure. 

P-p/E' 

A 1.17x10" 3 
B 1.14 

C UO 

D 1.00 

E .90 

F .70 

G .40 

H .10 


H 



(«) Dimensionless pressure. 

d 6 i, 5 inc^°m of 2 ^ dlm * nstonl *” n,m ' hlck "«* 



Dimensionless 
film thickness, 
H-h/R x 

A 6. 1x10' 6 
B 6,4 

C 6.8 

D 7.3 

E 7.9 

F 8.6 

G 9.4 

H 10.3 



(b) Dimensionless film thickness. 
Figure 6.5, - Concluded. 
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Dimensionless 
film thickness, 
H-h/R x 

B 5. 4x10' 6 

C 5.7 

D 6.1 

E 6.6 

F 7.2 

G 8.0 

H 9.0 



(b) Dimensionless film thickness. 
Figure 6.6. - Concluded. 



Dimensionless 
film thickness, 
H - h/R x 

A 4.&10’ 6 

B 4.2 

C 4.5 

0 4.9 

E 5.4 

F 6.0 

G 6.8 

H 7.8 



0 




(b) Dimension less film thickness. 
Figure 6. 7. - Concluded. 
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Dimension Itss 
pressure, 
P-p/P 

A l.toior 3 
6 1.6 

C 1.4 

D U 

E U 

F .9 

G .6 

H .3 



wimiiiuwiihh ensure, 


,nd dlm * ns,on, * 5s ,iim ,htc * n *” * 



Dlmansionless 
pressure, 
P-p IV 

A 1.8*10* 3 
B 1.5 

C 1.4 

D 1.3 

i 1.1 

F .9 

G .6 

H .3 



(at Dimensionless pressure, 

FlQurt 8. 9. * Contour ptota el dimensionless pressure anJ dimensionless film thickness tor dimensionless 
Inlet distance m el 4 and yeup 2 ol tobia k 2. 
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Dlm#nj|oftltss 
film »hicltn*«, 
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Ibl Dfnwnslontoi film thklrms. 


Flgurt 11 

- Conclude. 





Dimensionless 

pressure, 

P • pip 

A 1.8x10' 3 

B 1.6 

C 1.4 

0 1.3 

E 1.1 

E .9 

G .6 

H .3 



lal Dimensionless pressure. 


less Inlet dlstencem of^S ert'youp ^"rfttfle'sT' ,nd dlm * nslonl * SJ ,llm sickness lor dimension 



Dimensionless 
film thickness, 
H-h/R x 

A 26.5xl0' 6 

B 26.8 

C 27.4 

0 2B.4 

E 30.0 

F 32.0 

G 35.0 

H 39.0 


H 



(b) Dimensionless film thickness. 
Figure 6.11. - Concluded. 



Dimensionless 

pressure, 

P ■ p/E* 

A 1.8x10' 3 
B 1.6 

C 1.4 

D 1.3 

E U 

F .9 

G .6 

H .3 



(a) Dimensionless pressure. 

Figure 6. 12. - Contour plots of dimensionless pressure and dimensionless film thickness for dimen 
slonless inlet distance m of 2 and group 2 of table 6. 2. 
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Dimensionless 

pressure. 

P-p/E* 

A 1.8xl0‘ 3 
B 1.6 

C 1.4 

D 1.3 

E 1.1 

F .9 

G .6 

H .3 



(a) Dimensionless pressure. 

Figure 6.13. - Contour plots of dimensionless pressure and dimensionless film thickness tor dimensionless inlet dis- 
tance m of 6 and group 3 of table 6. 2. 


REPRODUCIBILITY OP Till* 

orkm.val r,\c,: is poor? 
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Dimensionless 

pressure, 

P - pfE' 

A 1.8x10' 3 

B 1.6 

C 1.4 

D 1.3 

E 1.1 

F .9 

G .6 

H .3 


/ 
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(a) Dimensionless pressure. 

Figure 6. 14. - Contour plots of dimensionless pressure and dimensionless film thickness tor dimensionless Inlet dis- 
tance m of 4 and group 3 of tele 6.?. 
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Dimensionless 

pressure, 

P • p/E' 

A 1.8X10' 3 

8 1.6 

C 1.4 

D 1.3 

F 1.1 

F .9 

G .6 

H .3 
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<e) Dimensionless pressure. 


Figure 6. 15. - Contour plots of dimensionless 
distance m of 3 and group 3 of table 6. 2. 


pressure and dimensionless film 


thickness for dimensionless 


Inlet 
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DlmtnsbnlMs 
film thkkntss, 
H-h/R, 

52.tolO‘ 6 
5?. 5 

53.0 

54.0 

56.0 

59.0 

63.0 

68.0 
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C 



4b) Dimtflsionitss film thickness. 
Floor# 6.15. - Concluded, 
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DImans tonltss 
prassura. 
P-p IV 

A 1. tolO* 3 
B 1.6 

C 1.4 

0 1.3 

1 U 

F .9 

G .6 

H .3 
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A 
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(6) Dlmanstonlw* prassura. 


Flgura 6. 16. - Contour plots of dim am ton lass prassura and 
InWdlstonc* m of 2. 5 and group 3 of tabla 6. 2. 


dlmanslonlass film 


thicknass tor dlmanslonlass 
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Otawnsbrifos 
film mtefcnws, 
H-h/R, 

A 47.CWO* 6 

B 47.$ 

C 48.0 

*) 49.0 

l 51.0 

F 54.0 

C 58.0 

H 610 


H 
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A 



(b> Olnwnslonim film mtdrms. 
Flgurt 4.16. * Coockxbd, 
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Dimensionless 

pressure, 

P-p/E’ 

A 1 . 8x10" ^ 

B 1.6 

C 1.4 

D 1.3 

E 1.1 

F .9 

G .6 

K .3 



(a) Dimensionless pressure. 

Figure 6. 17. - Contour plots of dimensionless pressure and dlmei, unless film thickness for dimen- 
sionless Inlet distance m of 2 and group 3 of table 6. 2. 







Dimensionless 

pressure, 

P - p/E 1 

A 1.8x10’ 3 
B 1.6 

C 1.4 

D 1.3 

E 1.1 

F .9 

G .6 

H .3 


H 



(a) Dimensionless pressure. 

Figure 6. 18. * Contour plots of dimensionless pressure and dimensionless film thickness for dl 
mensionless Inlet distance m of 1. 75 and group 3 of table 6. 2. 



Dimensionless 
film thickness, 

H - h/R x 

A 34. 8x10’ 6 

B 36.5 

C 36.5 

D 37.5 * 1 

E 39.0 

F 41.0 

G 45.0 

H 50.0 ‘3 



(b) Dimensionless film thickness. 
Figure 6.18. - Concluded. 
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Dimensionless 

pressure, 

P-p/E' 

A 1.8x10“ 3 
B 1.6 

C 1.4 

D 1.? 

E 1.1 

F .9 

G .6 

H .3 


H 



(a) Dimensionless pressure. 

Figure 6.19. - Contour plots of dimensionless pressure and dimensionless film thickness for 
dimensionless Inlet distance m of 1. 5 and group 3 o, table 6. 2. 











